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Abstract 

A  survey  is  given  of  the  theory  of 
equilibrium,  motion,  and  stability  of 
a  plasma  in  the  guiding  center  limit. 
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1 .   Introduction 

Perhaps  the  single  feature  which  dominates  the  landscape 
In  any  view  of  plasma  physics  Is  the  complexity  arising  from 
the  large  number  of  simultaneous  and  Interfering  phenomena 
[as  compared  to  the  relative  simplicity  of  the  related 
classical  disciplines  of  fluid  dynamics,  kinetic  theory, 
electromagnetic  theory].   We  cannot  expect  the  growth,  even 
ultimately,  of  a  significant  quantitative  theory  using  sub- 
stantial mathematical  tools  which  will  embrace  an  appreciable 
part  of  the  physics  and  at  the  same  time  Incorporate  a  non- 
trlvlal  geometry.   One  type  of  model  allows  the  use  of 
significant  linear  techniques  (essentially  Fourier  analysis 
enhanced  by  asymptotlcs)  to  describe  a  relatively  large 
amount  of  physics  by  remaining  within  the  realm  of  trivial 
geometries  such  as  an  Infinite  homogeneous  plasma,  possibly 
with  a  superposed  small  linear  macroscopic  gradient.   We 
shall  be  concerned  with  a  model  of  a  second  type  containing 
somewhat  less  physics  but  aimed  at  the  exploitation  of  more 
elaborate  mathematical  tools  In  more  practical  geometrical 
configurations . 

The  prime  example  of  the  latter  type  of  model  Is  given 
by  potential  theory  as  applied  In  fluid  dynamics.   There  are 
almost  no  fluid  flows  which  can  be  correctly  described  purely 
within  the  scope  of  the  Irrotatlonal  theory.   But  the  enor- 
mous mathematical  power  that  can  be  brought  to  bear  on  the 
potential  equation  gives  this  subject  Its  Importance.   In 
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effect  one  develops  a  foundation  for  the  study  of  fluid 
dynamics  through  a  profound  mathematical  Intuition  about 
a  nonexistent  fluid.   Whenever  possible,  a  real  fluid  flow 
is  thought  of  In  terms  of  deviations  from  this  Ideal  and 
completely  understood  model.   The  value  of  this  ideal  and 
physically  primitive  model  depends  first  on  the  extent  of 
the  mathematical  theory  and  second  on  a  knowledge  of 
exactly  where  this  theory  needs  radical  modification  (e.g., 
at  boundary  layers,  wakes,  shocks). 

Because  of  the  complexity  and  scope  of  plasma  physics, 
we  should  hope  to  be  able  to  develop  several  such  mathe- 
matically accessible  models  to  give  information  in  various 
parameter  ranges.   One  such  example  is  the  perfectly  con- 
ducting, scalar  pressure  macroscopic  theory  (MH).   This 
is  well  on  its  way  to  achieving  a  status  comparable  to 
classical  fluid  dynamics  [l],  except  that  the  second  pre- 
requisite, investigation  of  innumerable  boundary  layers, 
is  quite  incomplete  at  the  present  time. 

Another  candidate  for  this  role  as  the  source  of  a 
strong  mathematical  Intuition  about  plasmas  is  the  partly 
microscopic  and  partly  macroscopic  guiding  center  plasma 
(GCP)  and  the  associated  fully  macroscopic  but  anisotropic 
guiding  center  fluid  (GCF).   These  models  are  based  on  the 
simplified  asymptotic  description  of  a  charged  particle 
orbit  in  the  limit  of  small  gyro  radius.   We  recall  that 
there  is  a  standard  procedure  for  the  formulation  of  an 


exact  kinetic  description  (Liouvllle  equation)  when  given 
the  equations  of  motion  of  an  Individual  particle.   And 
there  are  standard  (though  heuristic)  procedures  for  reducing 
a  kinetic  model  to  a  macroscopic  one.   The  kinetic  result, 
when  carried  out  with  the  "exact"  particle  equations  is 
called  the  Vlasov  equations.   The  simplest  procedure  in 
deriving  the  GCP  equations  is  to  accept  as  given  the  orbits 
that  are  obtained  as  the  leading  term  in  the  guiding  center 
asymptotic  expansion  of  the  ordinary  differential  equations 
of  motion.   We  can  then  write  down  the  Liouvllle  equation 
for  these  artificial  orbits.   A  more  intricate  procedure  is 
to  introduce  the  guiding  center  expansion  directly  into  the 
Vlasov  equations.   This  amounts  to  repeating  the  guiding 
center  orbit  expansion  in  the  context  of  the  Vlasov  equa- 
tions, but  since  the  expansion  is  so  singular,  the  calcula- 
tions present  themselves  quite  differently. 

The  two  methods  of  obtaining  the  GCP  equations  were 
introduced  simultaneously  and  Independently,  the  simpler  by 
Grad  [2]  (later  by  Sagdeev  et  al.  [3]),  and  the  more  complex 
by  Goldberger  [^].   Extensions  of  the  GCP  theory  using  the 
second  method  were  made  by  Chew,  Goldberger,  and  Low  [5];  by 
Brueckner,  Chandrasekhar ,  Kaufmann,  and  Watson  [6],  [?]>  [8]; 
by  Rosenbluth  and  Longmire  [9],    and  Rosenbluth  and  Rostoker 
[10];  and  following  the  simpler  procedure  by  Grad  [11]  and 
Kulsrud  [12].   In  a  slightly  different  direction,  varia- 
tional and  quasl-varlatlonal  formulations  of  the  theory 
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have  been  made  by  Kruskal  and  Oberman  [13],  Kulsrud  [1'4-],  [12], 
Trubnlkov  [15],  and  Grad  [l6].   In  perturbations  about  a 
uniform  magnetic  field  the  GC  limit  Is  treated  In  [3],  by 
Chandrasekhar  et  al.  flY],  by  Rose  [l8],  and  by  Kadlsh  [19]- 

We  shall  give  a  brief  account  of  most  of  the  GOP  and  GCF 
theory,  for  which  there  Is  a  modicum  of  existence  theory  (and 
more  of  nonexistence.'),  paying  particular  attention  to  recent 
developments  In  stability  theory.   The  latter  Is  by  far  the 
most  subtle  and  difficult  part  of  the  subject,  as  witnessed 
by  the  fact  that  the  literature  Is  In  large  part  contradictory. 
For  further  details  In  the  variational  stability  theory  we  refer 
primarily  to  [l6],  [29],  [30],  also  to  [22],  and  [24]-[27]. 

Probably  the  most  Interesting  feature  of  the  GOP  Is  that 
It  Is  microscopic  In  only  one  dimension,  along  the  field  lines. 
This  Is  an  Immediate  consequence  of  the  singular  orbit  expan- 
sion which  yields  a  constraint  (first  order  system)  for  the 
perpendicular  motion  and  a  conventional  second  order  equation 
only  for  the  motion  In  the  parallel  direction.   A  further  source 
of  singular  behavior  In  the  GG  limit  Is  the  fact  that  momentum, 
mnv.  Is  multiplied  by  the  large  expansion  parameter  e/m  to 
become  current,  env .   As  a  result,  the  fluid  velocities  and 
electric  currents  that  appear  In  the  GCP  are  Independent  con- 
cepts, as  distinguished  from  the  Vlasov  theory. 

It  Is  not  surprising.  In  view  of  the  very  singular  GC 
limit,  that  the  GCP  and  GCF  equations  turn  out  to  be  mathe- 
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matically  unacceptable  In  some  parameter  ranges.   However, 
the  nonexistence  of  solutions  under  some  circumstances  can 
be  turned  to  advantage.   We  Interpret  this  phenomenon  as 
heralding  the  appearance  of  Instability.   This  Is  not  a 
logical  consequence  but  a  physical  Interpretation.   Moreover, 
If  we  credit  the  equations  with  complete  honesty,  they  are 
telling  us  that  the  growth  rate  of  these  Instabilities  Is 
too  rapid  to  be  described  by  the  theory  and  must  be  governed 
by  parameters  (such  as  gyro  frequency  and  plasma  frequency) 
which  have  been  set  equal  to  Infinity.   Such  an  Instability 
has  been  traditionally  called  a  micro -Inst ability. 

Another  distinguishing  feature  of  the  micro-Instabilities 
within  the  GCP  theory  Is  that  they  are  localized  physically. 
If  unstable,  they  can  be  excited  by  arbitrarily  localized 
disturbances.   If  micro-stable,  then  any  remaining  Instability 
Is  necessarily  global.   The  GCP  therefore  has  the  pleasant 
feature  of  distinguishing  categorically  between  micro-  and 
macro-lnstabllltles ,  both  In  localization  and  In  growth  rate. 

Mlcro-lnstablllty  makes  Its  appearance  In  three  distinct 
ways:  In  both  the  Initial  value  problem  and  In  the  boundary 
value  problem  for  static  equilibrium  by  violating  well-posedness 
or  continuous  dependence  criteria;  and  In  non-existence  of 
solutions  of  the  Integral  equation  which  governs  the  charge 
neutralizing  potential. 

There  are  several  sets  of  auxiliary  assumptions  and 
approximations  that  can  be  made  within  the  framework  of  the 
GO  theory.   For  example,  displacement  current  can  be  kept  or 
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dropped.   The  latter  yields  a  more  elegant  Galilean  invariant 
system,  is  simpler,  and  turns  out  to  be  not  more  restrictive. 
In  some  treatments  higher  order  terms  in  the  GC  expansion  are 
kept.   This  can,  in  principle,  extend  the  range  to  cover  new 
phenomena  (in  the  realm  of  micro-stability).  But  it  also  restricts 
the  range  of  validity  in  other  respects  (both  micro-  and  macro- 
stability).   If  done  without  care,  a  partial  contribution  from 
the  higher  order  GC  terms  can  narrow  the  domain  of  validity  of 
the  equations  in  one  direction  without  any  compensating  advan- 
tages in  others  (e.g.,  in  the  analyses  [21],  [23]  of  interchange 
stability) .   We  shall  distinguish  between  a  two-fluid  treatment 
and  a  one-fluid  treatment  which  does  not  differentiate  electrons 
from  ions.   The  distinction  is  subtle,  and  combining  the  two 
fluids  can  only  be  considered  a  marriage  of  convenience.   All  of 
the  references  cited  except  for  Grad  [11],  [l6],  [29]  and 
Andreoletti  [28]  are  single  fluid  (or  incorrectly  formulated  as 
two  fluids).   An  adiabatic  system  of  equations  is  Introduced  in 
Sec.  5  based  on  an  entirely  different  set  of  state  variables  than 
the  normal  GOP.   It  is  macroscopic  in  the  sense  that  magnetic 
field  displacements  and  motions  alone  enter;  but  it  is  distinct 
from  the  GOP.   These  equations  are  related  to  the  second  adia- 
batic invariant  and  they  can  give  information  about  some  classes 
of  instabilities  if  interpreted  carefully. 

While  our  main  purpose  is  to  present  a  coherent  account  of 
various  results  which  have  already  appeared  in  the  literature 
(but  in  conflicting  and  inconsistent  expositions),  we  also 


In  [3]  and  [17]  the  specialization  to  a  two-temperature 
Maxwellian  sidesteps  this  distinction. 
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present  several  new  results.   These  are  principally  the 
theory  of  the  neutralizing  potential  (Sees.  10  and  12), 
and  most  of  the  equilibrium  theory  (Sees.  9  and  10);  also 
the  Galilean  representation  of  the  Vlasov  equations  (See.  4), 
the  treatment  of  Interchange  stability  at  the  plasma  edge 
(See.  l4)  and  the  adiabatic  equations  (See.  5)- 


2 .   Plasma  Stability 

The  subject  of  plasma  stability  is  very  rich  in  the 
choice  of  alternative  procedures  and  even  equations  of 
motion.   Within  the  GO  theory  there  are  many  definitions  of 
stability,  some  explicit  and  some  only  implicit,  frequently 
noncomparable  and  often  contradictory.   To  some  extent  this 
proliferation  is  necessary  because  of  the  difficulty  of  the 
subject.   If,  for  example,  boundedness  of  solutions  of  the 
nonlinear  equations  of  motion  is  an  unattainable  goal,  some 
other  definition  of  stability  must  be  introduced.   Even  this 
"best"  definition  is  not  at  all  capable  of  answering  the 
questions  which  would  be  asked  by  the  designer  of  a  plasma 
containment  experiment.   But  totally  unnecessary  difficulties 
arise  when  the  formulations  and  even  the  equations  of  motion 
are  not  specified. 
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As  we  have  mentioned,  there  are  essentially  no  results 
known  with  regard  to  boundedness  of  solutions  In  the  non- 
linear GOP  theory.   Linear  boundedness  Is  not  comparable 
with  nonlinear,  so  It  must  be  taken  as  an  Independent 
definition  of  stability.   But  this  also  Is  frequently  too 
difficult,  and  we  retire  to  the  still  simpler  concept  of 
exponential  stability,  the  presence  or  absence  of  exponen- 
tially growing  normal  modes.   Explicit  counterexamples  can 
be  given  for  the  GGP  of  the  existence  of  growing  solutions 
when  there  are  no  eigenvalues  In  the  right  half  plane. 

A  still  more  special  concept  Is  that  of  marginal 

2 
stability.   This  requires  that  o)  have  a  real  spectrum,  that 

no  modes  te    can  arise  from  multiple  poles,  and  that  od  >  0 

2 
transfer  to  cu  <  0  only  by  passing  through  the  origin.   In 

other  words,  the  motional  operator  should  be  self  adjoint, 
the  continuous  spectrum.  If  any,  should  not  approach  the 
origin,  and  the  spectrum  should  be  bounded  and  depend  con- 
tinuously on  parameters.   Each  of  these  requirements  Is 
violated  within  the  GGP  theory. 

Almost  all  of  plasma  physics  Is  blessed  with  the  presence 
of  continuous  spectra.   The  common  streaming  operator,  ^Bf/5x, 
has  the  entire  Imaginary  axis  as  a  continuum.   The  Polsson 
bracket  [H,f]  which  describes  the  GC  motion  along  a  magnetic 
line  has  as  a  continuum  all  transit  frequencies  cd  of  particles 
between  turning  points.   This  continuum  covers  the  entire 
co-axis  (in  particular  Including  the  origin)  since,  In  the 


restoring  potential  |iB,  we  have  0  <  i-l  <  oo.   Although  the 
complete  GC  operator  couples  the  Poisson  bracket  [H,f] 
with  other  operators,  this  will  not  disturb  the  con- 
tinuum. 

As  to  continuous  dependence  of  the  spectrum,  consider 
the  anisotropic  wave  equation  (the  right  side  Is  self  adjoint) 

u,  ,  =  c-,u   +  c^u  (2.1) 

tt     1  XX     2  yy  ^         ^ 

In  a  rectangular  domain  with  u  -  0  at  the  boundary.   We  shall 
see  (Sec.  5)  that  this  example  Is  directly  related  to  stability 
of  a  uniform  GCP.   Suppose  that  Cp  approaches  zero;  Cp  =  0  Is 
marginal.   The  lowest  eigenvalue  Is 

2     2,  2/^2     2/  ?x  ,   2  2/^-d  ,^    _^ 

CD   =  TT  (c^/L^  +  Cp/Lp)  >  TT  c^/L^  (2.2) 

and  It  remains  bounded  away  from  zero.   It  Is  the  high  order 

2 
eigenvalues  that  change  sign  with  Cp  by  passing  through 

Infinity  rather  than  the  origin. 

We  should  distinguish  the  marginal  criterion  00  — >•  0 
(which  Is  required,  e.g.,  to  justify  use  of  the  second  adla- 
batlc  Invariant)  from  the  distinct  condition  oo/k  — >  0  which 
occurs  In  various  complex  analyses  of  stability  and  Is  more 
frequently  satisfied  In  a  transition  from  stability  to 
Instability.   In  the  uniform  plasma  [as  In  (2.1)]  the  disper- 
sion relation  Involves  only  the  ratio  cu/k.   Therefore  a 
transition  to  Instability  which  occurs  for  large  k  can  be 
discovered  by  setting  go  =  0,  even  though  the  transitional  mode 
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Is  not  slow.   We  shall  find  that  some  calculations  based  on 
the  concept  of  marginal  stability  give  correct  answers  and 
others  give  incorrect  answers.   To  distinguish  one  from  the 
other  it  is  necessary  to  abandon  this  concept  and  examine  the 
question  from  a  less  precarious  position. 

An  entirely  different  class  of  stability  analyses  is 
concerned  with  variational  methods  instead  of  boundedness  of 
solutions.   The  basic  argument  is  drawn  from  the  existence  of 
an  energy  integral,  K+^  =  const,  where  K  >  0  is  kinetic  and  ^ 
is  potential.   If  5   is  a  minimum  value  (subject  to  appropriate 
constraints),  then  any  perturbed  motion  which  has  an  energy 
constant  ^  +  6  will  satisfy  the  inequalities 


-o 


0  <  K  <  5  ,   J   <  5  <  J  +5  .  (2  3) 


-o    -^    -LO 

This  we  can  interpret  as  stability.   The  precise  consequences 
depend  on  what  departures  from  equilibrium  are  allowed  by  a 
small  value  of  K  and  a  prescribed  excursion  of  5  from  its 


minimum.   We  take  minimum  ®   (with  a  suitable  choice  of  neigh- 

j-o  ^  ^ 

borhood)  as  the  definition  of  variational  stability.   To 
correspond  to  the  classical  case  of  finite  degrees  of  freedom, 
stationary  ^  should  give  an  acceptable  definition  of  equilibrium. 
The  connection  with  boundedness  of  solutions  is  not  simple  (in 
either  direction,  necessary  or  sufficient),  and  the  results 
depend  crucially  on  completing  the  definition  by  a  specifica- 
tion of  the  class  of  admissible  functions. 

As  in  the  case  of  boundedness,  there  is  no  direct 
connection  between  the  linear  and  nonlinear  formulations. 
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here  between  existence  of  a  local  minim-urn  and  positivity 
of  the  second  variation.   We  adopt  6"^  >  0  as  the  definition 
of  linear  variational  stability.   We  shall  give  several 
examples  of  situations  which  are  linearly  stable  but  are 
unstable  with  respect  to  finite  displacements  of  arbitrarily 
small  amplitude  (Sees.  11  and  l4). 

There  is  a  direct  connection  between  the  sign  of  the 
second  variation  and  linear  boundedness  of  solutions  if  the 
system  is  self -ad joint  and  the  spectrum  is  discrete.   In 
this  case  the  normal  modes  are  obtained  by  minimizing  6  ^ 
with  respect  to  the  independent  variation  subject  to  suitable 
constraints  (Rayleigh's  principle).   This  connection  is  lost, 
even  with  a  self-adjoint  system,  if  there  is  a  continuum 
which  extends  to  the  origin.   One  can  have  unbounded  solutions, 
in  this  event,  even  with  a  strictly  positive  definite  5  ^. 
On  the  other  hand,  6  ^  >  0  is  equivalent  to  exponential  bound- 
edness when  the  system  is  self -ad joint ,  provided  that  we 
define  boundedness  to  be  in  Lp .   It  is  not  clear  whether  there 
is  any  definite  connection  between  5  ^  >  0  and  exponential 
boundedness  when  the  system  is  not  self -adjoint ;  and  one  can 

o  — 

show  by  explicit  counterexamples  that  6  _^  >  0  is  neither 
necessary  nor  sufficient  for  simple  boundedness  of  solutions. 
Although  self -ad jointness  has  not  been  established,  a 


There  is  reason  to  believe  that  such  degenerate  behavior 
at  the  origin  may  hide  within  the  GCP  theory  certain  types  of 
instabilities  (e.g.,  "drift")  which  might  seem  to  be  excluded 
from  a  lowest  order  GO  theory. 
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closely  related  formal  symmetry  of  the  single  fluid  GCP 
equations  of  motion  has  been  established  by  Kulsrud  for  an 
unperturbed  equilibrium  which  is  a  monotone  function  of 
energy,  Bf^/Se  <  0   [12].   From  the  associated  variational 
principle,  it  can  be  seen  that  the  continuous  spectrum 
extends  to  the  origin  under  most  circumstances.   And  when 
this  monotonicity  condition  is  violated  (as  it  often  is 
experimentally),  explicit  examples  can  be  given  of  normal 

■X- 

modes  with  complex  exponents. 

The  variational  approach  gives  an  improved  tool  to 
replace  the  discredited  notion  of  marginal  stability  in  the 
pessimistic  variation.   The  marginal  analysis  is  based  on 
the  second  adiabatic  invariant,  which  is  a  dynamical  constant 
for  slow  (marginal)  motions.   This  invariant  shares  two  dis- 
tinct properties  which  are  basic  to  the  entropy.   One  is  the 
dynamical  property  of  being  constant  for  a  slowly  varying 
system.   The  second  property  is  convexity  as  a  function  of 
energy  which  characterizes  equilibrium  as  the  state  of 
minimum  energy  at  fixed  entropy.   The  first  property  suggests 
the  marginal  stability  hypothesis.   The  second  leads  to  the 
pessimistic  variation.   The  pessimistic  variation  and  the 
second  adiabatic  argument  lead  to  identical  computations  in 
some  cases  and  to  quite  different  conclusions  in  others. 
Where  the  transition  to  Instability  does  occur  through  the 
origin,  the  adiabatic  calculation  will  be  correct.   Where 
the  adiabatic  calculation  gives  a  correct  answer  even 
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though  the  transitional  state  Is  not  In  slow  motion  [cf. 
(2. 2)], it  is  only  because  it  happens  to  agree  with  the 
pessimistic  calculation.   In  cases  where  it  is  obvious  that 
one  cannot  expect  adiabatic  behavior,  such  as  with  multiple 
turning  points  or  in  a  uniform  infinite  plasma  (the  invari- 
ant cannot  be  defined),  only  the  pessimistic  variation  is 
available.   But  there  are  less  obvious  cases  where  the  adia- 
batic invariant  can  be  defined  but  where  the  two  procedures 
disagree  violently,  e.g.,  with  non-monotone  f  (e). 

The  significance  of  the  monotonicity  criterion  is 
shrouded  In  ambiguity.   One  reason  is  that  there  are  two 
independent  stability  criteria  involving  partial  derivatives 
of  f  with  respect  to  energy,  Bf/5e,  holding  different 
variables  fixed.   The  two  separate  monotonicity  requirements 
cure  two  distinct  types  of  instability;  in  one  case   as  a 
sufficient  condition  for  interchange  stability  (Sec.  l4), 
in  the  other  as  a  necessary  condition  for  local  stability 
(Sec.  13)-   To  complete  the  confusion,  there  is  an  important 
special  case  where  the  two  partial  derivatives  coalesce. 
And  the  second  adiabatic  invariant  argument  which  greatly 
simplifies  the  analysis  of  one  instability,  completely  misses 
the  other. 

In  addition  to  the  cross-confusion,  each  of  the  partial 
derivatives  Sf/Se  is  subject  to  its  own  individual  misunder- 
standings.  With  regard  to  interchange  stability,  the 
appropriate  derivative  Sf/Se  <  0  is  correctly  taken  as  a 
sufficient  condition.   But,  even  though  the  literature  abounds 
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with  counterexamples,  it  is  Incorrectly  claimed  to  be 
sufficient  for  absolute  stability  (even  where  an  author 
himself  quotes  a  counterexample). 

Turning  to  local  stability,  we  find  three  distinct 
attitudes  toward  the  criterion  Sf/Se  <  0.   It  is  not  men- 
tioned at  all  as  a  requirement  for  stability  by  Rosenbluth 
and  Rostoker  in  a  marginal  analysis  [10]  and  by  Taylor  [21] 
and  Andreoletti  [26]  in  adiabatic  arguments  (but  the  latter 
two  require  the  other  derivative,  Sf/de  <  0,  for  inter- 
change stability).   It  is,  in  fact,  impossible  for  any 
motion  or  variation  which  preserves  the  area  of  each  energy 
contour  (adiabaticity )  to  excite  any  instability  which  is 
generated  by  an  inversion  of  f -values.   The  monotonicity 
condition  is  required  by  Kruskal  and  Oberman  [15]  but 
without  reference  to  stability,  simply  to  be  able  to  carry 
out  the  calculations;  they  are  neutral  parties  to  the 
dispute.   On  the  other  hand,  a  variational  analysis  (using 
the  pessimistic  instead  of  the  adiabatic  variation)  char- 
acterizes all  non-monotone  distributions  as  unstable  [l6]. 

An  accurate  study  of  boundedness  of  solutions  of  the 
equations  of  motion  should  give  a  result  in  between  the 
completely  optimistic  marginal  or  adiabatic  view  and  the 
completely  pessimistic  variational  opinion.   Partial  results 
are  known  for  non-monotone  f  only  for  the  infinite  uniform 
plasma  [19  ]  .   The  question  for  a  contained  plasma  (using 
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boundedness  as  the  definition  of  stability)  Is  not  likely 
to  yield  lightly  to  analysis.   However,  some  indications  are 
available.   Lack  of  existence  of  the  neutralizing  potential 
is  as  much  a  blow  to  the  equations  of  motion  as  to  the 
variational  formulation;  this  gives  certain  necessary 
criteria  on  the  degree  of  non-monotonlclty  of  f(e)  [Sees. 
10  and  12].   The  adlabatlc  equations  (Sec.  5)  are  strictly 
interpreted  as  governing  stability  only  for  monotone  f(e). 
But  they  describe  a  valid  approximation  to  an  actual  motion 
whenever  it  is  truly  marginal  (whatever  is  the  sign  of 
3f/Se).   Therefore,  a  prediction  of  adlabatlc  instability 
is  a  real  instability  if  we  independently  verify  that  cd  is 
zero  at  the  transition.   In  a  uniform  plasma  we  have  seen 
that  an  instability  that  transfers  at  large  k  can  be  found 
(without  any  marginal  connotation)  by  setting  ud  =   0 .      This 
indicates  that  a  transfer  of  stability  predicted  by  the 
adlabatlc  equations  for  small  wavelength  disturbances  may 
also  be  legitimate.   Even  in  the  infinite  plasma  there  are 
instabilities  associated  with  non-monotone  f  which  cannot 
be  found  in  this  way.   In  a  contained  plasma,  a  little 
reflection  shows  that  there  are  probably  many  more. 

A  variational  analysis  has  several  advantages  over  one 
involving  direct  examination  of  the  motion.   It  is  suffi- 
ciently simpler  so  that  some  nonlinear  results  become 
accessible  (cf.  Sees.  11  and  l4 ) .   It  is  also  more  complete 
in  that  a  natural  boundary  condition  is  supplied  if  one  is 
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Inadvertently  or  intentionally  omitted.   But  the  chief 
advantage  Is  that  we  have  a  certain  freedom  In  choosing 
the  admissibility  class  according  to  mathematical  exped- 
iency.  This  avoids  the  necessity  of  a  detailed  (and 
■usually  Impossible)  study  of  exactly  which  states  are 
accessible  via  an  actual  motion.   It  may  require  dropping 
a  constraint  which  is  known  to  govern  all  motions  in  order 
to  convert  a  lower  bound  to  an  attained  minimum  (cf.  the 
pessimistic  variation.  Sec.  11).   The  admissibility  class 
should  be  wide  enough  to  Include  all  states  which  can  arise 
in  a  motion,  and  narrow  enough  not  to  unduly  restrict  the 
class  of  equilibria. 

There  are  several  analyses  of  GCP  stability  which  lie 
somewhere  between  a  variational  and  a  full-blown  study  of 
actual  motions,  viz.,  by  Trubnlkov  [15]  and  by  Kruskal  and 
Oberman  [13]  •  Trubnlkov  considers  the  time  derivatives  K 
and  K  of  the  kinetic  energy  of  a  motion  as  it  passes  through 
a  state  of  equilibrium  at  t  =  0.   He  verifies  that  K  =  0 
and  takes  K  <  0  as  the  definition  of  stability.   At  first 
glance  this  looks  very  similar  to  the  condition  &  ^  >  0. 
One  basic  objection  in  principle  is  that  it  says  nothing 
about  the  stability  of  a  state  which  is  close  to  equilibrium 
but  which  is  not  an  Immediate  antecedent  or  successor  in 
time  to  the  exact  equilibrium.   It  is  exactly  this  question 
of  accessibility  of  states  during  a  motion  which  is 
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avoided  by  taking  the  more  contrete  admissibility  class 
in  a  true  variational  analysis. 

A  second  more  elementary  difficulty  is  that  a  growing 
oscillation  (complex  eigenvalue)  would  inevitably  be  clas- 
sified as  stable.   Thus  self -ad jointness  seems  to  be  a 
prerequisite . 

But  perhaps  the  most  interesting  difficulty  with  the 
Trubnikov  procedure  is  an  acute  sensitivity  to  how  one 
separates  the  variables  which  describe  the  plasma  state  into 
coordinates  and  rates .   In  a  classical  Hamiltonian  formula- 
tion there  is  an  unambiguous  one-to-one  correspondence 
between  coordinates  and  momenta.   The  coordinates  would  be 
fixed  at  an  equilibrium  and  K  would  be  evaluated  for  all 
initial  momenta.   We  describe  three  different  choices  of  the 
coordinates  that  can  be  taken  for  the  GCP,  to  each  of  which 
one  can  apply  the  Trubnikov  procedure. 

The  most  obvious  course  of  action  is  to  take  the  magnetic 
field  configuration  and  the  distribution  function  as  coordi- 
nates, leaving  only  the  magnetic  line  velocity  as  a  rate. 
This  is  an  unbalanced  but  natural  choice.   If  we  follow 
Trubnikov ' s  recipe  with  this  assignment,  we  do  not  get  his 
answer,  but  we  obtain  a  much  more  special  formula,  viz.,  the 
transverse  second  variation  of  Sec.  11.   This  is  a  special 
variation  in  which  f  is  not  varied  at  all;  it  arises  in 
this  context  because  not  only  K  but  even  f  (which  is  equivalent 
to  the  first  variation  of  f)  vanishes  at  t  =  0. 
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Trubnikov  actually  evaluates  K  using  a  special  approxi- 
mation Involving  the  second  adiabatlc  invariant,  and  in 
doing  so  obtains  the  pessimistic  second  variation  of  Sec.  11. 
But  this  is  not  an  incidental  approximation;  it  amotmts  to 
a  different  choice  of  the  state  variables  and  of  the  equa- 
tions of  motion.   Taking  the  adiabatic  invariant  as  constant 
implies  that  a  unique  distribution  function  is  assigned  to 
any  given  magnetic  configuration.   The  distribution  function 
is  no  longer  an  independent  coordinate,  and  the  resulting 
adiabatic  equations  of  motion  [pessimistic  if  f  (e)  is 
monotone]   involve  only  the  magnetic  field.   The  dis- 
placement and  velocity  of  field  lines  form  a  balanced 
set  of  coordinates  and  momenta. 

A  third  choice  of  coordinates  is  that  taken  by  Kulsrud 
in  his  symmetric  version  of  the  equations  of  motion  (see 
Sec.  5)-   Symmetry  is  found  by  splitting  f  =  f  +  f  into 
even  and  odd  components  in  the  parallel  component  of  velocity, 
taking  field  displacement  and  f"  as  coordinates  and  field 
velocity  and  f  as  momenta.   The  Trubnikov  method  applied 
here  gives  still  a  different  result;  the  full  second  varia- 
tion  6  ^  but  with  f  instead  of  f  as  the  argument. 

The  Trubnikov  recipe  is  very  flexible  but  it  seems  to 
require  independent  confirmation  of  the  calculation  by  a 
more  fundamental  method. 

Several  interesting  and  ingenious  definitions  of 
stability  have  been  introduced  by  Kruskal  and  Oberman.   The 
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first  criterion  is  simply  that  the  given  equilibrium  be 
the  unique  state  which  is  compatible  with  "all  regular 
time -Independent  constants  of  the  motion."   In  a  classical 
system  with  finite  degrees  of  freedom,  uniqueness  is  essen- 
tially equivalent  to  stability.   The  significance  with 
respect  to  more  complicated  systems  is  not  clear.   Unique- 
ness is  investigated  by  varying  the  total  energy,  kinetic 
plus  potential,  subject  to  the  stated  constraints.   At  a 
later  stage  this  criterion  is  altered  by  removing  the 
positive  kinetic  and  electrostatic  energy  and  examining  the 
sign  of  the  remainder,  called  5¥.   To  investigate  the  value 
of  6W  subject  to  all  constants  of  the  motion  means  to  take 
as  the  admissible  class  only  actual  motions.   This  would  be 
Impossible  to  pursue.   But  the  admissible  class  is  unwit- 
tingly widened  in  two  directions.   First  the  equation 
governing  the  motion  of  the  field  lines  (momentum  conservation) 
Is  Ignored,  making  the  field  displacement  and  magnetic  line 
velocity  independent  of  the  distribution  function  and  losing 
undetermined  constants  of  the  motion.   If  we  were  to 
literally  keep  all  remaining  constants,  the  admissible  class 
would  now  be  a  class  of  virtual  motions  in  which  the  field 
variation  is  chosen  arbitrarily  as  a  function  of  time,  and 
the  distribution  function  evolves  in  accordance  with  the 
correct  kinetic  equation.   This  procedure  Is  the  same  as 
Trubnlkov ' s  in  the  first  version  described  above  and  would 
lead  to  the  transverse  variation  for  5W.   But  Kruskal  and 
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Oberman  actually  constrain  f  by  a  very  specific  analytic 
expression  which  Is  far  from  all  Inclusive,  and,  In  so  doing, 
obtain  an  expression  comparable  to  the  full  second  variation, 
6  |.   A  final  formal  minimization  with  respect  to  the  distri- 
bution function  yields  a  formula  similar  to  the  pessimistic 
variation.   Analysis  of  the  minimization  shows  that  the 
minimum  can  be  achieved  In  general  only  by  relaxing  even 
the  explicit  constraint  formula  of  Kruskal  and  Oberman. 

An  unusual  formulation  has  been  Introduced  by  Taylor  [21] 
and  Andreolettl  [23],  [28].   They  take  particle  orbits  as 
governed  by  constant  second  adlabatlc  Invariant  as  well  as 
magnetic  moment  and  take  flux  line  coordinates  as  the  dynami- 
cal variables.   These  dynamical  equations  of  motion  describe 
a  higher  order  GC  "drift"  of  particles  which  are  not  frozen 
to  magnetic  lines  (as  they  are  in  the  GC  theory  described 
in  the  present  paper).   But  this  is  combined  with  a  flux 
equation  governing  the  magnetic  field  variation  [see  (3-9)] 
which  is  compatible  only  with  particles  which  are  frozen  to 
magnetic  lines.   Energy  variations  are  computed  for  stability 
purposes  assuming  that  the  particles  are  frozen;  in  fact, 
for  the  partially  higher  order  GC  approximation  used,  there 
is  no  relevant  energy  conservation  law.   Since  the  dynamical 
(drift)  equations  are  used  only  to  delimit  the  class  of 
equilibria,  not  to  govern  motion  or  stability,  all  results 
derived  from  this  procedure  are  equivalent  to  a  marginal  zero 
order  GC  theory  or  to  the  adlabatlc  equations  of  motion 
(Sec.  5)- 

It  should  be  clear  from  this  discussion  that  the 
unadorned  terms  stable,  unstable,  necessary,  sufficient 
are  subject  to  an  Infinite  variety  of  interpretations.   We 
consider  this  not  as  a  fault  but  as  evidence  of  the  Infinite 
richness  of  the  subject.   But  this  wealth  must  elude  us  if 
it  is  not  subject  to  acceptable  auditing  standards. 
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3-   Orbits  and  Electromagnetic  Potentials 

The  classical  equations  of  motion  of  a  charged  particle 


are 

m  ^  =  e(E  +  V  X  B)  (^-l) 

where  E(x,t)  and  B(x,t)  are  assumed  to  be  given.   The  particle 
motion  Is  closely  tied  to  Maxwell's  equations.   For  example, 
a  Lagranglan  exists  only  If  E  and  B  satisfy  the  two  homo- 
geneous Maxwell  equations, 

dlv  B  =  0 

(3.2) 

||-  +  curl  E  =  0 
In  this  case,  the  Lagranglan  and  Hamlltonlan  are  given  by 

ji_  =  p-  mv     +   ev  •  A   -   ej?5 


(3.3) 


where 


B  =   curl   A 

(3.4) 
E.-|A.^^         • 

It  will  be  convenient  to  discuss  the  electromagnetic 
field  before  turning  to  the  orbits.   First  suppose  there  Is 
given  an  incompressible  vector  field,  dlv  B  =  0,  which 
varies  smoothly  as  a  function  of  a  parameter,  B(x,t).   In  a 
domain  which  is  simply  covered  by  the  vector  field.  Fig.  1, 
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FIGURE  1 


we   can   Introduce   flux  coordinates  a(x),    P(x)    such  that 

B  =    Va  >c7p    .  (5.5) 

On  a  transverse  surface,  B-dS  =  dadp  ?=  f lux.   Since  a  and  p  are 

constant  on  each  line,  they  can  be  considered  as  coordinates 

for  the  field  lines.   There  is  considerable  flexibility  in 

choosing  a  and  (3.   For  example,  a  can  be  chosen  arbitrarily 

on  one  transverse  surface  (S '  in  Fig.  l)  and  p  still  has 

some  freedom.   This  arbitrariness  holds  for  each  value  of  t, 

and  we  only  require  a(x,t)  and  p(x,t)  to  be  smooth  in  t. 

After  making  a  choice,  we  adopt  the  expedient  of  identifying 

fixed  values  a  =  a  ,  P  =  f3   as  the  same  line,  moving  with  t. 

o       o        ° 

It  now  becomes  possible  to  speak  of  the  velocity,  U ,  of  a 
given  line,  and  we  may  clearly  set  U-B  =  0  to  fix  U  uniquely. 
To  be  specific,  we  introduce  a  Lagrangian  parameter  a  along 
each  line  (e.g.,  arclength  at  t  =  O)  such  that  the  curve 
o   =  constant  on  a  given  moving  line  is  orthogonal  to  the 
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field  lines.   The  complete  set  of  Lagrangian  coordinates 
with  respect  to  the  flow  U  is  (a,p,a)  and  we  have 


+  U-  Va  -  0 


Da  _  Sot 

Dt  -  "St 

(3.6) 

where  S/Bt  and  V  refer  to  the  Eulerian  variables  (x,t). 
A  little  vector  manipulation  yields  the  formulas 

B  X  U  =  ^  \7a  -  1^  VP  (3.7) 


B 


ll  +  curl(B  X'  U)  =  0  (3-9) 


We  repeat  that  the  velocity  U  which  describes  a  given  field 
variation  B(x,t)  is  far  from  unique  and  depends  on  a  quite 
flexible  assignment  of  the  coordinates  (a,p)  at  each  t. 

Now  suppose  we  are  given  a  pair  of  vector  fields  E(x,t), 
B(x,t)  which  satisfy  the  field  equations  (5-2)  and  also 
E-B  =  0  for  all  t.   We  define  a  vector  U  by  U  =  E  X  B/B   or 

E  +  UXB  =  0  ,U.B  =  0  .  (3.10) 

With  this  definition,  U  and  B  satisfy  (3.9).   It  is  a  well 
known  consequence  of  this  differential  equation  that 
magnetic  lines  are  mapped  onto  magnetic  lines  by  the  motion 
U^  so  that  coordinates  (a,p)  which  are  carried  by  the  flow 
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can  be  Introduced,  B  =  N/a  ><•   V^   ^-^^  Da/Dt  =  Dp/Dt  =  0. 
By  the  previous  analysis  we  conclude  that 

E  =  II  Va  -  ^  V-P  .  (3.11) 

This  relation  holds  not  for  arbitrary  flux  coordinates  but 
only  for  a  proper  choice  of  (a,p)  since  U  has  been  uniquely 
chosen  in  this  case. 

We  present  an  independent  derivation  of  (3-11) 
without  making  use  of  the  theory  of  the  equation  (3.9). 
Using  the  previous  analysis  of  B  alone,  we  introduce  arbi- 
trary parameters  (a,p")  associated  with  a  magnetic  flow 
velocity  U.   Since  U ,  defined  in  (3.10),  and  U  each  satisfy 
(3-9) J  we  have 

From  E-B  =  0  we  find  B-Vfi   =  0  and  conclude  that  /i  is  a 
function  of  (a,p,t)  only.   Thus 


St    Sa        St   Sp 
If  we  now  choose  a  and  p  to  satisfy 


Da.  ^  hjf),  a) 

Dt  S(a,  p") 

Dp  ^  ^{f),  p) 

Dt  S(a,  p) 


(3.1V 
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and  the  initial  condition 

%4l=  1  ,  (3.15) 

S(a,p) 

then  a  small  computation  shows  that  E  has  the  form  (3-11). 

Equation  (3>1^)  for  a  and  (3  is  a  Liouvllle  equation  in  the 

(a,p)-plane  with  j25(a,p  ,t )  as  Hamlltonlan.   The  derivative 

D/Dt  is  to  be  taken  with  (a,p)  fixed.   Equation  (3-1^) 

states  that  the  values  of  a  and  p  are  carried  with  the  flow 

velocity  (-Sj^^/BpT,  -S/)/Sa") ,  and  the  Jacoblan  relation  (3-15) 

is  easily  seen  to  hold  for  all  t. 

Finally  we  turn  to  a  pair  of  vector  fields  E(x,t)  and 
B(x,t)  which  satisfy  (3-2)  and  E-B  =C(5)  where  6  is  an 
extraneous  parameter  taken  to  be  small.   If  we  define 
U  =  E  ><  b/B  ,  then  the  flux  conservation  equation  (3-9)  is 
satisfied  only  to  within  0{5),    and  magnetic  lines  are 
carried  by  the  velocity  field  U  only  to  within  0(6).   It 
will  be  much  more  satisfactory  to  modify  the  definition  of  U 
by  a  small  amount  0(&)  in  order  to  preserve  (3-9)  exactly. 

To  do  this  we  construct  a  function  j6    (x,t)  with  the 
property  that  E-B  =  Vj6    -B.   To  find  /)    ,   we  merely  integrate 
/)     =   /E'dx  along  each  magnetic  line  from  an  arbitrary  trans- 
verse  surface  on  which  ^     =   0.      From  the  construction,  ^ 
Itself  is  0(5).   The  vector  field 


E*  =  E  -  Vj6*  (3.16) 
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differs  from  E  by  0(&)  and  satisfies  Maxwell's  equation 
(3.2)  as  well  as  E  -B  =  0.   We  construct  the  velocity  field 
U  from  E  and  observe  that  (3.6)-(3-9)  are  satisfied 
together  with 

E  +  U  X  B  =  '7'f^* 


^  =  5l^«  -Ir^P^^^^ 


0-17) 


y  .y. 

The  parallel  component  of  S/ p     is  unique  but  P     itself  is 
not  unique  to  within  an  arbitrary  added  function  of  (a,p,t) 
[which  we  restrict  to  be  0(5)].   Thus  U  is  unique  to  within 

0(6). 

After  these  preliminaries,  we  turn  to  the  particle 

orbits.   To  be  definite  we  consider  m  to  be  fixed  and  study 

the  orbits  in  the  limit  of  large  e.   We  assume  that  E  and  B 

are  bounded,  B  is  bounded  away  from  zero,  and  require  that  v 

remain  bounded.   The  last  requirement,  bounded  energy,  implies 

that  the  component  of  E  in  the  direction  B  approaches  zero 

with  l/e.   In  other  words,  there  is  no  satisfactory  limiting 

behavior  as  e  — >  00  unless 

E.B  =^0(|)  .  (5.18) 

In  order  to  study  the  behavior  of  individual  orbits,  we 
merely  postulate  this  relation.   In  order  to  apply  the  result 
to  a  plasma  it  will  have  to  be  established  that  the  collec- 
tive particle  motions  give  rise  to  such  an  electric  field. 


-  2( 


This  we  leave  to  later. 

In  the  limit  of  large  e,  a  particle  will  execute  high 
frequency  tight  oscillations  about  a  guiding  center  (GC). 
The  asymptotic  description  governs  the  motion  of  the  GC 
and  the  energy  of  the  oscillation.   As  is  usual  in  an  asymp- 
totic expansion,  the  parameters  which  appear  in  the  asymp- 
totic description  (e.g.,  the  GO  coordinates)  are  not  simply 
related  to  the  variables  (x,v)  taken  from  the  exact  equa- 
tions.  This  is  just  the  familiar  question  of  assigning  a 
definite  asymptotic  series  to  given  initial  values.   We  do 
have  an  asymptotic  formula  which  relates  the  asymptotic 
coordinates  to  the  instantaneous  state,  and  in  a  few  special 
cases  (primarily  a  constant  magnetic  field  region),  the 
relation  is  exact. 

The  motion  of  a  GC  transverse  to  the  magnetic  field 
is  given  by  the  velocity  E  X  B/B  to  lowest  order.   This 
is  true  for  all  particles,  independent  of  mass,  sign  or 
magnitude  of  e,  or  energy.   Consistent  with  this  description, 
we  can  introduce  the  velocity  U  of  (3-17)  (not  exactly 
E  X  B/B  )  which  carries  magnetic  lines  into  themselves  and 
state  that  a  GC  is  constrained  to  lie  on  a  given  (moving  and 
deforming)  magnetic  line.   The  motion  of  the  GC  along  the 
magnetic  line  can  be  described  as  that  of  a  classical  con- 
strained  particle  subject  to  the  additional  potential  |i.B  +  ep 
where  |J.  is  a  constant.   The  potential  fi      [as  in  (3-17)]  Is 
unique  within  an  added  constant  on  a  given  line,  and  from 

-  27  - 


(3- 18)  and  the  preceding  analysis,  e/i   Is  of  order  -unity. 
The  constant  \i    ( magne tic  mome n t )  is  related  to  the  oscilla- 
tory energy  by  the  asymptotic  expression. 


^L 


=  ^  mv^B  .  (3.19) 


The  fact  that  \i   is  constant  is  the  equation  that  governs 
the  variation  of  the  oscillatory  energy;  the  variation  of 
B  on  a  trajectory  is  determined  by  the  GC  motion. 

The  Lagrangian  of  a  freely  moving  constrained  particle 

12    12 
is  p-  mv  +  p-  mU  where  the  second  term  is  to  be  Interpreted 

as  the  constraining  potential.   Therefore,  for  the  GC  motion 

along  a  line  we  have 

2 

L  =  ^m^+^m\]^-[iBTj6  (3.20) 


2    ^   2   "     ^^  + 


H  =  ^  C^P^  -  ^  mn^   +   [IB  t  j6      .  (3.21) 


We  have  Introduced  o,  the  Lagrangian  coordinate  following  U, 
and  with  it 

y,        ho 
^   -   ^s" 

v  =  |=^P^.  (3.22) 

Also,  we  have  incorporated  the  factor  e  into  the  potential 
^    (which  is  now  of  order  unity);  ejd     is  written  +j6   for  ions 
and  -^   for  electrons. 
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The  perpendicular  and  parallel  components  of  the 
electric  field  enter  quite  differently  into  the  GC  formulas. 
In  the  decomposition  (3-l6)  E  ,  of  order  unity,  determines 
the  transverse  GC  motion  through  E  +  U  X  B  =  0.   The  poten- 
tial  p    ,    of  order  1/e,  enters  the  lowest  order  Hamiltonian 
(3-21).   It  is  therefore  necessary  to  consider  E  and  e/5  =  ±  j6 
as  independent  lowest  order  entities,  forgetting  their 
common  origin.   We  shall  use  E  henceforth  instead  of  E  ; 
thus  E  is  a  vector  which  is  orthogonal  to  B,  and  />  is  an 
unrelated   scalar  which  is  defined  only  to  within  an  added 
constant  on  each  line.   The  entire  collective  plasma 
analysis  to  follow  will  treat  E  and  ^   distinctly. 

To  summarize,  the  leading  term  in  the  asymptotic  expan- 
tion  for  large  e  describes  a  GC  which  is  constrained  to 
follow  a  given  magnetic  line,  moving  along  the  line  according 
to  the  Hamiltonian  (3.21).   The  fact  that  \i,    appearing  in  the 
Hamiltonian,  is  a  constant  governs  the  oscillatory  energy. 
To  complete  the  identification  of  the  asymptotic  coordinates 
we  should  also  introduce  the  oscillatory  phase 

0  =  ^  /  B  dt  ,  (3.23) 

oriented  relative  to  fixed  (a,p).  The  six  coordinates  (x,v) 

are  related  in  an  elementary  way  [using  the  gyro  radius 

1/2 
A  =:  (2mii./B)  ^    /e  ]  with  the  GC  coordinates  (a  ,p  ,a  ,v  ,|j.,  0  )  .   The 

parameter  9   will  not  enter  into  any  of  the  subsequent  analysis. 

Most  orbit  analyses  go  much  deeper  into  the  theory  than 
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the  leading  term  which  Is  all  we  require  for  this  theory; 
general  references  are  [51]  and  [32]  and  for  proofs  [35]- 


k .      The  Vlasov  Equatlon5--Gharge  Neutrality  and  Galilean 
Invarlance 

We  take  as  a  kinetic  equation  the  Llouvllle  equation 
formed  with  the  "exact"  Hamlltonlan  {J>.^), 

Here  P(x,£,t)  Is  the  particle  density  In  the  phase  space 
(x,£)  of  a  single  particle.   There  are  two  equations  (4.1) 
for  F  and  P_  In  terms  of  )-^     and  }f-_.      The  two  kinetic 
equations  (4.1)  together  with  Maxwell's  equations  (3-2)  and 


1  Be 

+  cui'x  a   -   \^   0 

(4.2) 


~2^+   curl  B  =  H^J 


K  dlv  E  =  q  , 
o        ^ 

and  the  identification  of  the  electromagnetic  source  terms 


q  =  e(n^  -  n_) 


J  =  e(n  u   -  n_u_) 


(^•3) 
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(the  density  n  and  the  flow  velocity  u  are  the  usual  moments 
of  p)  comprise  the  Vlasov  equations.   There  are  some  special 
existence  theorems  for  this  system  of  integro-dlfferential 
equations,   and  there  is  a  certain  amount  of  evidence  to 
indicate  that  the  system  is  mathematically  sound  and  none  to 
the  contrary  at  the  present  time.   On  the  other  hand,  the 
Vlasov  equations  are  inelegant  since  they  are  invariant  under 
neither  Galilean  nor  Lorentz  transformations.   The  system  can 
be  made  Lorentz  invariant  by  introducing  a  relativistic 
Hamiltonian.   Or  it  can  be  made  Galilean  Invariant  in  a 
charge -neutral  approximation  as  described  below. 

It  will  be  convenient  to  distinguish  the  charge  and 
current  that  appear  in  Maxwell's  equations  (q,yrj  J^)  from  the 
moments  of  F  which  are  obtained  from  the  kinetic  equations 
(q  ,  Ji<-)  •   Each  set  independently  [before  making  the  identifi- 
cations q  =  q  ,  J   =  J   implied  by  (4,3)]  satisfies  charge 


conservation. 


|2.  +  div  J  =  0  .  (4.4) 


From  (4.4)  we  see  that  the  identification  J   =  J   implies 
qp.  =  qjyj  provided  that  the  latter  is  satisfied  initially.   It 
will  be  even  more  convenient  to  identify  Bj  /St  with  Sj^St, 
in  which  case  both  the  q  and  J  identifications  become  initial 
restrictions  only.   From  Maxwell's  equations  we  compute 


D.  Gorman  and  H.  Weitzner,  private  communication. 
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[X     ^   =  -  ^  ^  -  curl  curl  E  (4.5) 

and  by  taking  moments  of  the  kinetic  equations, 

c 

where  Q,   is  the  plasma  frequency  and  X  does  not  contain  E 
explicitly.   Interpreting  X,  qualitatively,  as  an  inhomo- 
geneous  term,  the  identification  3j„/5t  =  Bj^VSt, 

2  2 

-^   ^-4  +  curl  curl  E+^E+X=0  ,       (4.7) 
c'^  at^  c 

describes  electromagnetic  waves  coupled  with  plasma  oscilla- 
tions.  We  can  take  as  equivalent  to  the  original  Vlasov 
system  the  two  kinetic  equations  (4.1),  Maxwell's  equations, 
and  equation  (4.7)  (plus  initial  data)  to  replace  the  source 
identification  (4.3). 

Before  continuing,  we  remark  on  a  qualitative  inter- 
pretation of  the  Debye  length 

d  =  c_/n  ,  (4.8) 

(C_  is  a  mean  electron  thermal  speed)  and  the  magnetic 
shielding  length 

5  =  c/n  .  (4.9) 

The  Debye  length  is  the  largest  distance  over  which  a  finite 
fractional  depletion  of  ions  or  electrons  will  give  rise  to 
a  finite  electrostatic  field.   Any  reasonable  electric 
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field  on  a  much  larger  length  scale  can  be  created  by  a 
small  charge  unbalance  which  leaves  n   essentially  equal 
to  n  .   Intuitively  (and  in  common  finite  difference  schemes) 
one  considers  the  charge  separation  q  to  give  rise  to  E. 
But  if  d  is  a  macroscopically  small  quantity  {D,   is  large), 
this  is  a  very  sensitive  operation;  numerically,  the  time 
step  would  have  to  be  very  small  (comparable  to  1/Q,)  .      The 
charge  neutral  (or  quasi-neutral)  approximation  consists 
of  inverting  this  process.   We  set  n  =  n_  as  a  constraint 
and  compute  what  E  must  be  in  order  to  maintain  neutrality. 
The  constraint  serves  to  determine  the  evolution  of  E  rather 
than  the  time  derivative  of  E  (displacement  current)  which 
is  dropped. 

A  similar  situation  holds  with  regard  to  the  magnetic 
shielding  distance,  6.   It  can  be  interpreted  as  the  largest 
distance  over  which  a  finite  relative  motion  J  between  ions 
and  electrons  gives  rise  to  a  finite  magnetic  field  B.   Any 
magnetic  field  which  varies  on  a  much  larger  length  scale 
than  5  can  be  produced  by  a  current  which  represents  a  small 
relative  motion  (compared  to  mean  thermal  speed)  between  ions 
and  electrons.   Since  magnetic  forces  are  smaller  than 
electric  by  a  relativistic  factor,  the  relation  6/d  =  c/C_ 
is  revealing. 

Formally,  the  charge  neutral  approximation  in  Maxwell's 
equations  consists  in  dropping  displacement  current, 
curl  B  =  M-_J|yr5  and  completely  ignoring  the  source  relation 
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K  dlv  E  =  q„„.   We  retain  the  pre -Maxwell  system  [3^]  {J-.,r-^J-n) 
O  M  l"i    r 


1^  +  curl  E  =  0  ,  div  B  =  0  (^-10) 

curl  B  =  ki^Jp   ,   div  Jp  -  0  (^-H) 


To  be  compatible  with  the  kinetic  equations,  E  must  be 
chosen  to  insure  q  ==  0,  div  J„  =  0.   This  is  accomplished 
by  setting  the  pre-Maxwell  value 


[^ 


-  =  -  curl  curl  E  (^-12) 


o  ^ 

equal  to  [i   Sj„/St  as  given  In  (4.6).   The  equation  for  E 
(heurlstically  uncoupled  by  looking  at  X  as  an  "inhomogeneous" 
term)  is  now  elliptic. 

curl  curl  E  +  E/6^  +  X  -  0  .  (4- 13) 

If  X  is  given,  E  can  be  solved  in  an  interior  domain  with 
a  boundary  condition  on  E, ,  and  in  an  exterior  domain  with 
suitable  regularity  at  infinity.   If  by  any  means  we  obtain 
a  solution,  with  q  =  div  J  =  0  initially,  of  the  system  of 
kinetic  equations  (4.1),  pre-Maxwell  equations  (4.10),  and 
constraint  (4.13),  then  q  =  div  J„  will  vanish  for  all 
time  since  curl  curl  E  has  zero  divergence. 

There  are  no  proofs  of  the  mathematical  legiti- 
macy of  this  charge -neutral  plasma  formulation,  but  there 
is  a  certain  amount  of  empirical  evidence.   For  example. 
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there  is  numerical  evidence  of  a  satisfactory  solution  in 
one  dimension  for  a  piston  moved  into  a  cold  plasma,  except 
that  under  certain  circumstances  one  is  forced  to  refer  to 
a  generalized  weak  solution  with  Internal  discontinuities 

[35]. 

A  further  reduction  is  obtained  by  considering  "magnetic 
neutrality"  with  6  small  (in  a  non-relativistlc  plasma  this 
is  always  more  restrictive  than  small  d).   We  can  interpret 
6  as  a  singular  parameter  in  the  elliptic  equation  (4,13). 
Except  near  a  boundary  or  in  a  shock,  we  may  set 

Making  the  charge  neutral  approximation  in  X,  this  takes 
the  more  explicit  form 

— i-Tv  (E  +  u  XT  B )  =  u  .  v'J  +  J  div  u  +  J  .  V'u  +  — J  X  B 

+  div  Q  (4.15) 

m   -  m 

Q.  .  =  —  PT.  -  ^  p.  . +— =^  J.J.  (4.16) 

^ij   m_^   ij    m_   ij    ne(m_^  +  m_)   i  j  ^  ^ 

[this  expression  is  also  correct  for  Insertion  into  (4.13)]. 

Each  of  the  forms  (4.7),  (4.13),  (4.15)  can  be  con- 
sidered to  be  a  colllsionless  version  of  Ohm's  law.   It  is 
evident  by  inspection  that  (4.15)  (and  the  entire  coupled 
system  of  equations)  is  Galilean  invariant  provided  that  B 
and  J  are  taken  to  be  invariant  under  a  change  of  frame 
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and  E  transforms  such  that  E  +  u  x  B  is  Invariant.   A  slight 
computation  shows  that  the  expression  (4.13)  Is  similarly 
Invariant.   Transfer  to  a  frame  with  a  constant  relative 
velocity  u   gives  a  term  In  curl  curl  E  which  Is  canceled 
by  u^-\7J  In  (4.15). 

In  the  GC  limit  we  let  e  become  large  (note  that  there 
are  factors  e  in  f2,  l/d ,  and  1/6).   The  speed  of  light,  c, 
is  an  Independent  parameter;  taking  it  finite  means  using 
Maxwell's  equations  and  letting  it  become  large  yields 
the  pre-Maxwell  equations  (for  definiteness  we  can  take 
M-  =  1  and  k  =  l/c  ).    But  examination  of  the  complete 
source  identification  relation  (4.7)  yields  a  result  which 
is  Independent  of  c  to  the  two  lowest  orders  in  l/e . 
Whether  or  not  c  is  kept  in  Maxwell's  equations,  it  does 
not  appear  in  the  source  compatibility  equation.   To 
lowest  order  we  find  E  +  u  x.  B  =  0  (consistent  with  the 
GC  orbit  analysis),  and  to  the  next  order  we  would  need 
to  include  only  a  Hall  term  J  K  B  and  a  term  in  div  P  from 
(4.15).   The  parallel  component  of  E  vanishes  to  lowest 
order;  this  confirms  from  the  plasma  equations  the  postu- 
late that  was  needed  in  the  GC  orbit  theory.   The  first  order 
parallel  component  of  E  is  given  by  the  term  in  div  P, 
about  which  more  later. 

We  have  based  this  discussion  on  a  nonsteady  plasma  motion, 
The  analysis  must  be  reexamined  in  a  steady  flow  since  the 
casually  dropped  "Initial  restrictions"  become  dominant,  and 
dominant  terms  such  as  curl  curl  E  vanish. 
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5.   Guiding  Center  Plasma  Equations  of  Motion 

The  dominant  feature  of  the  GC  motion  is  the  con- 
straint to  a  given  magnetic  line.   Evidently  the  transverse 
velocity  U  is  macroscopic  since  every  particle,  whether  ion 
or  electron,  and  without  regard  to  mass  or  energy  moves 
similarly.   Any  kinetic  behavior  occurs  one-dimensionally 
along  each  magnetic  line.   To  describe  this  motion  we 
introduce  the  one -dimensional  Liouville  equation  with  a  GC 
Hamiltonian  (3-21), 

Here  f  is  a  function  of  (a,p  ,a  ,p,M.,t )  where  (a,p)  are  the 
dynamical  variables  and  (a,p,M.)  are  parameters.   There  are 
two  such  kinetic  equations  for  ions  and  electrons.   The  GC 
distribution  f  is  normalized  as  a  number  density  in 
Q,  =    (a,p,a,p,M.)  and  is  related  to  the  Vlasov  distribution 
P(x,p)  by 

f  =  2Tm\F  (5.2) 

through  the  Jacobian  relations 

dx  =  dadpda/CB 

dp  -  27mi  CB  dpdix 

F  dxdp  =   f  dfi 

do,  =  dadpdadpdix 
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(5.3) 


An  obvious  way  of  interpreting  this  plasma  Is  as  a 
one -dimensional  gas  with  a  stress  tensor 


p.  .  =  p.b.b  .  (5.4) 


■  J 
\  =  ^  (5.5) 

and  with   "molecules"   which  have    Internal   magnetic   moments 

\i  and   Internal   energy  M-B  as  well   as   translatlonal   kinetic 

2  2  11 

energy  p  /2m.      We   have    <p  /2m>  =  ^   ^i    =  p"  P]/^-      This 

substance  Is  a  very  unusual  magnetic  material  whose  mag- 
netic polarization  Is  essentially  unrelated  (except  In 
direction)  to  the  applied  field  [2].  But,  by  reinterpreting 
the  Internal  molecular  energy  as  translatlonal  with  a 
corresponding  two-dlmenslonally  Isotropic  stress  component 

p„(6.  .  -  b.b  .), 
^2^  ij     1  j^' 

p..=p,b.b.+p„(6..-b.b.)  (5.6) 

we  obtain  an  equivalent  set  of  equations  with  a  more  familiar 
appearance . 

The  kinetic  equations  (5-1)  govern  the  gas  motion 
along  a  magnetic  line,  given  Its  constraining  motion.   But, 
since  U  Is  a  macroscopic  velocity,  the  evident  equation  to 
govern  Its  rate  of  change  Is  the  perpendicular  component 
of  the  macroscopic  equation  of  momentum  conservation  (which 
Is  universally  valid  for  microscopic  or  macroscopic  des- 
criptions) , 

b  <    [p(|^  +  u-v^u)  +  dlv  P  -  J  .X  B]  =  0  .     (5-7) 
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Here  u  Is  the  complete  macroscopic  velocity, 

u=U+V,bxV=0,  (5.8) 

and  p  and  P  are  the  sum  of  electron  and  ion  mass  densities 
and  stress  tensors  respectively.   The  macroscopic  moments 
for  electrons  and  ions  individually  are  defined  as 

n  =  p/m  =  CB  /  f  dpd^l 


mV  =  (C^B/n)  /  p  f  dpd^L 


(5.9) 


pe-^  =  I  nkT-j_  =  i  p^  =.  ^B  /  ^P~g^^   f  dpdn 
pe2  =  nkT^  =  P2  =  CB  /  ^lB  f  dpd^i 

Adding  the  electromagnetic  equations 
^  +  curl(B  X  U)  =  0 

div  B  =  0  (5-10) 

curl  B  =  IJ-  J 

to  the  kinetic  equations  (5-1)  and  momentum  constraint  (5-7) > 
we  observe  by  counting  variables  that  the  system  appears  to 
be  complete  with  the  exception  of  the  extra  variable  P   in 
the  particle  Hamiltonian.   We  have  already  committed  our- 
selves to  a  charge  neutral  description  [by  dropping  SE/St 
in  (5-10)  and  by  dropping  the  electric  part  of  the  Lorentz 
force  terms  in  (5-7)  ]•   This  also  determines  />.   In  the 
perpendicular  direction,  the  charge  neutrality  constraint 
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has  already  been  satisfied  by  our  choice  of  a  common 
velocity,  U,  for  Ions  and  electrons.   Taking  n  =  n  and 
V  =  V_  as  Initial  conditions,  we  determine  fi    (or  b-E)  by 
the  constraint  Sv  /St  =  5v_ySt.   These  derivatives  are 
obtained  as  moment  equations  of  the  kinetic  equations  and 
we  find  the  constraint, 

(m,  +  m  )n  2^  =  b-[m  dlv  P   -  m  dlv  P^]  .      (5-11) 
+     —    OS        +      —     —      + 

At  any  Instant,  jb   is  determined  on  each  line  within  an 
Irrelevant  added  constant. 

To  summarize,  the  GCP  equations  consist  of  kinetic 
equations  (5-1)  governing  the  parallel  motion  of  the 
particles  along  a  constraining  magnetic  line,  a  macro- 
scopic momentum  balance  (5- 7)  which  governs  the  transverse 
motion  of  the  constraint,  the  flux  equation  (5. 10),  and 
the  constraint  (5-11)  for  the  potential  f) .      If  we  wish, 
we  can  explicitly  eliminate  the  variables  J,  E,  and  ^   from 
the  equations  in  favor  of  curl  B,  U  X  B,  and  stresses. 
This  leaves  B  as  the  only  evident  electromagnetic  variable. 

A  more  primitive  version  of  the  GCP  theory  is  obtained 
as  a  one -fluid  formulation  Ignoring  the  potential  f) .      This 
is  strictly  correct  only  if  m  =  m_  and  the  ions  and 
electrons  are  given  Identical  initial  conditions.   Other- 
wise a  nontrivial  f)   will  generally  develop  during  the  motion. 
But  the  one-fluid  theory,  setting  /i  =  0  in  the  GC  Hamlltonlan, 
does  serve  the  purpose  of  a  simpler  theory  which  has  many  of 
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the  features  of  the  more  exact  theory. 

We  remark  that  the  charge  neutrality  condition  is  a 
necessity,  not  a  choice.   Within  the  formal  GG  expansion, 
charge  and  current,  q -^  ne ,  J '^  nev,  are  potentially  large 
compared  to  n  and  v.   In  order  to  avoid  infinite  charge 
and  current  as  e  — >  oo,  we  must  have  both  charge  and  current 
neutrality,  n  =  n_  and  u  =  u_,  but  the  resulting  finite 
values  of  q  and  J  are  indeterminate  from  the  kinetic  equa- 
tions.  The  choice  q  =  0,  i-l  J  =  curl  B  results  from 
completing  the  system  by  the  pre-Maxwell  equations.   One 
also  has  the  choice  of  retaining  displacement  current  and 
the  full  Maxwell's  equations.   This  depends  on  a  decision 
with  respect  to  the  speed  of  light. which  Is  Independent  of 
the  GG  parameter  e.   Both  constraints,  E  +  U  x,  B  and  (5. 11) 
for  /5  are  independent  of  the  disposition  of  c.   But  with 
finite  c,  the  momentum  equation  takes  the  full  Lorentz  force 

term,  qE  +  J  X  B.   In  this  expression  we  substitute 

o 
q  =  K^dlv  E  and  M-  J  =  curl  B  -  E/c  ,   remembering  that 

E  =  B  X  U.   In  the  complete  nonlinear  system  of  equations 

the  result  is  considerably  more  complex  than  the  pre-Maxwell 

system.   But  in  a  linearization  about  a  static  equilibrium 

we  lose  qE  as  a  second  order  term,  and  the  only  effect  of 

displacement  current  in  the  entire  system  of  equations  is 

to  increase  the  mass  transverse  to  B.   In  the  term 

b^X  p  Su/Bt  in  the  linearized  momentum  equation,  p   is 

replaced  by 

P^  =  P^  +  B^/m-^c^  =  p^(l+A^c^)   .  (5.12) 
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We  shall  find  that  stability  of  a  GC  plasma  Is  completely 
Insensitive  to  the  value  of  c  (i.e.,  to  the  choice  of  Maxwell 
or  pre-Maxwell  formulation)  In  the  nonlinear  as  well  as  the 
linear  case  even  though  the  equations  of  motion  are  altered. 

We  recall  that  the  CG  orbit  theory  requires  as  a  postu- 
late that  the  parallel  electric  field  be  small.   But  the 
charge  neutrality  condition  (5. 11)  yields  as  a  conclusion 
that  this  Is  so  for  any  GOP  solution  (the  factor  e  has  been 
removed  from  the  potential).   The  more  detailed  physical 
mechanism  involves  high  frequency  plasma  (and  modified  plasma) 
oscillations  which  are  not  explicit  in  the  GO  theory.   Any 
parallel  fields  larger  than  0(l/e)  give  rise  to  large 
velocities,  oppositely  directed  for  ions  and  electrons,  which 
act  to  restore  charge  neutrality  (except  possibly  in  certain 
unstable  situations.  Sec.  10).   More  precisely,  the  neutralizing 
mechanism  involves  oscillations  which  are  presumed  to  exist  at 
only  a  low  energy  level  since  their  presence  would  invalidate 
the  orbit  analysis. 

If  desired,  the  full  momentum  equation  [bracket  in  (5-7)] 
can  be  taken  as  a  governing  equation,  since  the  parallel  momen- 
tum balance  is  a  consequence  of  the  kinetic  equations.   In  this 
case,  the  identification  of  the  parallel  velocity  V  as  a  moment 
of  f  need  only  be  made  initially. 

In  addition  to  conservation  of  momentum  one  can  derive, 
after  considerable  manipulation,  a  relation  Implying  the 
conservation  of  energy.   In  almost  all  analyses  of  the  GCP, 
conservation  of  energy  is  taken  as  an  article  of  faith 
without  verification.   For  the  single  fluid  it  was  verified 
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In  [2],  for  two  fluids  in  [l6] ,  and  for  the  linearized 
equations  (in  which  even  the  formula  for  the  energy  is  quite 
nonintuitlve )  it  has  not  been  previously  verified  at  all. 

Energy  conservation  for  the  GCP  is  a  much  more  subtle 
affair  than  the  equivalent  relation  for  the  Vlasov  equa- 

tions.   In  the  latter  case,  taking  2"(^o^  "*"  ^   /^o^  ^^^  ^^^ 
Maxwell  and  B  /2[i     for  the  pre -Maxwell  electromagnetic 
energy  density,  we  obtain  an  electromagnetic  source  term 
J-E  and  a  balancing  kinetic  source  term  J  x.  B-U.   On  the 
microscopic  level,  a  particle  gains  energy  as  it  moves 
(creating  a  current)  in  the  direction  of  E.   But  a  guiding 
center  particle  always  moves  transverse  to  E;  It  is  a  small 
drift  velocity  of  order  l/e  (not  in  the  GC  theory)  which, 
when  converted  to  a  finite  current  by  multiplication  by 
the  large  parameter  e,  gives  the  conventional  energy  source 
E-J.   This  energy  source  is  missing  in  the  kinetic  part  of 
the  GC  theory.   But  there  is  an  entirely  different  energy 
component,  li-B,  which  varies  when  the  particle  is  carried 
to  a  region  of  different  B.   It  is  a  deeply  hidden  virtue 
of  the  singular  GC  limit  that  energy,  which  one  would 
imagine  is  essentially  concerned  with  higher  order  drifts, 
is  conserved  within  the  lowest  order  theory.   As  a  matter 
of  fact,  if  one  were  to  compute  the  higher  order  drifts  and 
include  the  corresponding  currents  and  energy  sources,  the 
plasma  energy  source  would  be  counted  twice.   Working 
backwards,  one  can  use  the  conservation  of  energy  alone  to 
compute  the  first  order  drifts  directly  from  the  lowest 
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order  GC  results  without  further  recourse  to  the  differ- 
ential equations. 

Specifically,  we  take 

1  ..2 


X  -  I  ^   PU  dx  (5.13) 

k  =-   /(V+  +  £_f_)d^  (5-1^) 

^H-I-^B^dx  (5.15) 


for  the  kinetic,  internal,  and  magnetic  energies  respec- 
tively; 


2 
)_ 

2m 


e  =  g-  +  HB  (5.16) 


1    2 
is  not  the  Hamiltonian;  the  term  p-  mU  must  not  be  kept,  and 

the  terms  t  /5  do  not  contribute  to  the  integral  {^.1^).      For 

the  rate  of  change  of  kinetic  energy,  one  can  show 

^  =  /  U.(J  ;.  B  -  div  P)dx  -  /  pv  ^  (|u^)dx   (5.17) 

where 

P .  .  =  p.,  b  .  b  .  +  Po  ( & .  .  -  b  .  b  . ) 

1.1    ^1  1  .1    ^2^  1,1     1  .1^ 


p^  =  p^  +  pv^  =  CB  /  (p^/m)  f  dpd^l 


For  the  total  plasma  energy  we  obtain 


(5.18) 


'^     {K  +  U)   =  I  \J  ■  J  KB  dx  (5.19) 


dt 
and  Including  the  magnetic  energy, 
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_d_ 
dt 


(/<  +  U+/>p  =  -^|EXB-dS.        (5.20) 


The  boundary  term  takes  this  form  only  If  the  boundary  Is 
fixed.   If  the  plasma  Itself  terminates  at  a  fixed  wall, 
E  =  B  K  U  In  this  Integral.   If  the  plasma  is  surrounded 
by  an  electromagnetic  field,  the  interface  conditions  are 
such  as  to  cancel  in  the  energy  balance,  and  E  is  evaluated 
at  the  fixed  outer  boundary  of  the  entire  plasma  plus  field 
domain  (e.g.,  E,  =  0  at  a  perfectly  conducting  wall). 

The  subdivision  of  plasma  energy  into  K.     and  Li,    is 
peculiar  to  the  GCP  formulation.   In  a  macroscopic  or 
Vlasov  version,  one  would  usually  include  a  component 

r     1  2 

J  p-  pV   in  kinetic  rather  than  Internal  energy. 

The  single  fluid  GCP  equations  linearized  about  an 
arbitrary  static  equilibrium  are 


Bf    r.j  ..     ,  r.  .O. 
O' 


^  =  [H^,f]  +  [H,f^] 


Po  If  +  (^^^  ^^±   =  ^   ^"^1  B  ^  ^o  (5.21) 


5B 

-f-    (J  ui'x  I  r, 
O 


■^  +   curl(B^  K  U)  =  0 


The  perturbed  Hamiltonian  and  stress  tensor  are  most  easily 
writ;;ten  as  rates, 

-^  +  U.  Vp°  =    -p°  div  U  +  Bq   /  m^v^  ||  dvd^i  (5.22) 

-^  +  U-Vp2   =   P2(^o'li-  -   2   div   u)    +   B^   /  mn  II  dvdn 
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These  can  be  integrated  with  respect  to  t  by  introducing 
the  Lagrangian  displacement  x  =  ^{a  ,fi  ,o  ,t)   with 
U  =  D^/Dt  ^  h^/ht, 

H  =  -[i|-^B   -  mv^(b  -ll)  +  iiB^(b  -li  -  div  I) 
'^  o       ^  o  as'  o^  o  OS         ^ 

o  2  2 

p-|  =  -div(^p-,  )  +  B  /  m  V  f  dvd|-L 

(5.25) 
P2  =  -div(^P2)  +  P2(^o"3s"  "  ^^^  ^^  +  ^o  -^  mtifdvd^i 

B/b^  =  -11  -  b^  div  I 

'   O    OS      o 

which  leaves  as  the  only  remaining  differential  equations 

11=  [H„,f]  +  [H,f°] 


(5.24) 

T  +  (div  P)  ,  =  T-r-  curl  B  x  B 
St^ 


Pq  — |-  +  (div  P)j_  =  -J-  curl  B  x  B 


These  linear  equations  have  an  energy  constant,  but  it  is 
much  more  formidable  in  appearance  than  for  the  nonlinear 
system.   The  simplest  computation  is  to  take  the  second 
variation  6  ^  of  the  nonlinear  potential  energy  from  Sec.  11 
This,  expressed  as  a  function  of  the  perturbed  distribution 
6f  ~  f  and  with  |  replacing  U  gives  the  conservation  law 
[cf.  (11.17)] 

/  I  P^U^dx  +  T(e)    +   2/gfdfi 

(5.25) 
-   /   f^(Sf°/Se)"^dn  =   const. 

where  T  is  defined  in  (11.  I8)  and  g  =  H+  [4  •  VB   is  given 
by  (11.11)  or  (5-23). 
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The  energy  constant  for  this  system  does  not  exist 
for  arbitrary  f  when  f   is  not  monotone.   But  if  initially 
f  =  0  where  ^f°/Se  =  0,  this  property  will  maintain  itself 
in  time.   Thus  the  energy  constant  exists  for  a  well- 
defined  special  class.   But  there  is  no  evident  singular 
behavior  of  the  equations  of  motion  where  the  energy 
constant  is  nonexistent. 

Kulsrud  shows  how  to  write  these  equations  in  a  sym- 
metric form  when  f   is  a  single-valued  monotone  function 
of  e .   We  first  illustrate  this  with  the  system  of  ordinary 
differential  equations  with  constant  coefficients 

w  =  Lw  (5.26) 

where  the  matrix  L  is  arbitrary.   It  is  assumed  that  every 
solution  satisfies  the  quadratic  "energy"  relation 

Q  =  p-  wAw  +  wBw  +  p-  wGw  =  const.  (5.27) 

We  may  take  A  and  C  to  be  symmetric.   Since 

Q  =  wAw  +  wL'Bw  +  wBw  +  wCLw  =  0 

for  all  initial  values  w  and  w,  we  have  wBw  =  0  setting 
w  =  0  and  wL'Bw  =  0  from  w  =  0.   This  leaves  w(CL+A)w  =  0 
whence  CL  =  -A  and  (5.26)  can  be  rewritten  as 

Cw  =  -Aw  .  (5.28) 

In  other  words,  multiplication  of  (5.26)  by  the  kinetic 
energy  matrix  C  makes  both  left  and  right  sides  of  the 
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equation  symmetric.   This  reduction  has  significant  con- 
sequences (principally  that  the  spectrum  Is  real)  only 
when  at  least  one  of  the  matrices  A  and  C  Is  definite. 
A  similar  analysis  (but  only  formal  and  Implicit 
since  no  properties  have  been  established  for  the  operators 
analogous  to  the  matrices  L,  A,  B,  C)  reduces  the  GCP 
equations  to  symmetric  form.   To  obtain  a  second  order 
equation  for  f,  Kulsrud  Introduces  even  and  odd  components 
of  f  with  respect  to  v, 

f  =  f  +  f  .  (5-29) 

Since  H  ,  H,  and  f  are  even,  and  the  operator  [H  ,  *] 
o  ^  o 

reverses  parity,  we  have 

(5.30) 

Sf    r..   xTi  ,  r.T  ^O. 

^O^ 


^  =  [H^,f ]  +  [H,f"] 


Only  f  Is  coupled  to  the  remaining  equations  through  the 
even  moments  n,  p.,,  Pp ;  In  the  linearization  the  odd  moment 
V  drops  out.   We  therefore  eliminate  f, 

^  -  [H^,[H^,f]]  =    [H^,[H,f°]]  .  (5.31) 

We  now  have  two  second  order  coupled  equations  In  the 
variables  ^  and  f.   But  they  cannot  be  made  symmetric  by 
only  algebraic  manipulations.   We  illustrate  this  with  the 
simple  example  of  a  perturbation  about  a  uniform  plasma 
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with  constant  B  ,  p-,  ,  Po  •   We  have 
at  ox 


B^  V^  /  ^L  f  dpdii  (5-32) 


Bt       Sx  Be   Sx 


The  "diagonal"  terms  (on  the  left)  coupling  ^  with  itself 
and  f  with  Itself  are  evidently  symmetric;  but  it  is  clear 
that  no  algebraic  manipulation  can  symmetrize  the  cross 
terms  which  Involve  a  third  derivative  of  |  and  a  first 
derivative  of  f .   But  symmetry  can  be  achieved  by  intro- 
ducing either  S^/5x  =   r\   or  /  f  dx  =  g  as  a  new  dependent 
variable,  then  dividing  through  (5.32)  by  v  Sf  /Se.   We 
see  that  the  operator  equivalent  to  C  in  (5. 28)  can  be  quite 
complicated.   Also  the  "kinetic"  energy  in  this  formulation 

is  more  complicated  than  the  natural  kinetic  energy.   The 

2r         _   ^ 

second  variation  6  ^  has  f  =  f  +  f  as  its  argument.   The 

present  interpretation  amounts  to  splitting  5  ^  into  two 
parts,  with  the  f  component  taken  as  kinetic  and  the  f 
component  as  potential.   One  can  expect  the  spectrum  to 
be  real  if  Sf  /Se  <  0,  but  there  will  always  be  a  con- 
tinuous spectrum  extending  to  the  origin  (this  is  easily 
seen  from  the  related  "Rayleigh's  principle"). 

The  symmetric  form  cannot  be  derived  if  f   is  not  a 

2 
function  of  v  ,  even  if  it  is  monotone  on  both  sides. 
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The  reason  is  that  the  Hamlltonian  couples  to  both  f  and  f 
In  (5.30).   In  a  mirror  configuration  f  must  be  even.   In 
a  torus  it  need  not.   Both  cases  are  represented  by  an 
infinite  uniform  plasma. 

We  may  ask  what  the  kinetic  equation  yields  if  we  take 
an  adiabatic  approximation,  Sf/Bt  ~  icof  —^0.   From  (5.24) 
we  derive 

0  =  [H^^f ]  +  [H,f°] 

-  [H  ,f]  -  (^)[H^,H] 


'o'''     -     ^^^Ln^, 


from  which 


?^f° 

^  =  H  l^+  "^(^o)  •  (5.33) 

o 

This  formula  associates  a  specific  function  f  with  any 
given  field  perturbation  ^  through  the  perturbed  Hamiltonian 
H(|).   We  shall  see  in  Sec.  11  that  this  formula  [with  the 
function  it {E    )  specified]  is  related  to  the  pessimistic 
variation. 

With  the  adiabatic  interpretation  of  this  formula,  we 
cannot  insert  this  expression  into  the  magnetic  field  equa- 
tion without  consistently  setting  cd  =  0  throughout,  e.g., 
to  seek  marginal  modes.   But  with  the  pessimistic  inter- 
pretation [which  gives  rise  to  the  same  formula  (5-33)  hy 
a  different  argument]  we  are  not  so  limited.   Inserting 
this  form  of  f  into  the  remainder  of  the  system  (5-24) 
gives  a  symmetric  second  order  system  in  ^  alone.   In  terms 
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of  the  matrix  analogue  (5.28), 

C^w^  =  -(^11^1  +  A^2^2^ 

(5.34) 

we  drop  Wp  from  the  second  equation  and  eliminate  Wp  to 
obtain  the  symmetric  "adiabatic"  system 

^l^'l  "  ^^12''^22^21  "  ^ll^^l  (5-55) 

The  inversion  App  is  explicit  through  the  formula  (5-33) 
in  the  actual  plasma  case.   Monotonicity  of  f   is 
irrelevant  in  deriving  the  adiabatic  equations  since  this 
enters  only  in  the  analogue  of  the  matrix  Cp  which  has  been 
dropped . 

We  can  carry  out  this  procedure  explicitly  for  the 
uniform  plasma.   In  this  case  the  adiabatic  approximation  is 

f  =  -HB^div  ^  |§^  (5.36) 

and  the  adiabatic  equation  is 

^■=  oy   ^-+  c2  V;  div    ^  (5.37) 

5t^  ^    ax^  ^ 


where 


^?  =  p-^  [bX  -  4  -  PiJ 

^2   =   r   ^^>o  +    2P2   -^   B^   /   ^'   If  dpd^x] 


(5.38) 


Po 
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2 
We  shall  come  across  the  stability  criteria  c.  >  0  and 

2 
Cp  >  0  again.   The  anisotropic  wave  equation  (5- 37)  does 

not  properly  describe  any  actual  motion  of  the  plasma  since 
to  is  bounded  away  from  zero,  cf.  (2.2).   But  its  connec- 
tion with  the  pessimistic  variation  establishes  that  when 
Sf  /be   <  0  the  solutions  of  this  wave  equation  are  stable 
or  unstable  together  with  the  actual  motions.   On  the 
other  hand,  the  connection  with  the  adiabatic  variation 
implies  that  even  when  Sf  /be   <{;  0,  the  adiabatic  equation 
(5-37)  [or  its  equivalent  in  the  general  case]  would  give 
information  about  stability  in  any  case  where  one  verified 
that  CD  ~  0  or  k  ~  00  for  the  modes  in  question.   This  adia- 
batic equation  is  interpreted  as  an  actual  equation  of 
motion  by  Andreoletti  [23]. 

The  equation  (5-37)  is  a  macroscopic  one.   We  shall 
find  that  the  macroscopic  GCF  equations  predict  a  more 
complicated  anisotropic  behavior  in  the  uniform  plasma. 
The  reason  is  that  the  GCF  equations  do  describe  actual 
motions,  albeit  in  a  crude  macroscopic  approximation. 

It  is  instructive  to  mention  the  reason  why  the  pessi- 
mistic variation  gives  a  correct  stability  criterion  even 
though  it  does  not  approximate  the  actual  marginal  motion 
(which  is  not  adiabatic).   The  pessimistic  variation 
minimizes  the  numerator  alone  in  a  Rayleigh's  principle  based 
on  the  symmetric  equations;  the  correct  marginal  eigenvalue 
and  eigenmode  arise  from  minimizing  the  entire  fraction  which 
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contains  a  complicated  "kinetic  energy"  contribution  from  f 
in  the  denominator.   But  the  sign  of  the  minimum  value  is 
the  same  in  both  cases. 


6.   The  Macroscopic  Guiding  Center  Fluid 

There  are  standard  expedients  for  converting  a  micro- 
scopic model  to  any  number  of  approximating  macroscopic 
models  [3^].  The  simplest  is  to  write  the  complete  equa- 
tions of  conservation  of  mass,  momentum,  and  energy 
(these  are  exact  but  not  a  closed  system),  and  then  system- 
atically remove  quantities  such  as  the  stress  deviator 
and  heat  flow  which  are  not  represented  in  the  fluid  state; 
another  is  to  add  new  moment  equations.   In  the  GCP  the 
fluid  state  has  two  pressures  (p-,,Pp)  instead  of  one.   After 
writing  the  stress  tensor  as  in  (5.6)  and  dropping  heat 
flow  from  the  equation  of  conservation  of  energy,  there  is 
still  one  equation  missing.   Clearly  this  should  relate  to 
the  conservation  of  magnetic  moment.   We  easily  establish 
that  the  mean  magnetic  moment  per  particle  is  given  by 

<^L>=  kTp/B  =  Pg/nB  =   mp2/pB  .  (6.1) 

Thus  we  are  led  to  adjoin  the  equation  (d/dt  =  B/St  +   u-SZ  ) 

ro  0  =  0  (6-2) 


-  53  - 


Pp/pB  is  constant  on  a  particle  path.   This  equation 
would  also  result  from  the  exact  moment  equation  for 
<  |i  >  after  dropping  a  term  analogous  to  heat  flow  (net 
magnetic  moment  flow  with  respect  to  the  mean  flow). 
As  in  ordinary  fluid  dynamics  or  magneto-fluid 
dynamics,  the  total  energy  equation  is  equivalent  to  an 
entropy  equation.   But  the  entropy  in  this  case  turns  out 
to  be 

log  T]  ^  ^   log  p^  +  log  P2  -  f  log  P         (6.3) 

Instead  of  the  usual 

log  T]  -- c^  log  P  -  Cp  log  p  (6.4) 

This  can  be  obtained  by  manipulation  of  the  energy  equation, 
and  it  also  arises  as  /  f  log  f  for  an  equilibrium  f  which 
is  a  product  of  two  Maxwelllans ,  one  degree  of  freedom  at 
T-,  and  two  at  Tp . 

Our  final  system  is  therefore 

d£ 


dt 


+  p  dlv  u  =  0 


p  ^  +  dlv  P  =  7^  curl  B  X  B  (6-5) 

o 

—  -  B-  Vu  +  B  dlv  u  =  0  , 

plus  two  energy  or  entropy  equations.   For  the  total  energy 
per  mass 
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e  =  ^  RT-^  +  RTg  =  Pi/2p  +  Pg/p  (6.6) 

we  may  take  the  conservation  equation 

or  the  equivalent  entropy  equation  [cf.  (6.3)] 

^  (PlP^P^)  =  0  .  (6.8) 


For  magnetic  moment,  we  may  take  the  conservation  equation 

r\         P?      Pp 

^  (-^)  +  ^  dlv  u  =  0  (6.9) 


or  the  particle  path  Invariant 

^  (p-|)  =  °  ■  (6.10) 

The  two  particle  path  Invariants  can  be  chosen  as  any  two 
of  the  following: 

P2/PB  ,    p^B^/p^  ,  PlP2/p^  >    p|/P3_b5  .        (6.11) 

The  stability  of  these  GCF  equations  for  small  per- 
turbations about  a  constant  state  was  first  studied  by 
analyzing  the  growth  of  plane  waves  (normal  modes)  [37], 
[38],  Instability  Is  found  whenever  the  unperturbed 
constant  state  Is  too  anisotropic.   But  a  more  careful 
analysis  of  these  results  showed  that  they  are  techni- 
cally not  Instabilities;  they  describe  breakdown  of  the 
equations  of  motion  [11].  In  the  unstable  regime  the  equations 
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become  elliptic  and  the  initial  value  problem  is  not  well 
posed.   We  choose  to  interpret  this  breakdown  as  an 
instability,  but  this  is  a  physical  interpretation  and 
not  a  mathematical  deduction. 

Since  the  GCF  system  (as  distinguished  from  the  GCP) 
consists  of  partial  differential  equations,  we  can  investi- 
gate well-posedness  in  complete  generality  by  computing 
the  characteristics.   The  most  convenient  representation  is 
in  terms  of  the  normal  speed  locus.   Taking  c  as  the 
characteristic  speed  relative  to  the  mean  motion  u,  the 
characteristic  equations  for  (6.5),  (6.8),  (6. 10)  become 


-c6p   +  p6u    =0 

-pc6u   +   n-5P   =  -^    (n  >^   oB)  X  B 
^o 

-c6B    -   B    6u   +  B5u^    =0  (6.12) 

n  n  '  ' 

^Pl  ^Pp  bo 

— i.  +  2  — -  -5^=0 

P^  P2  P 

6p. 


;;2   _   6£  _   &B 


Pq         P  B 


0 


where  n  is  the  normal  to  the  wavefront.   The  normal  speed  c 
is  the  eigenvalue  of  this  homogeneous  linear  system. 
Ignoring  particle  path  invariants,  c  =  0,  we  find  the 
modified  Alfven  or  transverse  wave 

3^  o  -    o   o  _     Po  ~  P 


2 
c 


=  ^  (^  +  Pp  -  P.  )cos^e  =  a2cos20[1  +  -^ i-]     (6.13) 
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and  the  modified  compressive  (or  slow  plus  fast)  wave  as 
the  solution  of  the  biquadratic 

2  1      2 
p  c   -  Xpc   +  Y  =  0 

X  =  B^/H^  +  2p2  +  (2p^  -  P2)cos^0 

Y  =  cos^e(Y^sln^e  +  Ygcos^e)  (6.11) 

Y^  =  3P;l(^^/%  +  ^^2  "  P2/5P1) 
Y2  =  3P;j_(B^/m.^  +  P2  -  P^) 

Here  cos  9   =  b-n  Is  the  angle  between  B  and  the  wave  normal. 

Although  the  normal  speed  c  as  a  function  of  6   gives 
the  characteristic  cone  only  after  some  manipulation,  loss 
of  hyperbolicity  follows  Immediately  if  c  is  not  real  for 
all  9.      Prom  the  modified  Alfv^n  wave  (which,  as  in  MH, 
represents  strictly  one -dimensional  propagation  along  field 
lines  [  39] )  we  obtain  the  necessary  condition  for  "stability" 

B  /[i^   +  P2  >  Pi  •  (6.15) 

2 
Since  X  and  the  discriminant  X   -  1Y  are  always  positive, 

the  single  condition  Y  >  0  is  necessary  and  sufficient  for 

2 
c  >  0  in  (6.14).   Since  Y2  >  0  is  Implied  by  (6.I5),  we 

need  adjoin  the  single  additional  condition  Y-,  >  0, 

B^/[x^   +  2p2  >  P2/3P1  .  (6.16) 

The  two  conditions  (6.I5)  and  (6.16)  are  necessary  for 
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the  system  to  be  hyperbolic;  to  show  that  they  are  also 
sufficient  requires  verification  that  the  fast  cone 
(giving  the  domain  of  dependence)  is  convex.   The  Inequality 
(6.16)  is  a  consequence  of  (6.15)  and  can  be  dropped  if 
desired  when 

-2  <  ^  (3 +y2T)  -  3.79  .  (6.17) 

P-|       C- 

These  conditions  are  necessary  for  stability  in  an 
arbitrary  nonlinear,  nonuniform,  nonsteady  flow.   The  equa- 
tions are  Improper  if  hyperboliclty  is  violated  anywhere 
in  the  flow.   The  enormous  generality  of  this  as  a  stability 
criterion  is  a  result  of  the  disastrous  nature  of  the 
instability.   With  any  more  conventional  interpretation  of 
stability  it  is  not  clear  how  to  even  define  stability  in 
nonsteady  flows.   There  is  a  further  deduction  that  follows 
from  our  Interpretation  of  non-posedness  as  instability, 
viz.,  that  such  instabilities  are  micro-instabilities;  we 
shall  return  to  this  point  in  the  discussion  of  the  GCP . 

A  similar  situation  occurs  in  the  common  Raylelgh- 
Taylor  instability,  e.g.,  of  an  inverted  glass  of  water. 
Here  too,  the  linear  equations  of  motion  are  not  well 
posed.   The  difficulty  is  resolved  either  by  examining  the 
nonlinear  equations  or  by  inserting  surface  tension  in  the 
linear  equations.   For  the  GCF,  we  have  seen  that  the 
difficulty  is  not  removed  in  the  nonlinear  version.   We 
shall  see  that  it  also  appears  (although  quantitatively 
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modified)  in  the  microscopic  GCP.   It  is  apparently  a 
consequence  of  the  singular  GO  limit,  but  this  remains 
only  a  plausible  conjecture  until  better  existence  theorems 
are  established  in  the  finite  gyro  radius  and  Debye  radius 
Vlasov  formulations. 

We  conclude  this  very  brief  review  of  GCF  theory  with 
the  remark  that  the  classical  scalar  pressure  MH  results 
do  not  follow  from  setting  p.,  =  Pp  in  these  equations. 
The  Alfven  wave  does  formally  reduce  to  the  MH  result,  but 
this  is  not  true  of  the  compressive  waves.   The  reason  is 
that  an  isotropic  initial  state  will  not  remain  so,  even 
for  the  linear  equations.   The  constraint  (6.2)  is  quite 
distinct  from  anything  which  appears  in  the  MH  theory.   Two 
simple  cases  which  are  easily  compared  with  the  MH  theory 
are  propagation  perpendicular  to  the  field  (cos  9=0)  and 
parallel  (cos  9=1).   In  the  first  case  we  have 

c^  =  0  ,  c^  =  B^/n^p  +  2P2/P   .         (6.18) 

The  zero  root  agrees  with  the  slow  wave,  and  the  finite 
root  corresponds  to  the  fast  wave  in  a  two-dimensional  gas 
(y  -   2).   For  parallel  propagation  we  have 

c  =  B  /n^p  +  (P2  -  p^)/p  ,  c  =  3p-j_/p   .     (6.19) 

The  modified  Alfven  root  is  repeated,  as  in  the  MH  case, 
and  there  is  an  acoustic  speed  corresponding  to  a  one- 
dimensional  gas,  Y  =  3-   Thus  the  anlsotropy  Introduced  by 
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the  magnetic  field  in  MH  is  compounded  by  the  anisotropy 
of  the  gas  in  the  GCF;  even  with  p,  =  Pp ,  the  various 
degrees  of  freedom  are  uncoupled  and  have  individual  gas 
constants  7=2,3  instead  of  a  combined  y  =  5/3-   It  is 
also  an  interesting  fact  that  the  transverse  (Alfve'n) 
speed  does  not  always  separate  the  slow  and  fast  speeds 
as  it  does  in  MH;  other  conclusions  await  a  detailed 
calculation  of  the  GCF  characteristic  cones. 


7.   The  Infinite  Uniform  Plasma 

The  linearized  motion  about  a  constant  state  is  well 
suited  for  analysis  by  Fourier  and  Laplace  techniques  [17], 
[^0],  [3],  [18],  but  has  been  evaluated  only  recently  by 
Kadish  using  the  strict  criterion  of  boundedness  of  solu- 
tions [19]-  First  of  all,  it  is  found  that  "instability"  is 
equivalent  to  non-well  posedness  for  every  discrete  normal 
mode.   This  is  easily  seen  qualitatively.   Since  there  is  no 
length  or  time  parameter  in  the  equations  of  motion,  the 
dispersion  relation  is  a  function  of  (oo/k-,  ,  a/kp).   If 
this  ratio  takes  any  complex  or  imaginary  value,  the  growth 
rate,  Im  go,  becomes  arbitrarily  large  for  small  wavelength 
disturbances.   On  the  other  hand,  for  any  stable  set  of 
parameters  (described  below)  the  existence  and  uniqueness 
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of  a  solution  to  the  initial  value  problems  is  established. 

The  constant  unperturbed  state  is  defined  by  an 
equilibrium  distribution  function  f  (e,l-L)  [  e  =  p-  mv   +  [i-B], 
and  a  constant  vector  B  ,  or  by  two  distributions  f    {e,\x); 
in  the  uniform  plasma  there  is  no  alternative  to  j?!)  =  0 
in  equilibrium.   There  is  a  continuous  spectrum  on  the 
entire  real  co-axis.   This  continuum  gives  a  non-exponential 
decaying  contribution  and  does  not  give  rise  to  any 
instability  if  Sf  /'be   is  Holder  continuous.   For  such 
smooth  functions  f  ,  Kadish  establishes  the  following 
necessary  conditions  for  stability  for  the  single  fluid 

b^/Hq  +  P2  >  Pi  (7.1) 

B^/[i     +   2p2   >    -mB^   f  tx^(SfySe)dvd^x  (7-2) 

The  first  is  the  same  as  in  the  GGF  (6.15),  and  the 
second  becomes  comparable  to  the  second  GGF  condition 
(6.16)  when  we  remark  on  the  inequality  (for  Sf  /he   <  0 
only) 

-mB^  /  [i^(afySe)dvdii  >   P2/3p^      .  (7-3) 

The  two  conditions  (7.I)  and  (7-2)  become  sufficient  as  well 
as  necessary  for  stability  if  f   is  monotone  in  e  (bell 
shaped  in  v).   But  there  are  definite  counterexamples  that 
show  lack  of  sufficiency  when  f   is  not  monotone. 

In  the  two-fluid  case  the  second  condition  (7.2)  is 
replaced  by 
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where  C^  is  a  certain  complicated  combination  of  moments  of 
Bf°/5e  [see  (10. l6)  and  (10.2^)],   C^  also  satisfies  the 
Inequality  (7o)  for  the  total  pressures  Pp  =  Pp  +  p,J, 
P-,  =  P-|  +  p.  If  both  distributions  are  monotone. 

The  continuous  spectrum  can  give  rise  to  unstable 
solutions  If  Sf  /Se  Is  not  sufficiently  smooth;  the  growth 
is  not  exponential  for  such  solutions.   With  f  not  monotone 
In  e,    the  situation  Is  not  completely  known.   But  explicit 
examples  can  be  given  of  complex  (not  pure  Imaginary)  poles 
for  CD;  and  examples  can  be  given  of  such  poles  making  their 
appearance  at  some  point  removed  from  cd  =   0-    It  has 
not  been  established  whether  there  Is  any  continuum  In 
addition  to  the  one  on  the  real  oj-axls;  but  there  Is  no 
continuum  In  the  upper (exponentially  unstable)  half  plane. 

Writing  Sf  /Se  =  (mv)    Bf  /Sv  we  see  that  the 
criterion  (7'2)  depends  sensitively  on  the  behavior  of 
f  (v)  near  v  =  0.   This  observation  will  take  on  much  more 
significance  In  a  nonuniform  contained  equilibrium  (See.  15)- 

With  the  variational  definition  of  stability  to  be  des- 
cribed, the  results  turn  out  to  be  the  same  as  for  boundedness 
of  solutions  except  for  the  monotonlclty  criterion  and  the 
contribution  from  the  continuum.   No  instability  from  the 
continuum  can  be  found  varlationally  because  the  continuum 
is  on  the  real  co-axis.   Also,  the  variational  analysis 
categorizes  all  non-monotone  distributions  as  unstable. 


A.  Kadlsh,  private  communication 
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Thus  the  necessary  and  sufficient  condition  for  variational 
stability  of  the  Infinite  uniform  plasma  is  [29] 

Sf7Se  <  0 

B^M-Q  +  P2  -  Pi  >  0  (7.5) 

B^^L^  +  2p2  -  B^C^  >  0 

This  is  also  the  exact  necessary  and  sufficient  condition 
for  variational  stability  of  a  nonuniform  plasma  with 
unidirectional  B  and  transverse  variation  subject  to 
Pp  +  B  /2ii   =  const . 


8.   Equilibrium- -MH 

We  are  concerned  primarily  with  GO  equilibria  but 
present  a  brief  summary  of  the  scalar  pressure  theory  for 

•X- 

comparison.    The  equations  are 

VP  =  J  X-  B 
div  B  =  0  (8-1) 

|j.  J  =  curl  B 

o 


■X- 

For  a  more  complete  account  of  the  material  in  this 
section,  including  toroidal  geometries  which  we  do  not  dis- 
cuss here,  see  [4l].  The  natural  boundary  conditions  and 
potentials  1^   and  q  are  discussed  in  [42]. 
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By  inspection  we  see  that  the  magnetic  lines  and  the  current 
lines  lie  In  constant  pressure  surfaces.   From  dlv  B  =  0  and 
J-Vp  =  0  we  can  conclude  that  the  magnetic  field  Is  a  two- 
dimensional  weighted  harmonic  vector  field  In  each  pressure 
surface . 

B=  |Vp|nxV"^ 

(8.2) 
B  =  V'j^  =  Vj25  -   n{n-Vj6)    . 

(without  the  weight  |^p|,  the  scalars  ^   and  ud   would  be  exactly- 
conjugate  surface  harmonics).   Given  the  set  of  pressure  sur- 
faces, the  vector  field  B  within  each  surface  is  subject  to  an 
elliptic  equation;  thus  B  is  determined  either  by  boundary 
conditions  or,  on  a  topologlcally  complex  surface  (e.g.,  a 
cylinder)  by  periods. 

The  real  characteristics  of  the  system  (8.1)  are  the  mag- 
netic lines  counted  twice,  and  there  is,  in  addition,  a  second 
order  elliptic  characteristic  cone.   Identification  of  the 
characteristics  allows  us  to  guess  at  plausible  well-posed 
problems.   These  guesses  are  refined  by  the  variational  analysis 
to  follow  and  by  the  recognition  of  B  as  a  surface  harmonic 
within  a  pressure  surface.   In  special  cases  where  the  real 
characteristics  are  explicitly  integrated  out,  we  can  refer  to 
actual  existence  theorems  for  purely  elliptic  equations. 

Consider  the  tubular  domain  of  Fig.  1.   The  boundary-  con- 
ditions B   =  given  (e.g.,  B  >  0  on  S ' ,  B  <  0  on  S" ,  B  =0 
n*^     ^^n  n  n 

on  S  )  are  appropriate  to  the  elliptic  cone.   The  real  char- 
acteristics require  two  initial  values  at  one  end  of  the  tube 
or  one  at  each  end.   We  may  specify  p  and  J   (the  normal  current 
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flow)  on  S '  or  specify  p  at  both  ends.   In  the  latter  case, 
the  two  p  specifications  must  be  compatible  with  the  B 
specification 


I   B^dS'  =  J   B^dS"  (8.3) 


P<Po        P<Po 

and,  in  addition,  a  period  of  B  must  be  given  on  each 
p-surface  if  it  is  closed  (this  is  not  Immediately  evident 
by  counting  characteristics). 

A  variational  formulation  allows  more  precise  evalua- 
tion of  permissible  boundary  conditions  than  is  given  by 
the  mere  counting  of  characteristics.   If  we  take 


(8.4) 


div  B  =  0 
B-Vp  =  0 
as  admissibility  conditions,  variation  of 

|(B,p)  =  /(B^/2n^  -  p)dx  (8.5) 

gives  the  first  line  of  (8.1)  as  an  Euler-Lagrange  equation. 

The  admissibility  conditions  (8. ■4)  can  be  more  con- 
veniently parametrized  by  taking 

B  =  Va  X  VP 

(8.6) 
p  =  p(a,p) 

where  p(a,p)  is  a  given  function.   This  is  essentially  a 
new  formulation  with  the  variational  function 

J(a,p)  =  [[2^  iVa  X  VPI^  -  p(a,P)]  dx        (8.7) 
■^         o 

An  interior  variation  of  (8.7)  again  gives  Vp  =  J  x  B. 

The  boundary  variation  vanishes  if  a  and  p  are  fixed  on 

S'  and  S"  and  B  =  0  on  S  .   Note  that  specification  of 

no  ^ 

a  and  p  is  stronger  than  just  giving  B  .   For  example,  if 
the  p-surfaces  are  tubular,  fixing  the  ends  of  the  magnetic 
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lines  determines  the  "twist".  I.e.,  the  period  of  the 
harmonic  vector  B  In  the  pressure  surface. 

Relaxation  of  a  boundary  condition  will  give  rise  to 
a  corresponding  natural  boundary  condition.   For  example, 
suppose  the  pressure  contours  assigned  on  S'  are  concen- 
tric curves  and  take  a  -   const,  on  these  contours, 
p  =  p(a).   We  fix  a  on  S '  but  drop  the  boundary  condition 
on  p.   The  natural  boundary  condition  turns  out  to  be  zero 
twist, 

§   B-dx  =  0  (8.8) 

P 

on  each  pressure  contour. 

The  more  drastic  relaxation  of  both  a  and  p  on  S ' , 
keeping  B  fixed,  leads  to  the  natural  boundary  condition 
§   B-dx  =  0  on  any  contour  within  S';  In  other  words, 

J   =  0  on  S'  .  (8.9) 

This  result  gives  the  connection  between  (8.5)  and 
the  classical  Dlrlchlet  principle.   If  we  choose  the 
function  p(a,p)  =  0,  the  variational  equation 

curl  B  X  B  =  0  (8.10) 


describes  a  force-free  field  or  pressureless  plasma.   The 

r    r  1   2 
variational  function  Is  Dlrlchlet 's  Integral,  £  =  J  2"  B  > 

but  the  boundary  conditions,  a  and  (3  fixed,  are  more 

restrictive  than  the  proper  boundary  condition,  B  =  given, 

for  a  Dlrlchlet  problem.   The  equations  (8. 10)  can  be 


66 


written 

curl  B  -  aB  (8.11) 

where  a(a,p)  is  constant  on  each  line.   The  natural  boundary 
condition  (8.9) j  obtained  by  relaxing  the  boundary  condi- 
tions from  (a,p)  =  given  to  B  =  given,  yields  a  =  0;  thus 
curl  B  =  0,  and  B  is  harmonic  as  in  the  Dirichlet  problem. 
If  a  and  p  are  relaxed  at  both  ends,  S'  and  S" ,  we  of 

course  obtain  J  =  0  at  both  ends.   A  somewhat  weaker 
n 

relaxation  of  the  original  boundary  condition  is  the 
interchange .   This  allows  (a,p)  to  vary  arbitrarily 
(subject  to  fixed  B  )  at  one  end,  but  forces  an  identical 
variation  of  (a,p)  at  the  other  end.   To  be  precise,  we 
note  that  any  admissible  (a',P')  carrying  the  same  B  as 
(a,p)  is  obtained  by  an  incompressible  mapping. 

We  require  the  transformation  (a,p)  — >-  (a',f3')  to  be  the 
same  on  S'  and  S" ,  but  it  is  otherwise  arbitrary.   For  a 
given  admissible  field  B,  any  transformation  (8.12)  is 
trivial  and  leaves  B(x)  unchanged.   What  it  does  is  alter 
the  relative  assignment  of  p-values  to  magnetic  lines. 
Thus  an  interchange  does  not  alter  B  but  maps  p(a,p)  incom- 
pressibly  in  the  (a,p)  plane.  The  resulting  natural  boundary 
condition  is  that  the  current  input  and  output  are  equal 
at  the  two  ends  of  each  magnetic  tube,  specifically 
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where 


j'  +  j"  =  0  (8.15) 


j  =  J  /B  (8.1^) 


is  the  current  density  with  respect  to  flux, 

J^dadp  ^   J^dS  (8.15) 

Just  as  B  can  be  parametrized  by  B  =  Vp  "<  Sl'^i 
J  can  also  be  written 

J  =  ^C  ^  ^P  (8.16) 

where  J-dS  =  di^dp.   The  stream  function  ^  (which  is 
constant  on  a  J-line)  is  closely  related  to  the  Indefinite 
integral 


o 


S 

B(a,p,sj 


q(«^P'S)  =  f  ^,%      .  (8.17) 


1 


s 

along  a  magnetic  line;  the  lower  limit  is  taken  on  S ' .   It 
is  easily  shown  that  (^  differs  from  q  by  a  function  of  a 
and  p  alone,  i.e.,  they  differ  by  a  constant  on  each  mag- 
netic line  [42].  With  closed  p-contours,  ^  can  be  multivalued; 
the  "no-twist"  condition  (8.8)  makes  t,    single-valued.   If 
J  =  0  on  S '  we  can  take  t   -   (\   and  write 

J  =  S/q  X  \7p  (8.18) 

In  general,  ^  differs  from  q  by  the  initial  value  of  C  at 
one  end;  this  is  closely  related  to  the  initial  J   on  S'. 
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If  we  have  J  =  0  on  S '  and  S" ,  we  conclude  that 
n 

s" 

q(a,p)  =  q(a,p,s")  =  j   ds/B  (8.19) 

s  ' 

is  a  function  of  p(a,p);  ( \7q  x  Vp)  "^S  =  J-dS  =  0.   Thus  a 
special  equilibrium  with  J  =  0  at  both  ends  Is  character- 
ized by  p  contours  which  coincide  with  q  contours.   This  is 
also  true  of  the  more  general  interchange  equilibria  in 
which  the  (a,p)  specification  is  only  partially  relaxed  as 
described  above.   Allowing  a  and  p  to  be  free  at  both  ends 
requires  p  to  readjust  so  that  it  corresponds  to  q.   The 
natural  boundary  condition  j   +  j"  =  0  is  equivalent  to 
the  global  condition  p  ^  q. 

The  boundary  condition  that  (a,p)  be  fixed  at  an  end 
wall  S'  has  the  physical  Interpretation  of  a  perfectly 
conducting  wall;  the  perfectly  conducting  boundary  condi- 
tion E,  =  0  implies  no  tangential  motion  u,  =  0  in  the 
presence  of  B  /^  0.   An  insulating  end  wall  will  allow 
the  ends  of  the  magnetic  lines  to  move.   The  two  physical 
facts  associated  with  an  Insulator,  motion  of  magnetic 
lines  and  inability  to  draw  current,  are  seen  to  be  related 
mathematically  as  conjugate  boundary  conditions.   The  con- 
straint B  =  fixed  is  a  mathematical  convenience  to  isolate 
the  given  domain  from  the  external  magnetic  field  which 
can  then  be  considered  to  be  unaffected  by  any  internal 
motions.   It  could  be  realized  by  a  perfectly  conducting 
wall  faced  by  a  thin  insulator  in  contact  with  the  plasma. 
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The  Interchange  condition  (as  distinguished  from  both 
ends  moving  independently)  has  no  apparent  physical 
realization.   It  is  related  to  the  mathematical  realiza- 
tion of  a  toroidal  system  as  an  open  one  with  identified 
ends.   The  condition  that  q  be  a  function  of  p  is  a  restric- 
tion on  the  class  of  open-ended  equilibria,  but  it  is  a 
property  of  all  toroidal  equilibria  [-^3]  as  an  alternative 
statement  of  the  fact  that  J   is  single  valued. 

In  a  two-dimensional  or  axially  symmetric  problem, 
the  real  characteristics  can  be  integrated  out  explicitly, 
leaving  a  purely  elliptic  problem  with  known  theory.   We 
Introduce  the  stream  function  if   for  B  in  an  axially 
symmetric  geometry  and  consider  p  =  p(''/^)  as  a  given  function. 
The  equation  for  i/   is 

^  A%  =  -r^p'(-^) 

(8.20) 
A*  =  S^/Sz^  +  ^^/br^    -  p  S/Sr  . 

For  sufficiently  well-behaved  functions  p{if)    and  a  domain 
which  is  not  too  large,  solutions  of  this  nonlinear  elliptic 
equation  exist  and  are  unique.   The  more  general  axially 
symmetric  problem  with  Bn  ^   0   exhibits  a  similar  reduction. 
Writing 

B  =  ^  n  X.  ^^  +  nBg  (8.21) 
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where  n  Is  the  azlmuthal  unit  vector,  we  discover  that 

2  2 
r  Bq  Is  constant  on  flux  tubes.   Setting 


we  find 


-^  A*^  =   -r^p'iif)    -   g'(^)  (8.23) 

o 

with  the  same  theory  as  the  special  case  (8.20).   The 

function  g('^)  Is  closely  related  to  the  boundary  value 

of  J  ;  thus  the  two  functions  p(i/')  and  giif)    amount  to  a 

direct  Integration  of  the  two  boundary  conditions  p  =  given, 

J  =  given  along  the  real  characteristics.   With  this 
n   •=         to 

degree  of  symmetry,  p  given  at  one  end  determines  p  at  the 
other  end,  so  this  alternative  formulation  is  lost. 

The  differential  equation  (8.23)  can  frequently  be 
solved  by  Iteration.   This  suggests  a  procedure  for  solving 
problems  in  a  general  geometry  [^1].  Consider  the  tubular 
domain.  Fig.  1,  with  B  given  and  p  and  J  given  at  one  end. 
Taking  an  initial  approximation  to  B,  which  satisfies 
dlv  B  =  0  and  the  boundary  condition  B  ,  we  find  p  in  the 
entire  domain  by  extending  the  values  along  B-lines.   Next 
we  compute  J  from  the  system 

Vp  =  J  X  B 

(8.24) 
dlv  J  =  0 

taking  p  and  B  as  given.   To  do  this  most  directly  we  use 
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(8.16),    J   =3  VCxVp,   where 

s 

C   -   C    +   q  =   C(a,P,s')   +  I     ds/B    ;  (8.25) 

s' 

^ '  Is  determined  by  the  boundary  value  of  J  .   Having 

obtained  J,  we  compute  the  next  approximation  to  B  from 

curl  B  =  J/ml 

(8.26) 
div  B  =  0 

and  continue. 

An  iteration  process  should  preferably  have  the 
property  of  smoothing  with  regard  to  derivatives.   The 
iteration  described  does  not  seem  to  be  smoothing,  although 
this  is  so  in  axial  symmetry.   The  construction  of  p(a,p) 
is  such  as  to  plausibly  imply  that  ^p  is  comparable  to 
B  (=  Va  x^p).   But  VC  also  involves  a  derivative  of  B 
(in  the  direction  normal  to  b).   Thus  inversion  of  (8.26), 
which  gains  one  derivative,  seems  to  preserve  but  not  improve 
the  level  of  differentiability  in  a  complete  iteration. 

A  possibly  superior  form  for  the  iteration  would  be 
to  solve  for  B  from 

^  curl  B  -  cB  =  ^^ 
°  ^  (8.27) 

div  B  =  0 

assuming  that  the  right  side  is  given  and  both  B  and  a  are 
to  be  found  on  the  left.   This  is  a  nonstandard  differential 
system  which  seems  to  gain  a  derivative  for  B  over  Vp  and  B 
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on  the  right.   This  differs  from  the  previous  Iteration  In  that 
the  qualifications  parallel  and  perpendicular  refer  to  the  new 
B  rather  than  to  the  old. 

There  are  a  number  of  plausible  Iteration  processes  based 
on  the  property  that  B  satisfies  an  elliptic  two-dimensional 
equation  in  any  flux  surface;  but  in  the  absence  of  a  conver- 
gence proof  there  is  little  point  in  describing  these.   In 
summary,  we  can  refer  to  actual  theorems  on  nonlinear  elliptic 
partial  differential  equations  to  prove  existence  in  the  axially 
symmetric  (and  the  two-dimensional)  problem;  but  we  can  only 
rely  on  plausibility  arguments  for  the  general  case. 


9-   Equillbrlum--GCF 

For  the  GCF  we  have  the  equilibrium  equations 

dlv  P  -  J  K  B 

div  B  =  0  (9.1) 

\i  J   =   curl  B 
The  parallel  and  perpendicular  components  of  the  pressure 
balance  take  the  form 

S^'tPl-Ps'li  (9.2) 

(Pl  -  P2)k  +  y^Pg  =  J  K  B  (9-3) 

where  k  =  Sb/5s  is  the  magnetic  line  curvature  and 

ViPp  ^  (Bx\7p2xB)/B  .   Taking  B  as  the  independent  variable 

along  the  line  instead  of  arclength,  (9-2)  takes  the  form 
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?^  Pi     p? 


Expanding 


curl  B  K  B  =  B-VB  -  V  ^  B^ 


(9.5) 


=  KB^  -  V^  I  B^ 
and  Introducing 


P*  =  Pp  +  B^/2n   ,  (9.6) 


the  perpendicular  equation  (9-3)  becomes 


V^p^  -  k(b^/h^  +  P2  -  P^)  .  (9-7) 


Since  this  system  Involves  two  pressures  Instead  of 
only  one  as  In  the  MH  version,  we  will  expect  to  supply 
an  additional  relation.   This  system  reduces  exactly  to  the 
MH  system  If  we  take  the  additional  relation  as  the  state- 
ment that  p.  Is  constant  on  a  magnetic  line.   Prom  (9.2) 
we  Immediately  conclude  that  p^  =  Pp j  and  everything  else 
follows.   If  Instead  we  assume  that  Pp  Is  constant  on  a 
line,  we  conclude  from  (9.'^)  only  that 

p^  =  P2  +  c(a,p)B  (9.8) 

which  Is  slightly  more  general.   In  the  next  section  we 
shall  see  that  specification  of  a  distribution  function,- 
as  Is  appropriate  to  a  GGP  equilibrium,  will  assign  both 
P-,  and  Pp  as  functions  of  B  along  each  line,  but  In  such  a 
way  that  they  are  automatically  compatible  with  (9-2).   Thus 
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we  can  consider  that  either  p-,(B,a,p)  or  PpCBjQjpi)  is 
given  arbitrarily  and  we  can  then  find  the  other  from 
(9.4).   [in  an  equilibrium  contained  by  mirrors,  the  inte- 
gration constant  in  p.,  is  determined  by  the  fact  that  p., 
and  Pp  vanish  simultaneously  at  the  mirror.]   ¥e  shall  show 
that  a  solution  can  be  expected  to  exist  given  Pp(B,a,p) 
and  boundary  conditions  on  (a,p).   One  special  case  (MH) 
is  if  the  given  function  Pp(a,p)  is  independent  of  B.   Two 
other  interesting  special  cases  which  we  shall  consider  are 
Pp(B)  independent of  (a,p)  and  p^(a,p)  independent  of  B. 
Another  way  of  completing  the  system  is  to  adjoin 
one  of  the  particle  path  invariants  appropriate  to  a  motion. 
Since  p  does  not  appear  in  (9-1)  we  adjoin  the  relation 
(6.11) 

p|/P3_b5  =  a(x)  (9.9) 

where  a(x)  is  an  arbitrary  given  function.   We  can  incor- 
porate both  types  of  supplementary  relation  in  the  very 
general  form, 

g(p^,P2,B,x)  =  0   .  (9.10) 

An  elementary  computation  of  characteristics  shows 
that,  in  addition  to  real  characteristics  along  the  field 
lines,  we  obtain  an  elliptic  characteristic  cone  provided 
that 

B^/H^  +  P2  -  Pi  >  0  (9.11) 
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and 

|^[!^lf^-(P2-pi)|§7-Bii>°-     (9-12) 

The  first  relation  has  already  been  fo-und  as  a  condition  for 
the  Initial  value  problem  to  make  sense.   Although  not 
compelling,  it  Is  natural  to  require  that  a  static  equilib- 
rium formulation  be  elliptic  (e.g.,  to  be  able  to  impose 
B  as  a  boundary  condition).   We  therefore  choose  to 
Interpret  the  equilibrium  problem  as  being  well-posed  only 
when  it  is  elliptic  (that  is  in  addition  to  the  twice  covered 
real  magnetic  lines),  and  we  further  characterize  non- 
ellipticlty  as  instability. 

First,  taking  the  particle  path  invariant  (9-9)? 
5  log  p,^  -  log  P-]_  -  5  log  B  -  log  a  =  g(p-L^P2^B,x) ,  we  have 

^g-_^  ^&     -  J_       ^S-_5  (Q  1^) 

3^  "   P^  '  'S^  ~  v^   '  'Sb  b         i^-^^; 

and    (9-12)    takes   the   form 

bVHq  +   2p2   >   Pg/^P-^  (9.1^) 

which  is  exactly  (6.l6).   In  other  words,  with  the  supple- 
mentary choice  (9-9) J  the  same  conditions  are  found  for  an 
elliptic  equilibrium  as  for  a  hyperbolic  motion.   This  is 
not  surprising. 

Next  take  the  supplementary  condition 

p,  =  P(B,a,p)  (9.15) 
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we  have 


and  the  stability  condition  Is 


B^/^^o  +  B  II  >  0  (9.17) 


which  can  be  rewritten  as 


>  0  .  (9.18) 


This  "stability"  condition  Is  completely  different  from 
(9.1^)  and  Is  apparently  unrelated  to  stability  of  a  motion; 
but  the  two  concepts  will  be  tied  together  In  the  next 
section. 

The  condition  (9.I8)  can  be  obtained  for  a  larger 
class  of  supplementary  relations,  viz.,  if  in  (9. 10)  we 
assume  that  g  is  not  a  general  function  of  x  but  only  a 
function  of  (a,p), 

g(p^,P2,B,a,p)  =  0  .  (9.19) 

Differentiating  (9.19)  "totally"  with  respect  to  B  (i.e., 
keeping  only  a  and  p  fixed) 

and  inserting  5p-,/SB  =  (p-,  -  Po  )/B  yields  (9.I8)  again. 
The  constraint  (9-9)  can  be  brought  into  the  form 
(9.19)  if  a(x)  is  taken  to  be  a  function  of  (a,p)  alone; 
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the  particle  path  Invariant  Is  constant  on  each  magnetic 
line.   In  this  case  a  simple  computation  shows  the  equiva- 
lence of  (9.1^)  and  (9-l8).   But  this  special  case  Is  very 
uninteresting  since  It  Is  never  stable.   Prom 

-2  =  a(a,p)  |-  (9.21) 

Pi  ^2 

we  see  that  Pp/p-,  Is  unbounded  at  the  mirror  (p„  — >  O) 

which  violates  (9-l'^)-   In  other  words,  the  two  constraints 

(9-9)  and  (9-19)  are  Independent,  as  are  the  corresponding 

stability  criteria  (9-1^)  and  (9'l8),  except  In  an 

uninteresting  special  case. 

Just  as  In  MH,  we  can  obtain  an  explicit  elliptic 

system  by  specializing  to  axial  symmetry.   We  suppose  that 

Pp(''/',B)  Is  a  given  function.   We  have 


J  K  B  = ^  V^(A%) 

^v   r 

^      ^Pq  --*     ^Pq  _ 


^P2  —r       ^Pc 


The  term  A  i/   does  not  determine  the  nature  of  the  differ- 
ential system  because  there  are  second  derivatives  implicit 

In  Vj_-^  ^^  well  as  In  k.   A  more  appropriate  form  for  the 

2 
equation  Is  (9-7)  •   The  function  ^^{1/ ,B)   ==  Pg  +  B  /2m-q  Is 

known,  and  we  write  (9-7)  as 

k(b2/h^  +  P2  -  Pi)  -  ^  V,B  =  ^  V^  .      (9.23) 
The  right  side  of  the  equation  and  the  coefficients  of 
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Vj_B  and  of  k  Involve  only  first  derivatives  of  ii . 
Expanding  k  =  5b/Ss  we  find 

K  =  -i-^l  +  i'l)~^^^[i'li'^^    -   2^rV^z  +  ^z^rr^       (9-24) 

The  differential  operator 

2  2  2 

^2  S  _  2^  ^   ^^  +  ^2  a                 (5,25) 

r^2  rz  oroz     z  ^v  2              \^      ^i 

dz  or 

« 

is,   within   first   derivative   terms,    the    square    of 

■^     ^   -  if     ^  =  27rrB  4,-    •  (9-26) 

r   oz  z    or  as  ^  ' 


Insofar  as  characteristics  are  concerned,  the  dominant 

second  order  contribution  in  k  is  the  degenerate  hyperbolic 

2    2 
operator  -S  /Ss  .   A  similar  computation  shows  that  the 

2    2 
dominant  term  in  VlB  is  essentially  S  /Sn  where  n  is 

in  the  direction  VV'.   The  equation  (9.23)  is  therefore 

elliptic  as  a  linear  combination  with  positive  coefficients 

of  5  /Ss   and  S  /Sn  . 

Another  special  case  which  turns  out  to  be  purely 

elliptic  is  when  Pp(B)  is  independent  of  a  and  p  (Northrop 

and  Whiteman  [44]);  the  same  will  be  true  for  p^(B).   An 

elementary  computation  [using  (9-4)]  shows  that 

2 


where 


B  X  curl(aB)  =  P^  V^B  -  B  ok  (9-27) 


a(B)  =  (B^/H^  +  P2  -  Pi)/b2  .  (9.28) 
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We  can  therefore  write  (9- 7)  in  the  form 

B  X  curl(aB)  =  0  .  (9-29) 

Outside  the  plasma,  where  p-,  =  Pp  =  0,  we  have  a  =  ^/U-^- 

A  boundary  condition  J  =  0  therefore  Implies  that 

n 

curl(aB)  =  0  (9-30) 

which,  together  with  dlv  B  =  0  Is  a  nonlinear  elliptic  system 
for  the  vector  B. 

A  final  special  example  which  simplifies  In  an 
entirely  different  way  Is  the  marginal  case  with  p^(a,3) 
constant  on  each  magnetic  line .   One  of  the  criteria  for 
stability  Is  just  met  (the  significance  will  be  discussed 
In  Sec.  1^).   It  Is  convenient  to  Introduce  the  flux  coor- 
dinate I/'  on  which  p^(a,p)  =  const.   Thus  we  can  take  p^('i/^) 
as  given,  but  with  if   Itself  unidentified.   The  equation  Is 

vp*  =  pIv^  =   k(b^/m-q  +  P2  -  Pj)  •        (9.31) 

We  see  that  k  Is  orthogonal  to  a  i/'-surface ;  the  magnetic 
lines  are  geodesies  on  these  surfaces.   The  differential 
equation  Is  degenerate  hyperbolic  and  one  can  formulate  an 
Initial  value  problem,  given  the  vector  B  as  well  as  p^ 
on  an  Initial  transverse  surface,  say  midway  between  S' 
and  S"  In  Fig.  1  (the  domain  D  of  Fig.  1  Is  not  known 
beforehand  In  this  problem).   First  we  note  that  p-|(B,i/')  and 
Pp(B,T^)  are  both  explicit  once  p^{i')    Is  known  [using  (9-^)] 
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P2(B,^)  =  P^  -  B  /2H^ 
p^(B,^)  -  P^  +  B^/2n^  -  B(2pyn^)^/2 


(9.32) 


On  the  given  transverse  surface  S,  the  initial  values  of  the 
vector  B  and  of  p^>B  /2ia  determine  p-,  and  Pp  as  well  as  i^ . 
We  compute  VP^  a-s  being  perpendicular  to  B  from  its  projection 
on  S.   Thus  we  have  the  vector  k  on  S .   We  can  move  to  a 
neighboring  surface  S-,  with  a  corrected  direction  for  B  and 
employ  div  B  =  0  to  compute  the  increment  in  the  magnitude  of  B. 
Although  the  hyperbolic  operator  is  degenerate  and  involves 
only  derivatives  along  a  line,  the  term  Vp^  provides  coupling 
between  lines. 

A  very  degenerate  special  case  occurs  as  a  combination  of 
the  two  previous  examples,  viz.,  p^  =  constant  in  the  plasma 
domain.   The  equilibrium  equation  says  only  that  k  =  0;  all 
magnetic  lines  are  straight.'   In  both  this  case  and  the  previous 
case,  p^(a,p),  where  one  cannot  give  boundary  values,  the  plasma 
domain  can  only  be  joined  to  an  external  field  at  a  boundary 
with  a  surface  current.   A  more  appropriate  formulation  is  as 
a  very  unusual  free  boundary  problem  in  which  continuity  of 
the  vector  B  eliminates  the  surface  current.   For  example,  with 
axial  symmetry  one  would  look  for  a  solution  as  in  Fig.  2. 


Figure  2 
-  81  - 


A  solution  will  presumably  be  defined  uniquely  by  specifying 
P:^{i')    In  the  plasma  region  and  magnetic  field  boundary  condi- 
tions (vacuum  or  force  free  field)  at  the  exterior  boundary. 


10.   Equlllbrlum--GGP 

The  equations  to  be  satisfied  are  the  macroscopic 
momentum  balance  (9- 3),  the  two  Llouvllle  equations  for 
the  distribution  functions,  and  the  constraint  equation 
for  the  potential  j6.      Solution  of  this  problem  can  be 
separated  Into  an  analysis  of  the  distribution  functions 
and  potential  ^   followed  by  a  macroscopic  analysis;  the 
latter  has  already  been  done  In  the  previous  section. 

Consider  first  the  one-fluid  theory.   A  time -Independent 
solution  of  Llouvllle 's  equation  Is  f  =  f  [e ,[i,a ,^)   where 

£  -  ^  mv^  +  ^iB   .  (10.1) 

If  the  energy  contours  are  not  simple  In  the  (s,v) -plane, 
then  f  need  not  be  a  single -valued  function  of  e.   The 
macroscopic  connection  Is  given  by 

p^(B,a,p)  =  m^B /v^f°(|-  mv^ +M.B,ix,a,p)dvd^L        (10.2) 

Pp(B,a,p)  =  mB  /  ii.f  (^  mv  +  ii.B,iJ.,a  ,p  )dvdiJ.       (10.3) 
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THE  GUIDING  CENTER  PLASMA  NYO-l480-50 

"by  H.  Graa 

MF-^8 

Insert,  page  82  (after  Sec.  9) 

In  the  general  case  of  an  arbitrary  GCF  equilibrium, 

introduction  of  the  quantity  a(B,a,p.)  =  ^/^q   +  (Pg  "Pi)/^ 

as  in  (9.28)  yields  the  equation 

Sp,      Sp 
B  X  curl(aB)  =  "  ^  ^«  "  ^  "^p  (9-33) 

where  a,  p.,,  Pp  are  given  functions  of  (B,a,p).   Since  the 
highest  order  derivatives  enter  only  on  the  left  side  of 
this  equation,  plausible  iterations  can  be  devised  as  in  the 
scalar  pressure  case. 

For  a  generalized  axially  symmetric  problem  we  Introduce 
a  stream  function  if   and  Bg  as  in  (8.21).   Instead  of 
rBn  =  constant  on  i/^ -surfaces ,  we  now  find  crrBg  -   constant. 


Setting 


I  a^r^B^  =  g(-^)  (9.34) 


we  find  an  equation  for  i/ , 

A*^  +  v^.V  log  a  =  -  -^  g'  (•^)  -  ^  ^      (9.35) 


1   W.,^    r^  SPl 


The  left  side  is  second  order  and  elliptic  while  the  right 
side  is  first  order  through  the  dependence  of  a  and  p.,  on  B. 


At 


1    2 
The  limits  of  Integration  are  0  <  ii  <  0°  and  0  <  p  mv 

<  V-{B-.    -   B)  where  B-,  is  the  largest  value  of  B  taken  on 

the  given  line.   In  a  plasma  contained  within  mirrors, 

f  (e,i-i)  takes  values  in  a  sector  of  the  (£,|-l)  plane, 

e/l^  <  B-,  ,    and   f     must   vanish  before   reaching   e/[i  -  B-,  . 

a  given  point  on  the  line,  the  range  of  integration  is 

B  <  e/[x  <  B^. 

To  establish  that  p-,  and  Pp  are  compatible  with 

equilibrium  along  a  magnetic  line,  we  compute 

^  (^)  -  n?    /  v^H  ^   dvdM-.  (10.^) 

There  is  no  contribution  from  the  variable  limits  of 
integration  since  f  vanishes  there.   Noting  that 

7^=—^  (10.5) 

OE    mv  ov  ^ 

and  integrating  by  parts  (the  limits  again  vanish), 

^  (^)  -  -m  /  nf°dvdn  =  -  pp/B^  (10.6) 

which  is  the  required  relation,  (9.4).   The  remaining  analy- 
sis of  the  equilibrium  equations  follows  as  in  Sec.  9  with 
P-,(B,a,3)  and  p„(B,a,p)  given  compatibly. 

We  can  obtain  another  interesting  piece  of  information 
by  differentiating  Pp , 

^  (^)  =  m  /  n^  ^^  dvd^x  =  -Op  (10.7) 

B 
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But  this  moment  of  Sf^/Se  Is  exactly  what  appears  in  the 
stability  criterion  (7-2)  for  an  infinite  homogeneous 
plasma  (we  shall  also  confirm  it  as  a  necessary  condi- 
tion in  a  general  static  equilibrium  in  Sec.  13)-   Thus 
we  can  rewrite  the  stability  criterion  (7-2)  as 

B^/p.^  +  2p2  +  B^  ^  (4)  >  0  (lO-«) 

B 


or 

Bb~ 


>  0  .  (10.9) 


We  now  summarize  the  complex  situation  with  regard  to 
well-posed  problems.   The  ubiquitous  condition 
B  /m-  +  Pp  -  p.,  >  0  is  necessary  for  well-posed  equilibrium 
or  well-posed  motion,  for  both  the  GCP  and  the  GCP.   The 
microscopic  GCP  theory  has  a  second  stability  condition 
which  is  again  the  same  in  equilibrium  and  in  a  motion, 
but  it  can  be  written  in  two  different  but  equivalent 
forms,  (10.9)  and 

B^/^L^  +  2p2  -  B^G^  >  0  (10.10) 

respectively.   In  the  macroscopic  GCF  theory  there  are  two 
similar  criteria  for  equilibrium  and  motion,  viz.,  (10. 9) 


again  for  equilibrium  and 

'o  +  ^'^2   ^2^  ^^1 


B^/H^  +  2p„  -  p^/3p.  >  0  (10.-11) 


instead  of  (10. 10)  for  motion;  but  these  two  criteria  are 
independent  in  this  theory.   There  is  no  way  to  relate  the 
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2 
derivative  Bp^BB  to  the  quantity  Po/^P-,  (even  as  an 

inequality)  without  recourse  to  a  distribution  function.   On 

the  other  hand,  a  different  constraint  in  the  macroscopic 

theory,  viz.,  (9-9),  brings  the  static  and  dynamic  conditions 

together  In  the  common  form  (10.11). 

The  GCF  equilibrium  equations  are  essentially  incom- 
plete.  When  they  are  completed  in  one  way  they  yield  well- 
posedness  criteria  which  agree  with  the  GCF  equations  of 
motion.   When  completed  in  another  way  the  conditions  agree 
with  the  GCP  equations  of  motion  (cf.  Sec.  13). 

We  now  turn  to  the  two-fluid  theory.   Again,  solution 
of  Llouville's  equation  gives  f   =  ^.(e, )  where 

e+  =  I  "i+v^  +   \i-B  ±  jb    .  (10.12) 

We  must   choose  j6  to   Insure   charge   neutrality.      By   inspection 
of 

"+   =   m^B    /  f^(|  m^v^   +    iiB  +  ^)dvdM,  (10.13) 

we  conclude  that  n   is  a  monotone  decreasing  function  of  ^ 
(with  B  fixed)  if   f   is  monotone  in  e,  5f  /^e  <  0.   Both 
the  integrand  and  the  domain  of  integration  decrease  with 
an  increase  in  ^.      The  opposite  is  true  of  n_.   Since 
n,  — >  0  as  /i  — >  +  00  ,  there  is  a  unique  value  of  f)   at  which 
n  =  n_.   This  value  depends  on  the  parameter  B,   fi   =   />(b). 
Inserting  this  neutralizing  value  of  the  potential  makes 
n   =  n  a  function  of  B,  and  the  same  is  true  for  p,  and  Pp. 
We  may  now  continue  with  the  macroscopic  determination  of 
the  equilibrium  configuration  exactly  as  before.   This 
argument  depends  on  the  fact  that  the  determination  of  j25(b) 
is  entirely  local  (on  each  magnetic  line  and  at  each  location 
B)  and  can  be  done  first  before  analysing  the  global  pressure 
balance  as  in  Sec.  9- 
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The  situation  is  more  complicated  If  f   Is  not  monotone 
in  e.   For  example,  as  shown  in  Fig.  3j  there  could  be 


<t> 


Figure  5 


several  solutions  for  j6   which  yield  n  =  n_.   If  we  look 
at  the  original  physical  picture  from  which  the  quasi- 
neutral  formulation  was  extracted,  we  see  that  the  middle 
value  in  Fig.  3    is  unstable  (the  restoring  electric  field 
enhances  any  unbalanced  charge),  and  both  extremes  are 
stable.   In  an  infinite  homogeneous  plasma,  there  would 
be  no  difficulty  in  selecting  either  of  the  stable  choices. 
But  in  a  more  general  equilibrium  with  varying  B  along  a 
line,  we  can  find  a  situation  in  which  there  is  no  con- 
tinuous solution  jt){B)  .      There  are  two  choices  for  inter- 
mediate B,  and  a  unique  choice  for  large  and  small  B,  but 
there  is  no  choice  over  the  whole  range  which  is  physically 
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satisfactory.   Physically  one  would  expect  plasma  oscilla- 
tions of  finite  amplitude,  at  least  in  the  intermediate 
section.   Again  we  conclude  mathematically  that  the  problem 
is  not  well  posed  and  we  interpret  this  physically  as  an 
Instability.   We  note  that  not  every  non-monotone  profile 
f  (e)  produces  a  mon-monotone  r\{f)) ,    not  every  non-monotone 
n(/i)  produces  non-unique  values  f) ,    and  even  if  f)   is  multi- 
valued, limited  ranges  of  B  can  give  reasonable  (even 
though  non-unique)  solutions. 

There  is  a  great  deal  of  interesting  information  that 
can  be  obtained  by  differentiating  the  expressions  for  n,  p-,, 
and  Pp ,  just  as  in  the  one-fluid  case.   We  have 

n      =   m  B   /   f  °  dvdiJ. 

^  2„     ,      2„o 


p^   -  m^B   /  V   f^  dvdii  (10.14) 


and 


p^   =  m,B    /  l-Lf ,    dvdi-i. 


^    ^"B^    =  "   °1    +  ^o   Bb 

+ 

5        P2  t    _  t 

^    (-2)    =  -    ^2   +  C^ 


where 
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+ 


o, 


C^  =  -  m^  Jihfybe)^   [i'    dvdn 


+ 


+ 


(10.16) 


We  recall  that  h/hB   is  a  derivative  along  the  line  keeping 
(a,p)  fixed.   We  set 

n  =  n  +  n 


Pi  -  P^  + 

Pi 

^         + 

Pi  =  Pi 

P2  =  P2  + 

p" 

^           + 

P2  =  P2 

0^  =  01-, 

^; 

•^                   4- 

c    =  c 

r          r 

p- 


-  c. 


(10.17) 


Adding  and  subtracting  the  equations  (10.15)^  and  setting 


n  =  n_  =  p-  n  ,  we  obtain 


S  ,n. 


Sj? 


3b  ^B^  ~  "  ^1  ~  ^o  dB 
0  =  ^1  +  ^o  if 


S   Pi. 


P^ 


B 


^    (^)    =       £2  _  n  S^ 
^  ^  B  ''    ~  „2    B  ^ 


B 


P. 


3b  ^12^  -  "  ^2  ~  ^1  3b 


^  ^^2''  -  -  ^ 


2    ^1  Bb 


(10.18) 

(10.19) 
(10.20) 

(10.21) 

(10.22) 

(10.23) 


First  we  note  that  (10.20)  is  the  compatibility  relation 
between  p.,  and  Pp  for  the  macroscopic  equations.   Any  of  the 
other  equations  can  be  considered  as  nonlinear  integral 


equations  for  the  determination  of  /5(b)  (the  parameter  f) 
is  implicit  In  the  moments  n,  p-,  ,  Pp ,  C  ).   The  relation 
(10.19)  Is  equivalent  to  the  one  which  was  used  to  prove 
existence  of  /5(b).   Formula  (10.21)  is  the  equivalent 
(in  equilibrium)  of  the  general  constraint  equation  (5. 11) 
for  P   during  a  motion.   The  simpler  appearance  here 
(absence  of  factors  m  )  results  from  the  equilibrium 
relation  b-div  P  +  b-div  P~  =  0.   Eliminating  S/i/Sb  between 
(10.19)  and  (10.22)  gives  the  two  fluid  analogue  of  (10. 7), 


(10.24) 


The  inequality 


dB 

Pp 

-  -c* 

c* 

=   ^2    ■ 

^2 

E 

i^C^  > 

3Pn 

(10.25) 


,0  . 


for  monotone  f   is  stronger  than  the  similar  single  fluid 
result  with  Cg •   It  can  be  proved  by  integrating  the  explicitly 
negative  expression 

.  Pp     p    BG.  2  Sf° 
I  =   S^m_^  /  (^iB  -  J  —  m^v^  +  ^)  (^)  dvdl^  <  0  .  (10.26) 

The  constant  C^  plays  the  same  role  in  two-fluid  stability 
that  Cp  does  in  one  fluid. 
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11.   Variational  Formulatlon--Slngle  Fluid 

Our  starting  point  Is  the  observation  that  the  total 
energy  /(_  +  ^  Is  a  constant  of  the  motion  subject  to 
appropriate  boundary  conditions  (5. 20).   If  ^  =  ^^  Is  a 
minimum  of  the  potential  within  a  certain  specified  class 
of  neighboring  states,  we  conclude  that  there  is  stability 
in  the  form 

0  <  K  <  & 

(11.1) 

L  <  5  <  ?o  +  ^ 

for  a  perturbation  with  energy  constant  K  +  £  =  ^q  "^  ^' 
The  relation  between  existence  of  a  minimum  ^  and  a  score 
of  other  definitions  of  stability  is,  as  we  have  seen  in 
Sec.  2,  a  subtle  and  complex  affair.   We  therefore  accept 
the  minimum  property  as  the  definition  of  variational 
stability. 

The  utility  of  this  concept  depends  entirely  on  the 
precise  class  of  admissj.ble  functions  which  are  allowed 
Into  competition.   The  purpose  of  a  variational  approach 
is  to  replace  the  complex  study  of  those  states  which  can 
arise  during  a  motion  with  a  simpler  class  that  is  mathe- 
matically accessible. 

We  consider  a  tubular  domain  (Fig.  l)  within  which 
a  plasma  is  contained  by  mirrors,  i.e.,  there  is  a 


For  details  of  this  theory,  see  [I6] 
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relatively  large  value  of  B  at  each  end  of  the  tube.  There 
Is  not  necessarily  a  single  minimum  of  B  on  each  line;  this 
means  that  the  energy  contours  in  (s,v)  are  not  necessarily 
a  family  of  simple  closed  curves.  In  any  event,  even  with 
a  single  m.inimum  of  B,  the  electric  potential  f)  can  be  such 
as  to  complicate  the  topology.  And  in  a  torus,  the  energy 
contours  are  never  simple. 

We  introduce  the  admissible  magnetic  fields  and 
admissible  distribution  functions  independently.   For  B  we 
require  that  it  cover  the  domain  simply  and  satisfy 

div  B  =  0  ,  B  =  Va  x  V  P  (11-2) 

and  a  boundary  condition;  either 


BCI:   B  =  0  on  S  ,  a  and  p  given  compatibly  on  S' 

n  O  r-   o 

and  S" 
BGII:   B  ==  0  on  S  ,  compatible  B  >  0  on  S '  and 
B  <  0  on  S"  . 

The  compatibility  is  taken  care  of  (as  is  conventional  in 

all  variational  formulations)  by  requiring  that  there  exist 

at  least  one  admissible  vector  field  B.   In  BGII  this  implies 

only  that   /  B  +  /  B  =0,  while  in  BCI  there  is  a 
g,  n   g„   n 

restriction  on  the  topology  of  the  curves  a  =  const., 

P  =  const,  as  well  as  on  the  total  flux,  /  dadp  =  /  dadp . 

S'        S" 
BCI  corresponds  to  conducting  end  walls  with  tied 

lines  and  BCII  is  appropriate  for  insulating  walls  or  a 
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vacuum  surrounding  the  plasma.   For  a  contained  plasma, 
the  pressure  Is  required  to  approach  zero  before  reaching 
any  physical  walls.   The  distinction  between  BCI  and  BCII 
Is  between  a  pressureless  plasma  (force  free  field)  or  a 
vacuum  acting  as  a  buffer  zone  around  the  primary  plasma. 
Whether  a  given  physical  plasma  is  a  vacuum  or  a  force 
free  field  depends  on  the  ratio  of  two  small  quantities, 

the  plasma  energy  compared  to  magnetic  energy, 

1  2 

P   =  (p"  Pi  +Pp)/(B  /2m.  ),  and  the  gyro  radius  compared  to 

a  plasma  dimension  t^  =  'K/r.      If  pV^  »  1,  then  sufficient 

current  can  be  carried  to  appreciably  alter  B  from  a 

vacuum  field  (curl  B  =  O)  to  a  force  free  field 

(curl  B  X  B  -  0).   Since  r\   is  the  GCP  expansion  parameter, 

it  is  implicit  that  PV^l  »  1  wherever  the  GCP  equations 

are  used.   A  vacuum  domain  must  be  treated  separately, 

using  a  form  of  Maxwell's  equations.   On  the  other  hand,  a 

vacuum  can  arise  as  a  special  case,  curl  B  =  0,  even  when 

the  fluid  itself  is  considered  to  be  a  pressureless  plasma, 

e.g.,  as  a  result  of  the  boundary  condition,  J  =  0. 

The  distribution  function  f  is  fixed  in  its  dependence 

on  (a,f3,[i.)  and  is  varied  only  in  Its  dependence  on  the 

dynamical  coordinates  (a,p).   This  is  consistent  with  the 

equations  of  motion  which  preserve  (a,p,ii.).   Any  motion 

subject  to  Llouvllle's  equation  is  incompressible,  and 
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preserves  area  In  (a,p)  or  (s,v).    On  a  given  line  (a,p) 

and  for  a  given  value  I-l,  we  shall  consider  as  admissible 

all  functions  f(a,p)  which  are  compatible,  through  an 

incompressible  mapping,  with  a  given  reference  function. 

To  be  precise,  take  a  fixed  monotone  function  f*(A)  and 

consider  as  admissible  any  function  f(a,p)  which  is  equi- 

measurable  with  f*:  if  A   is  the  area  of  the  domain  on  which 

o 

f(o,p)  >  f  ,  then  f  =  f*(A  ).   The  domain  f  >  f  need  not 
\     '^  I    -^       q'  o       ^o'  o 

be  connected,  and  two  admissible  functions  need  not  be 
topologically  similar.  More  completely,  the  reference 
function  f*  is  a  function  of  (A,|i-,a,p). 

The  restrictions  on  B  and  f  are  not  entirely  indepen- 
dent; the  field  must  be  such  that  f  does  not  "spill"  over 
the  mirrors . 

For  some  purposes,  e.g.,  to  obtain  first  and  second 
variations,  it  is  more  suitable  to  consider  deformations 
of  B  and  f  rather  than  mappings.   This  gives  a  more 
restricted  admissibility  class,  e.g.,  the  topology  of 
f -contours  is  fixed.   Although  any  increased  flexibility 
introduced  by  mapping  as  compared  to  deformation  is  irrele- 
vant physically  and  incompatible  with  solutions  of 


For  purposes  of  defining  the  admissible  class  the 
distinction  between  s  and  a  (3-22)  can  be  ignored,  ^  =  1. 
The  parameter  (^  is  essential,  however,  in  performing  a  con- 
tinuous variation  of  an  admissible  function  depending  on  a 
parameter. 
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Llouvllle's  equation,  we  shall  on  occasion  find  this  flexi- 
bility to  be  -useful,  and  we  shall  even  generalize  the 
admissible  deformations  to  correspond.   The  reason  is  that 
certain  auxiliary  minima  which  are  not  attained  by  any  member 
of  the  lesser  admissible  class  are  attained  (without  altera- 
tion of  the  minimum  value)  in  the  greater  admissible  class. 

The  variation  of  B  and  f  is  most  conveniently  described 
in  terms  of  mappings  of  the  underlying  space,  x  — >x'  carrying 
the  magnetic  field  lines  and  (a.p)  — >  (o',p')  carrying  f.   We 
have  seen  in  Sec.  3  that  any  deformation  of  a  solenoidal 
field,  div  B  =  0,  can  be  represented  by  a  "velocity"  U  such 
that 

^  +  curl(B  X  U)  =  0  ,  U-B  =  0  .  (11-3) 

This  is  a  variational  equation  and  t  is  merely  the  parameter 
which  describes  the  family  B(x,t);  SB/St  at  t  =  0  would  be 
more  conventionally  represented  as  5B  and  U  is  a  displacement 
Bx.   We  similarly  represent  a  deformation  of  f  in  terms  of 
the  incompressible  velocity  field  w  in  (a,p)  =  z 


Bt  +  ^  •  ^  =  ° 

(11.4) 

-^  •  w  =  0 
dz 

A  variation  of  (B,f)  is  described  by  a  pair  (U,w). 
U  is  arbitrary  subject  to  U-B  =  0  and  a  boundary  condition 
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BCI:   U   =  0  on  S^ 

U  =  0  on  S' ,  S" 

BCII:   Uj^  =  0  on  S^ 

^  U  X  B-dx  =  0  on  S' ,  S" 
C 
for  every  closed  curve  C  on  S '  and  S" .   Note  that  U^  may 

be  different  from  zero  on  S'  using  BCII;  U  does  not  repre- 
sent motion  of  plasma,  only  of  field  lines.   The  end  of 
the  plasma  Is  determined  by  f  =  0;  this  point  Is  not 
carried  perpendicular  to  B  with  U.   The  velocity  w  Is 
arbitrary  subject  to(S/Bz)-  w  =  0. 

For  the  one  fluid  theory  we  take 

J[B,f  ]  -  /  2^  B^dx  +  /  e  f  dQ  (11-5) 

as  the  variational  function.   Performing  a  first  variation 
with  respect  to  f  alone,  holding  B  fixed,  yields  the 
variational  condition  that  f  be  a  function  of  e, 

f  =   f°(£,n,a,p)  .  (11.6) 

This  Is  a  local  condition,  and  f  need  not  take  the  same 
value  on  disjoint  e-contours.   But  a  very  simple  analysis 
shows  that  f  Is  an  absolute  minimum  If  and  only  If  f   Is 
single-valued  and  a  monotone  decreasing  function  of  e, 

i^  <0  (11.7) 

Examination  of  the  second  variation  of  ^  (also  varied  with 
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respect  to  f  alone)  reveals  a  mlnlm-um,  5  ^  >  0,  even  If 
f°(£)  Is  not  single-valued,  provided  that  each  branch 
Is  locally  monotone.   Such  a  multivalued  equilibrium  f 
provides  a  minimum  in  the  linear  theory  but  is  not  a 
minimum  with  respect  to  finite  perturbations,  no  matter 
how  small. 

We  can  now  state  that  a  necessary  condition  that  any 
state  [B,f]  be  stable  is  that  f   be  monotone  in  e;  single- 
valued  and  monotone  in  a  finite  theory,  and  piecewise 
monotone  in  an  infinitesimal  theory. 

Performing  a  first  variation  with  respect  to  B  yields 
the  result 


.(div  P  -  rr^  curl  B  X.  BJ  -  0  (H-S) 


for  an  interior  variation  U.   The  boundary  variation 
vanishes  identically  under  BCI,  tied  ends,  and  yields 
the  natural  boundary  condition 

Jn  -  0  (11-9) 

under  BCII,  free  ends.   Thus  we  obtain  variationally 
exactly  what  was  termed  an  equilibrium  state  in  Sees.  9 
and  10,  together  with  a  natural  boundary  condition  when 
it  is  appropriate.   The  natural  boundary  condition  J  =  0 
implies  that  the  buffer  zone  surrounding  the  plasma 
satisfies  curl  B  =  0  even  if  it  is  a  force-free  field  and 
is  capable  of  carrying  current. 

The  previous  analysis  of  the  variation  of  ^  with  respect 
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to  f  alone  holds  not  only  for  an  equilibrium  magnetic 
field  but  for  any  given  admissible  B.   We  can  therefore 
minimize  ^  in  two  steps.   For  any  fixed  admissible  B  we 
have  a  unique  function  f  =  f  ,  depending  on  B,  which 


provides  an  absolute  minimum.   If  we  associate  this  f 


IT 


W 


ith  B,  then  $  becomes  a  functional  of  B  alone.   We  call 


this  special  class  of  admissible  pairs  [B,f  ]  a  pessimistic 
variation.   Even  though  this  class  of  variations  is 
special,  it  is  evident  that  all  minima  of  ^  can  be  found 
within  this  class. 

The  existence  of  a  pessimistic  variation  depends 
vitally  on  allowing  mappings  rather  than  deformations  in 
the  (s,v)-plane.   In  a  complex  topology  with  islands  (Fig.  k), 
a  variation  of  B  may  require  the  transfer  of  f-values  from 
the  center  of  one  island  to  the  center  of  the  other  without 


Figure  4 
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disturbing  the  intermediate  values.   Of  course,  such  a  map- 
ping could  be  approximated  by  a  flow.  Pig.  5-   Restriction 


Figure  5 


of  the  admissible  class  to  the  more  physically  relevant 
class  of  deformations  would  yield  a  functional  ^  which  has 
only  a  lower  bound  but  no  minimum.   The  wider  class  of 
mappings  Is  what  allows  the  simple  characterization  of  the 
minimum  with  respect  to  f  and  the  Introduction  of  the 
pessimistic  variation  In  a  complex  geometry.   It  is  evident 
that  these  mappings  cannot  be  represented  by  an  incompres- 
sible flow  w.   A  suitably  generalized  class  of  flows  can  be 
introduced,  but  we  shall  not  require  it. 

It  is  not  immediately  evident  that  one  can  Introduce 
the  concept  of  a  pessimistic  variation  in  an  infinitesimal 
theory  (first  and  second  variations)  where  one  is  necessarily 


restricted  to  deformations,  not  mappings.   In  any  finite 
pessimistic  variation  it  is  only  the  sum  of  the  areas  of 
distinct  f-contours  which  is  preserved  rather  than  each 
individual  area  as  would  follow  from  an  incompressible  w. 
But,  if  we  examine  the  second  variation  6  ^  we  are  able  to 
construct  an  explicit  and  unique  minimizing  w  in  terms  of 
the  given  variation  U  simply  by  completing  squares.   This 
pessimistic  w  =  w   turns  out  to  be  incompressible,  except 

^  IT 

for  a  dipole  type  singularity  at  the  critical  point  where 
two  islands  join.   The  reason  for  this  is  geometrically 
clear.   For  an  Infinitesimal  variation,  the  interior  of 
each  island  is  isolated  from  the  others  and  can  be  assigned 
an  incompressible  w.   This  preserves  the  area  of  each 
individual  interior  f-contour.   Since  a  variation  of  B  can 
transfer  area  between  islands,  the  singularity  in  w  is 
required.   Although  the  area  of  every  interior  f-contour  is 
preserved,  the  areas  of  the  critical  contours  are  not.   We 
remark  that,  although  the  transferred  area  is  of  first 
order,  the  change  in  energy  induced  by  the  transfer  of  f 
values  between  islands  is  of  order  higher  than  second 
(i.e.,  negligible  In  the  second  variation);  otherwise  no 
minimizing  w   could  be  found  for  a  piecewlse  monotone 

T 

f  (e).   Even  if  Sf  /5e  is  discontinuous  between  two  branches, 
f   Is  stationary  as  a  function  of  (s,v)  at  a  critical  point. 
There  is  another  problem  where  the  Hamlltonian  require- 
ment of  constant  area  of  each  f-contour  must  be  dropped 
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to  obtain  a  pessimistic  variation.   In  the  infinite  homo- 
geneous plasma  the  energy  contours  are  horizontal  straight 
lines  in  equilibrium.   For  finite  energy  perturbations  we 
consider  variations  of  B  which  are  bounded  in  extent  (compact 
support).   The  rule  for  a  pessimistic  variation,  placing  a 
value  of  f  on  the  energy  contour  of  appropriate  area,  cannot 
be  applied  because  the  areas  are  infinite.   Formal  minimization 
of  the  second  variation  with  respect  to  w  gives  an  explicit 
result,  but  one  which  violates  the  area  criterion.   To 
interpret  this  result,  we  choose  a  large  region  surrounding 
the  perturbed  domain,  minimize  within  this  region,  and  then 
let  the  region  grow.   Within  the  large  region  the  minimum 
is  given  by  placing  the  largest  f-values  on  the  smallest 
e-contours,  preserving  area.   The  minimizing  f  will  be 
discontinuous  on  the  boundary  of  the  large  region  (Fig.  6). 


Figure  6 


As  this  region  grows,  the  value  of  f  which  Is  placed  on  a 
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given  e-contour  will  approach  the  original  unperturbed 

value,  f  (e).   But  the  limiting  result  f  =  f  (e)  (i.e., 

the  original  values  of  f  applied  to  the  perturbed  e-contours) 

violates  the  area  constraint.   The  area  discrepancy  which 

Is  caused  by  the  perturbation  of  B  on  a  given  contour  Is 

finite,  but  It  is  distributed  over  a  longer  and  longer 

length.   Area  is  not  preserved  in  a  nonuniform  limit,  but 

the  energy.  Jet,    can  be  seen  to  converge  uniformly; 

(since  /  ef  is  stationary  for  f  =  f  ,  the  energy  discrepancy 

is  of  second  order  in  the  length  of  the  region). 

An  Infinitesimal  analysis  gives  the  same  result;  the 
pessimistic  f  is  the  unperturbed  function,  f  ,  of  the 
perturbed  argument  e.   There  is  no  proper  minimum  within 
the  original  admissibility  class,  but  the  value  of  the 
minimum  is  unchanged  when  we  extend  this  class  to  Include 
the  pessimistic  variation. 

Although  minimization  with  respect  to  f  is  essentially 
explicit,  the  subsequent  minimization  with  respect  to  B  is 
not  and  requires  intensive  study.   We  list  several  forms 
for  the  second  variation,  6  ^,  both  for  a  general  variation 
and  a  pessimistic  variation.   The  pessimistic  variation  w 

IT 

is  given  in  terms  of  U  through  the  formula 

where 

g  =    -mv2(b  •  ||)  +  nB(b  .  Il  -  dlv  U)  (11.11) 
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and 

J  -^[e   -^iB(s)]"^/^  ds 
<i'>=   -^72 (11.12) 

is  a  mlcrocanonlcal  average.   Although  the  formula  (ll.lO) 
gives  only  the  component  of  w  normal  to  the  energy 
contour,  Incompresslbillty  serves  to  complete  the  definition; 
moreover,  the  normal  component  is  all  that  appears  In  any 
of  the  later  formulas.   In  the  infinite  uniform  plasma,  the 
pessimistic  variation  merely  drops  the  term<g>^  in  (11. 10). 
But  in  a  uniform  plasma  with  finite  (e.g.,  periodic) 
boundary  conditions,  the  full  expression  (11. 10)  would  be 
kept . 

The  formula  (11.10)  is  easily  interpreted  once  we 
identify  g  as  the  pessimistic  variation  of  e  (at  fixed  area) 
Induced  by  the  displacement  U.   A  pessimistic  variation 
requires  that  an  e-contour  map  into  another  e-contour, 

||  +  „.||=*{.)   .  (11.15) 

But   a   simple   computation   shows   that  <w  •  Se/Sz>=   0;    therefore 
■i/^(e)    =<5e/St>and    (11.10)    follows. 

For  the    second  variation  we   find 

&^^  =  T(U)    -   Q(w    )   +   Q(w    -   w    )  (11.1^) 
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where 


Q(w)  =  -{(w  -If)   ^dQ  (11.15) 

and  we  shall  define  T(U)  below.   Q  is  positive  if  hf  /he   <  0, 
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and  w  =  w   clearly  minimizes  6  ^  at  fixed  U.   If  f   i£ 


not  monotone,  we  can  take  U  =  0  and  easily  find  a  varia- 
tion  w  which  makes  6  ^  negative.   The  pessimistic  second 


o, 


V 


ariation  [defined  only  for  plecewise  monotone  f  (e)]  is 

6^$^  =  T(U)  -  Q(w^) 


,2  Sf° 


=  T(U)  +  /(g  -  <zy)      ^   d^      (11.16) 
=  T(U)  +  /(g^-<g>^)  11^  df]  . 


An  expression  for  6  ^  directly  in  terms  of  5f  rather 
than  w  is 

5^J  =  T(U)  +  2  /g(5f)dn  -  /(6f)^(Sf°/5e)"-'-  df^     (H-l?) 

The  transverse  variation  T(U)  is  a  special  variation 
taking  w  =  0;  in  other  words,  f  as  a  function  of  the  canonical 
coordinates  (a,p)  is  frozen  on  a  given  line.   It  happens  to 
be  identical  to  the  macroscopic  GCP  transverse  variation  (but 
in  the  GCP  a  general  allowed  displacement  U  would  also  have  a 
component  parallel  to  B).   We  have 

T(U)  =  /[P3_G^  +  P2G2  +  (U-VU)-div  P  +  G^]dx     (11. I8) 
=  /[PlG3_  +  PgGp  +  ^m^^^  (11-19) 

"  2U~  f  (B-Vu^)B-dS  -  2^  f  B^(U.'7u).dS 
^o  ^  o  / 

where 
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(11.20) 


G  =  —  ([curl(U  X  B)]^  +  U-curl(U  X  B)  x;  curl  B  j 

Boundary  terms  have  been  dropped  whenever  they  vanish  for 
both  BCI  and  BCII.   The  first  boundary  term  In  (11.19) 
vanishes  on  the  lateral  boundary,  S  ;  both  terms  vanish 
at  the  ends  S'  and  S"  for  BCI.   But  at  least  one  boundary 
term  is  always  present,  except  for  a  local  disturbance  in 
which  U  becomes  zero  before  reaching  the  boundary. 

The  second  form  (11.19)  is  useful  for  making  certain 
estimates  since  it  contains  no  derivatives  of  the  equilib- 
rium parameters  p-,  ,  Pp ,  and  B.   Interchanges  are  more 
easily  studied  with  the  first  form  (11.18).   It  is  also 
fortunate  that  (11.18)  is  insensitive  to  whether  the  ends 
are  tied  or  free;  this  is  what  allows  the  Kruskal  and 
Oberman  formula  (originally  derived  for  tied  ends)  to  be 
used  by  Taylor  and  others  for  interchange  stability. 

Another  useful  form  is  obtained  for  the  pessimistic 

2 
variation  by  expanding  the  term  in  g  , 
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(11.21) 


2 

+  /  (B^/2H^  +    2p2    -   B^C2)(b  •  ||   -   div  U)^dx 

2    Sf° 
+  /  Vp^-  (U  dlv  U    -  U-  VU)dx  +  /<g>      3^  dn 

_      1      (    (B-VU^)B-dS    -   p^  I  B^(dlv  U)U-dS 

In  this  form  the  boundary  terms  appear  only  at  the  ends  S' 

and  S"  and  only  for  BCII. 

We  have  introduced  a  variational  formulation  as 

originally  suggested  by  the  equations  of  motion.   But  this 

procedure  can  be  reversed.   The  pessimistic  variation  6  ^^ 

is  a  quadratic  functional  of  the  magnetic  displacement  U. 

Prom  this  we  can  obtain  a  symmetric  system  of  pessimistic 

equations  of  motion,  taking  5  ^   as  the  potential  and 

/  i  p  U  dx  as  the  kinetic  energy  [the  notation  is  more 

conventional  if  we  denote  the  displacement  by  ^  and  take 

U  =  C  as  the  velocity]. 

The  pessimistic  variation  of  f. 


5f   -  -w_.^-  -w^-lll?^-  [g  -<g>]  1^     (11.22) 


is  clearly  the  same  as  the  adlabatic  variation  (5-33)  (g  Is 
the  perturbed  Hamiltonlan) .   Thus  the  pessimistic  equations 
of  motion  agree  in  form  with  the  adlabatic  equations  of 
motion  which  were  introduced  in  Sec.  5j  even  though  the 
motivation  and  the  range  of  application  differ  widely. 
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Prom  their  derivation,  the  pessimistic  equations  of 
motion  are  limited  to  monotone  f  (e).   Variational  bounded- 
ness  as  we  have  defined  it  is  necessary  and  sufficient  for 
exponential  boundedness  of  these  equations;  any  nonexponen- 
tlally  growing  solutions  are  extraneous.   There  is  no  claim 
that  pessimistic  motions  approximate  actual  motions,  only 
that  the  stability  criteria  are  related. 

With  the  adlabatlc  interpretation  of  these  equations 
(e.g.,  as  used  by  Andreolettl  [23]),  we  are  not  limited  to 
monotone  f  (e).   On  the  other  hand,  they  give  information 
about  stability  only  if  an  unstable  adlabatlc  motion  is 
verified  to  be  slow;  we  have  seen  that  this  can  be  easily 
violated  even  when  f  (e)  is  monotone. 

Choosing  the  most  useful  Interpretation  in  each  case, 
exponential  stability  of  the  adlabatlc  equations  is  equivalent 
to  variational  stability  when  f  (e)  is  monotone;  and  exponen- 
tial stability  of  the  adlabatlc  equations  is  necessary  for 
exponential  stability  of  the  actual  equations  of  motion 
for  general  f  (e),  but  only  if  it  is  further  verified  that 
the  marginal  motion  is  slow. 
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12 .   Two-Fluid  Varlatlon--Existence  of  the  Potential 

For  the  two-fluid  version  we  specify  two  reference 
functions  f*(A,^L,a,^)  subject  to 

/  f*  dAdii  =  /  f^  dAdM.  .  (12.1) 

In  other  words,  there  are  equal  numbers  of  electrons  and 
ions  on  each  line.   Admissible  functions  f   are  chosen  to 
be  compatible  with  f  and  also  subject  to  charge  neutrality 
at  each  point  on  the  line 

/  f  d^idp  =  /  f_d^Ldp   .  (12.2) 

Note  that  the  electrostatic  potential  is  not  mentioned  in 
this  formulation  of  a  charge  neutral  admissibility  class. 
The  charge  neutrality  requirement  (12.2)  means  that 
the  two  variations  (w  ,w  )  are  not  independent;  they  are 
subject  to  the  constraint 


1 


Sf^      Sf 
(w+  -^  -  w_  •  ^)  dpdix  =  0  (12.3) 


or,  setting  w  =  (w  ,w  ), 


/  w'^f^dpdM-  =  /  w  f  dpdM-  .  (12.^) 


Writing 


—    12 
e,=2   i^4.v  +  M-B 


e^  =  ^  m^v^  +   [iB  ±  jt) 


(12.5) 


we  have 
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$  =  /  2^  B^dx  +  /  (e^f^  +  e_f_)dfi  (12.6) 

o 

=  I  2k-   B^^^  +  /  (^+^+  +  e.fjdfi  (12.7) 

o 

The  terms  Involving  ^   Implicit  in  £   cancel  out  of  (12. 7) 
because  of  charge  neutrality,  (12.2),  for  any  choice  of  the 
function  j6{o)  .      It  is  an  easy  matter  to  verify  that  ^  is 
stationary  if  f  and  f_  are  functions  of  e  and  e_  respec- 
tively.  The  converse,  that  all  charge  neutral  stationary 
states  take  this  form,  is  more  difficult;  this  amounts  to 
justifying  the  Lagrange  multiplier  rule  (which  yields  this 
result  instantly).   To  verify  this,  it  is  sufficient  to  show 
that  given  an  arbitrary  w  ,  we  can  find  a  related  w_  satis- 
fying (12.3)  which  makes  the  second  term  vanish  in 

Sf^         Sf 
5(5  -yr/)    =   -  f  (£_^w_^  •  ^  +  e_w_.  ^)dn  .       (12.8) 

Alternatively,  it  suffices  to  find  a  class  of  functions 
w  (a,p)  which  is  general  enough  to  allow  us  to  assign  an 
arbitrarily  given  value  to  the  second  term  in  (12.8)  while 
satisfying  /  w  f  dpd^.  =  0.   The  trial  function  w  =  [a'(p)b(a), 
-a(p)b'(a)]  is  easily  seen  to  satisfy  these  requirements. 

We  note  that  j6{o)    is  arbitrary  in  this  entire  analysis 
of  the  first  variation.   The  function  ^{o)    is  determined  by 
the  condition  that  the  equilibrium  be  charge  neutral:   to 
any  given  pair  f  (A),  a  potential  jt)   can  be  found  which  will 
make  f, (e.)  charge  neutral.   We  recall  that  we  were  able  to 
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prove  a  similar  theorem  In  Sec.  10.   There  the  functions 
f  (e, )  and  f_(£  )  were  assumed  to  be  given,  and  the 
neutralizing  ^(s)   was  found  locally,  at  each  point  s.   But 
If  f  are  fixed  as  functions  of  A  Instead  of  e,  variation 
of  ^{s)    changes  n  (s)  along  the  entire  line,  and  j6{s)    must 
be  found  from  a  nonlinear  Integral  equation  which  we  shall 
present  later. 

Just  as  It  Is  necessary  to  prove  the  existence  of  a 
neutralizing  potential  In  equilibrium.  It  Is  necessary  to 
find  the  varied  potential  in  a  pessimistic  variation.   Note 
that  the  potential  j6   has  no  significance  In  a  generel 
variation,  but  It  must  be  defined  In  a  pessimistic  varia- 
tion where  f  Is  constrained  by  the  auxiliary  minimization 
to  be  a  function  of  e.   The  relation  (5-11)  which  determines 
j6   In  a.   motion  has  no  relevance  to  the  variational  analysis 
which  does  not  employ  the  dynamical  equations  of  conservation 
of  momentum.   The  varied  potential,  i^   =  bjt) ,    Is  found  to 
satisfy  the  Integral  equation  [l6] 

L[^]  =  X  (12.9) 

where 

L  =  L_^  +  L_ 


Lji/]    =  m^  !  {if    -  <-^X  )(Sf7Se),    dvd^i  (12.10) 

'T  'T  T"  -T 
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and 


X  -   L_[gJ  -  L^[g^]  (12.11) 


and  the  g  are  given  functions  (ll.ll)  in  terms  of  U.   The 

mean  values  <  >   are  taken  on  distinct  contours  e  =  const. 

The  term  ii   can  be  taken  outside  the  Integral, 

+ 
L^[-^]  =  -Cq  1^  +  K_^[l^]  (12.12) 

+ 
where  C~  is  the  previously  defined  moment  (10.16)  and  K 

is  a  symmetric  integral  operator  with  the  kernel 
k^(s,s')  = 
J[£-HB(s)+j2^(s)]"^/^[£-HB(s>);/5(s')]-^/2  T^(e,n)ded^i 


The  quantity  bA/'de   is  essentially  the  orbital  period  of  a 
particle  with  given  e. 


M|)'^'|[s-.B,;.]-V2ds.|[i^     (12. U) 


The  integral  equation  is  to  be  solved  for  j6{s)    on  a  given 
line;  a  and  p  are  parameters. 

We  consider  only  monotone  f  (e),  since  this  is  a 
necessary  condition  for  stability.    This  implies  that  the 
kernel  k(s,s')  is  a  positive  function.   It  also  follows 
easily  that  both  -L  and  K  are  positive  operators. 
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(^,L^)   =  Ji/L[^]ds   =  m  J  {i/    -  <if>)      ^  dvd^ids  <   0 

2  (12.15) 


L  Is  not  definite;  Li^  =  0  if  and  only  if  V'  =  <i>,    i.e.,  i/ 
is  constant.   We  verify  that  the  inhomogeneous  term  X  1^ 
(12.11)  is  always  orthogonal  to  a  constant.   The  solution  i/ 
will  therefore  be  unique  only  up  to  an  added  constant. 

The  spectrum  (and  in  particular  any  discrete  eigen- 
value) of  K  relative  to  C   does  not  exceed  unity, 

(12.16) 
0  <  A  <  1 

This  is  an  immediate  consequence  of  the  positivity  of  K 
and  -L, 

{i',Ki')    <  {i',G^1')       .  (12.17) 

It  is  easily  verified  (this  amounts  to  A  =  1  being  an  eigen- 
value) that 

/  k(s,s' )ds'  =  C  (s)  .  (12.18) 

Therefore   the    operator 

M[^]    =  ^  Kiif]  (12.19) 

o 

is  a  mean  value  with  positive  weight,  and  the  norm  of  M  is  not 

greater  than  unity  taking  a  maximum  value  norm  on  i^ .      Prom 
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this  fact  alone  It  follows  that  Iterations  of 

^  =  ^  K[^^]  -  ^  (12.20) 

o        o 

wll,l  converge  in  absolute  value  if  "X/C   is  bounded.   To 
estimate  the  speed  of  convergence  we  must  show  that  the 
eigenvalue  A  =  1  in  (12. 16)  is  isolated.   In  terms  of  the 
new  variables 


.1/2 


o 


C-V2  X   =  X 


12.21) 


we  have  the  integral  equation 

?  =  K[l!]    -X  (12.22) 

[with  the  same  spectrum  as  ( 12.16)]  where 

k(s,s')  =  [C  (s)C^(s')]-l/2  k(s,s')        (12.23) 


Examination  of  the  kernel  (for  example,  taking  t  to  be 
bounded)  shows  that  it  is  bounded  except  near  s  =  s'  and 
near  s  =  s-,  or  s '  =  s,  where  s,  is  the  limit  of  the  plasma, 
f  (s-|)  =   0.   Near  s  =  s'  there  is  a  logarithmic  singularity 
k  ^  log|s  -  s' I ;  k  is  bounded  at  the  boundary  of  the  domain 
except  at  the  corners,  s  =  s'  =  s-,  where  k  grows  at  most 
as  (ss')  ^     .      Thus  k  is  square  integrable,  and  the  Neumann 
series  of  (12.22)  converges  uniformly  in  Lp •   Given  the 
existence  of  an  Lp  solution  J   of  (12,22),  it  follows  that 
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K[if]    Is  continuous  and  bounded,  and  V'  Is  also  continuous 
provided  that  this  Is  true  of  the  Inhomogeneous  term. 

The  existence  of  a  ;25  to  neutralize  f_j_(A,U.)  in  equilib- 
rium depends  on  the  solution  of  a  nonlinear  Integral  equa- 
tion which  is  closely  related  to  the  linear  equation  just 
studied •   We  have 

n  /b  =  m  /  f  °  dvdii. 

-  (I  ni^)^/2  /f°[£  -  HB  +  jdr^^^dedix 

5fO         (12.24) 

=    -(2m^)^/2  /[e^-ixB  +  ^]^/^  {^)dkdv. 

Here  n  depends  on  f)   not  only  explicitly  but  in  e  when 
expressed  as  a  function  of  A  and  \i.      The  charge  neutrality 
condition  is  G[jb]    =  0  where 

gW    =    (m^)^/2  /[e^  -  HB  -  ^]^/^   ^   dAdu 

^f*         (12.25) 

-  (m_)^/^  /[e_  -  ixB  +  J25]^/^  ^  dAd^L 

If  we  set  /i  =  0  in  this  expression  we  obtain  a  definite 
function  G[0]  =  N[ s ] ,  essentially  the  residual  difference 
n   -  n_  when  the  potential  is  ignored.   We  consider  ^{s,t) 
to  be  the  solution  of  the  integral  equation 
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G[j6]    =  (1  -t)N(sj  .  (12.25) 

Att  =  0,j?5=0  1sa  known  solution;  we  wish  to  extend  the 
solution  to  t  =  1.   Differentiating  with  respect  to  t 
(i.'e.,  performing  a  variation)  we  obtain 


L[^]  =  2^/^N(s)  .  (12.27) 


This  is  a  linear  integral  equation  (the  one  we  have  already- 
studied)  for  h^/bt   with  a  kernel  which  depends  on  j25.   In 
other  words,  it  is  an  ordinary  differential  equation  for 
j6{t)    in  a  function  space.   We  shall  not  carry  the  analysis 
further,  but  it  seems  quite  likely  that  existence  can  be 
proved  in  this  way. 

The  analysis  of  the  second  variation  is  similar  to 
the  one-fluid  case,  except  that  a  pessimistic  variation 
Involves  the  variation  5/5  which  must  be  obtained  from  the 
integral  equation  (12.9).   The  pessimistic  variation  is 
given  by 

w     •  ^r —  =  <6e  ,>  ,    -   5e  . 
IT      oz  +   +  + 

(12.28) 
-  <g+>+    -   g+  t    (<6^>+    -    5^) 

where  g  is  given  by  (11. 11)  as  before.   For  a  general 
variation  (b/)   does  not  appear),  formulas  (11. l4),  (11.15)5 
(11.17)  are  merely  summed  over  the  two  fluids,  and  (11. I8), 
(11.19),  (11.20)  are  identical.   For  a  pessimistic  variation 
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we  have 


5^$^  =:  T(U)  -  Q^(U) 


(12.29) 


where  T(U)  is  identical  to  the  one-fluid  expression  and 

,      =    -E       I  Iff       -  <R    •>. 

IT 


+ 


-\  j(g+    -  <S+>+|   (Sf7^0^.dfi 

Z^iibjd   -<bj6>^^    (Sf7ae)^dn  (12.50) 


To   5   $      in    (11.21)    summed   over  +   and    -  we   must    only  add 
the   term 

y    \  [dj6  -  <6/5>]^(Sf7S£)  dfi 
tJ  t  t  (12.31) 

=  j  (6/5)^C^B  dx  +  Z^   [<5j23>^(SfySe)^dfi 

For  a  sufficient  condition  we  can  drop  this  entire  positive 

term,  taking  (11.21)  as  given;  for  a  necessary  condition 

2  2 

we  can  drop  the  <6/!i>  term  and  keep  only  the  term  in  (6/))  . 
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1J>.      Local  Stability 

It  would  seem  plausible  to  expect  that  the  stability 
of  a  general  contained  plasma  with  respect  to  localized 
(or  small  wavelength)  disturbances  should  be  somehow  related 
to  stability  of  the  infinite  homogeneous  plasma  with  the 
local  value  of  the  distribution  function  and  magnetic  field. 
To  be  more  precise,  the  nonuniform  plasma  can  be  expected 
to  be  unstable  if  the  infinite  plasma  is  unstable  for  small 
wavelengths;  and  when  the  infinite  plasma  is  stable,  the 
general  plasma  should  be  stable  with  respect  to  all  suffi- 
ciently localized  disturbances. 

For  the  scalar  pressure  plasma,  this  is  a  known  result. 
The  infinite  homogeneous  plasma  is  always  stable;  it  is  in 
absolute  thermodynamic  equilibrium.   And  there  is  a  theorem, 
proved  by  H.  Rubin  [45], that  any  smooth  equilibrium  is  stable 
with  respect  to  sufficiently  localized  disturbances. 

For  the  GCP  there  is  no  natural  length  scale.   There- 
fore stability  of  the  infinite  plasma  is  independent  of 
wavelength.   It  turns  out  that  local  stability  of  a  general 
contained  plasma  is  essentially  identical  to  variational 
stability  of  the  Infinite  plasma.   We  recall  that  for  the 
infinite  plasma  variational  stability  disagrees  with  bounded- 
ness  of  solutions  with  respect  to  the  continuous  spectrum  and 
with  respect  to  monotoniclty  of  f  (e).   We  cannot  expect 
anything  different  for  the  contained  plasma,  so  we  compare 
only  variational  stability  for  the  two  cases.   Specifically, 
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the  following  conditions  (at  each  point  of  the  plasma)  are 
necessary  for  stability  in  an  arbitrary  contained  plasma 

Sf7Se  <  0 


B  /l^^  +  P2  -  Pi  >  0  (13.1) 


B^/H^  +  2p2  -  B^C^  >  0 

where  C^  is  defined  in  (10.2^1-)  and  (10. 16).   The  same  three 
conditions  are  also  sufficient  for  stability  of  every 
localized  disturbance  in  an  arbitrary  contained  plasma  if 
we  add  a  smoothness  condition;  e.g.,  bounded  V Pp  is 
adequate  [29].   The  degree  of  localization  required  for 
stability  can  be  estimated  from  the  equilibrium  parameters. 

As  immediate  corollaries  of  this  theorem  of  local 
stability,  any  short  plasma  with  tied  ends  is  stable  if 
(13-1)  is  satisfied;  and  any  long  and  only  slightly  curved 
plasma  with  tied  ends  will  also  be  stable. 

The  proof  of  the  local  stability  theorem  depends  on 
the  particularly  apposite  form  (11.21)  of  the  second 
variation  and  on  the  possibility  of  solving  the  integral 
equation  for  5/)  for  a  localized  disturbance.   It  is  shown 
that  the  first  two  terms  of  (11.21)  and  the  first  term  in 
(12.27)  dominate  in  any  local  disturbance. 

Just  as  in  the  MH  theory,  one  may  expect  a  discontinuity 
surface  to  provide  additional  sources  of  local  instability. 
But  in  the  GC  case  the  question  is  open.   It  is  not  sur- 

-  117  - 


prising  that  In  the  two-fluid  theory  the  constant  C^ 
replaces  Cp .   The  GCP  equilibrium  theory.  Sec.  9,    Is  Iden- 
tical for  one  and  for  two  fluids.   And  the  condition  for 
the  equilibrium  to  be  well  posed  has  been  given  In  terms  of 
Bpp/^B  which,  by  direct  calculation,  gives  Cp  for  one  fluid 
(10.7)  and  C^  for  two  (10.24). 

Although  this  has  not  been  proved.  It  seems  likely 
that  violation  of  (13-1)  will  result  In  non-posedness  of 
the  equations  of  motion  as  It  does  for  the  equilibrium  equa- 
tions and  as  It  does  for  motion  In  the  Infinite  uniform 
plasma.   Assuming  this,  we  can  Interpret  violation  of  (13-1) 
as  a  mlcro-lnstablllty .   The  equations  refuse  to  predict  a 
growth  rate;  we  conclude  this  must  be  related  to  omitted 
parameters  such  as  the  plasma  frequency  and  gyro  frequency. 
The  Identification  of  exactly  which  mlcro-lnstabllltles  are 
Incorporated  In  the  GCP  theory  Is  not  complete.   The  advan- 
tage of  this  theory  Is  that  It  gives  a  very  broad  criterion 
(13.1)  applicable  to  arbitrary  geometries,  whereas  all 
other  mlcro-lnstablllty  analyses  are  made  In  restricted 
geometries  allowing  at  most  small  linear  gradients.   The 
disadvantage  is  that  the  GCP  theory  yields  no  Information 
about  the  Instability  other  than  Its  presence  or  absence. 

There  Is  another  advantage  to  this  theory.   When  (13.I) 
Is  satisfied  any  Instability  requires  Interaction  between 
distant  parts  of  the  plasma.   We  are  Justified  In  calling 
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these  macro-Instabilities.   Thus  the  GCP  theory  provides 
a  clear-cut  separation  between  micro  (local)  and  macro 
(global)  instabilities.   We  repeat  that  although  it  is 
not  proved,  it  seems  very  likely  that  recourse  to  the 
equations  of  motion  would  yield  catastrophic  growth  for 
the  local  and  finite  growth  for  the  global  instabilities. 
It  is  interesting  to  note  that  micro-instabilities  appear 
even  in  the  GCF  which  is  a  totally  macroscopic  theory. 

The  significance  of  the  second  condition  (13-1)  needs 
no  elaboration;  it  is  present  in  exactly  the  same  form  in 
both  GCF  and  GCP  theories  and  is  required  for  a  well-posed 
equilibrium  formulation  as  well  as  initial  value  problem. 
The  significance  of  the  third  condition,  the  velocity 
gradient  instability,  is  deeper  and  lies  hidden  in  the 
combination  C„  of  moments  f , .   We  set 

,3^     ^2 


B^C^  =  T]p^/3P3_  (13-2) 


where  r;  is  a  dimensionless  constant  which  depends  on  f  . 
We  recall  that  for  n  =  1  this  stability  criterion  becomes 
the  same  as  in  the  GCP  theory  (6.l6)  and  for  monotone  f 
we  always  have  r]  >  1 .   Thus  local  stability  of  the  GCF  is 
a  necessary  condition  for  local  stability  of  the  GCP.   This 
is  a  special  case  of  a  general  theorem  which  states  simply 
that  GCF  stability  is  necessary  for  GCP  stability  (linear 
theory  [15],  nonlinear  [l6]). 
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The  minimum  value  tj  =  1  can  be  approximated  arbi- 
trarily closely,  but  it  is  not  attained  by  any  stable 
contained  distribution.   In  this  sense,  the  GCP  theory 
gives  a  best  possible  approximation  to  the  GCP  theory.   We 
have  T]  =  1  only  if  f  is  singular  in  m-  and  a  step-function 
in  e;  f  =  6(m.  -  |i,  )H(e   -  e),  where  (i  >  0  and  H  =  1  for 

e  <  e^,  H  =  0  for  e  >  e  •   But  for  this  distribution,  C>, 
o  o  * 

becomes  unbounded  near  the  edge  of  the  plasma;  thus  it  is 
always  unstable,  whatever  the  configuration  of  B.   It  is 
exactly  this  special  limiting  case  which  made  the  two 
constraints  (9-9)  and  (9-19)  coalesce  in  the  equilibrium 
GCF  theory.   (in  the  Infinite  homogeneous  case,  the  limiting 
distribution,  11  =  1,  is  not  unstable.) 

For  an  Isotropic  distribution,  f  (e)  independent  of  \i 
(this  cannot  be  contained  in  an  open  ended  mirror  system), 
we  have  rj  =  6  Independent  of  what  function  f   is  of  e. 
Also,  for  a  two-temperature  Maxwelllan  [f  =  a  exp(-be  -cm-)], 
which  can  be  anisotropic  and  can  be  contained  by  mirrors, 
we  also  have  r)  =  6.   The  moment  Cp  (and  therefore  C^)  is 
unbounded  if  f  (e)  has  an  infinite  slope  f  ~  (e  -e  )  ' 
on  the  left  (e  >  e^)  or  f°^  (e  -e)"'"/^  on  the  right  (e  <  e  ), 
Pig.  7-   This  Is  to  be  compared  with  the  infinite  plasma, 
where  the  stability  criterion  is  sensitive  to  bad  behavior 
of  f   only  at  V  =  0.   But  in  a  contained  plasma,  any  bad 
behavior  in  f   is  reduced  to  v  =  0  at  some  value  of  B.  This 
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Figure  7 

is  the  reason  why  the  "optimal"  distribution  with  r)  =  1  is 
allowable  in  the  infinite  case  but  not  in  any  contained 
geometry. 

These  restrictions  on  f  are  quite  relevant  prac- 
tically.  Stated  more  precisely,  instability  results  from 
a  large  gradient  of  f  at  a  given  value  of  e/li,  i.e.,  near 
a  certain  turning  point.   Almost  any  laboratory  injection 
of  plasma  will  yield  an  initial  distribution  with  discon- 
uities  of  f   or  very  large  phase  space  gradients  across 
some  critical  turning  point  B  =  e/iJ-.   Whether  the  plasma 
is  instantly  lost  or  quickly  shifts  to  a  more  congenial 
contained  situation  depends  on  the  circumstances;  both 
cases  can  be  observed  experimentally. 

Poor  values  of  ^  (i.e.,  values  in  excess  of  six)  can 
only  arise  together  with  anisotropy.   And  a  very  special 
case  of  this  velocity  gradient  instability  is  the  mirror 
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Instability  which  arises  from  large  values  of  Pp/p-,  •   But 
the  coefficient  C^  can  tiecome  unbounded  even  at  a  location 
where  Pp  =  p-,  ,  so  the  measure  of  anisotropy  which  is 
implicit  in  C^  is  a  much  more  subtle  one  than  just  the 
ratio  Pp/p-,. 

We  recall  from  (10.2^)  that  C^  can  be  eliminated  from 
the  stability  criterion  in  favor  of  a  "total"  pressure 
derivative , 


^B" 


>  0 


(15.3) 


The  indicated  derivative  is  with  B  as  the  parameter  along 
the  magnetic  line.   This  form  of  the  stability  criterion 
has  the  simple  interpretation  that  not  too  much  plasma 
can  be  placed  on  the  line  depending  on  the  strength  of  the 
mirrors  (Fig.  8).   The  maximum  amount  of  plasma  that  can  be 
placed  on  a  line  is  obtained  by  filling  in  Pp  up  to  the 
maximum  value  B  =  B-,  at  the  mirror, 


Figure  8 


-  122  - 


There  are  several  disadvantages  to  examination  of  Sp^SB 
rather  than  C^.   First  Bp^^B  is  completely  unobservable 
experimentally  whereas  the  distribution  function  is  both 
measurable  and  is  a  primary  input.   In  an  approximately 
uniform  field  Sp^^B  =  (Bp^Ss )/ (5b/Ss  )  becomes  indeter- 
minate while  the  limitation  on  the  distribution  function 
remains  viable.   Second  there  is  much  more  information 
revealed  by  examining  the  dependence  of  C^  on  f  .   And  third, 
the  quantity  Spp/SB  is  not  local.   The  derivative  Bpp/SB  at 
one  point  on  the  line  is  related  to  its  values  all  along  the 
line  implicitly,  through  the  common  distribution  function. 
For  example,  (13.^)  gives  a  legitimate  upper  estimate 
to  the  amount  of  plasma  that  can  be  contained  (taking  into 
account  that  the  field  will  be  altered  by  addition  of  plasma) 
[22];  but  it  is  not  a  very  good  estimate.   The  reason  is  that 
assigning  a  specific  function  Pp(B)  on  a  magnetic  line  is 
very  restrictive  of  the  distribution  function.   We  can  write 
the  moment  Pp  as 

P2(B)  =  B^  J   (x  -B)-^/2  C(x)dx  (13.5) 


B 

where 


C(x)  =  (2m)V2 


H^/^  f°(M.X,pL)d^i  (13.6) 


•J 

o 


But  the  integral  equation  (13-5)  can  be  inverted  uniquely 
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for  i^(x)  in  terms  of  p^(B), 

B 


'2 


C(^)  =  -^J   (B-x)-^/2  ^  (%)  dB  (13.7) 


X 

In  other  words,  not  only  does  the  distribution  function 
determine  Pp(B,a,p),  but  the  pressure  f-i;:ncblon  p2(B,a,p) 
determines  the  distribution  function  i!^(x,a,f3)  [and  greatly 
restricts  f  (m-x^I-l)].   In  particular,  the  special  function 
Pp  =  p^(a,p)  -  B  /2ii,  yields  the  explicit  distribution 

C(x)  =  ^  ((B^  -  x)V2^-^  ^    3_) 


4xB-^    8x  B-, 


8x^ 


tan-l(4  -  1)^/2J  _  ^^3^8) 


Unless  one  is  fortunate  enough  to  inject  this  precise  dis- 
tribution of  particles  on  each  line,  the  optimum  amount  of 
plasma  Indicated  by  the  formula  [Y^ A)    cannot  be  contained. 
In  plasma  stability  criteria  in  which  a  maximum  stable  amount 
of  plasma  proportional  to  the  well  depth  is  obtained, 
bp-^const.bB    [20],  [25]j  the  constant  is  more  sensitive 
to  the  distribution  function  than  to  the  geometry. 

In  the  GCF  theory,  local  instability  implies  non- 
posedness  in  an  arbitrary  nonsteady  flow.   We  conjecture 

that  the  same  is  true  of  the  GCP  local  instability  criterion. 

2  3 

In   this   case   Sp^SB  <   0  and  B  /[x     +  2p2  <  B  C^  are   no   longer 

equivalent.   It  is  evident  that  the  C^  formulation,  in  terms 

of  the  distribution  function,  will  be  the  correct  one. 

Further  details  of  this  theory  of  local  stability  will 

be   found   in  Ref .    [29] . 
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l4.   Mostly  Interchange  Stability 

An  Interchange  Is  a  special  variation  (allowable 
only  under  BCII  with  Insulated  ends)  In  which  the  dis- 
placement U  Is  equivalent  to  an  Incompressible  flow  In 
the  A  =  (a,P)  plane.   Writing  u.  for  the  velocity  In 
(a,P ) ,  we  have 

u-V  =  u.^ 

(14.1) 

The  constraint  curl  (U  x  b)  =  0  Is  equivalent  to 
(h/'b'K)  •  U   =   0.      In  an  Interchange  the  magnetic  field  Is  not 
varied  at  all  In  physical  space.   Only  the  distribution 
function  f°(e,|i.,A)  Is  redistributed  among  the  lines,  sub- 
ject  to  a  given  f  (a,^,A).   As  In  the  general  case,  there 
Is  no  loss  of  generality  In  restricting  ourselves  to  the 
pessimistic  variation.   This  means  that 

where 

j(e,^,A)  -  (  pda  =  2(2m)^/^  f  /e  -M.B(s)  ds       (1^.3) 

and  f  (J,|-L,A)  Is  carried  unchanged  by  the  motion  U   . 

As  a  special  variation,  stability  with  respect  to 
Interchanges  can  only  give  a  necessary  condition  for 
stability.   But  this  is  a  sufficiently  Important  concept 
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to  warrant  singling  It  out  as  a  definition  of  interchange 
stability.   Note  that  an  Interchange  is  a  variation  and 
not  a  motion;  it  is  not  clear  that  in  any  but  very  excep- 
tional circumstances  can  an  Interchange  be  an  actual  motion. 
In  an  Interchange,  the  Invariant  function  f  (J)  changes 
as  a  function  of  e  through  the  variation  of  the  dependence 
of  J  on  e  when  A  is  varied.   In  a  one -fluid  formulation 

J  {e  ,[i  ,'h)   depends  only  on  the  fixed  magnetic  field  geometry 

1/2 
and  not  on  f ;  more  precisely,  j/|-l  ^       is  a  function  of  the 

single  variable  e/ia-.   But  in  the  two-fluid  theory,  the 
potential  j6    (contained  in  e  )  depends  on  f  as  well  as  on 
the  magnetic  field.   Thus  j(e,M.,A)  depends  on  f  and  not 
only  on  the  field  geometry.   Almost  all  the  simplification 
of  an  interchange  analysis  over  an  analysis  of  the  full 
second  variation  is  lost  in  the  two-fluid  theory. 
Almost  all  interchange  results  are  restricted  to  the  one- 
fluld  theory,  although  this  is  not  usually  recognized  in  the 
literature.   We  do  obtain  a  sufficient  condition  for  inter- 
change stability  by  ignoring  ^   in  the  two-fluid  theory. 
This  only  has  significance  if  we  take  Interchange  stability 
as  an  Independently  defined  concept;  otherwise  a  sufficient 
condition  for  another  condition  which  is  only  necessary  for 
stability  yields  no  definite  conclusion. 

The  analytic  manipulations  (for  finite  Interchanges 
as  well  as  infinitesimal)  are  most  easily  carried  out  In 
the  form  [30] 
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J  =  /  e(J,^,A)f  (j,n,A^)dJdp.dadf3  (l^.^) 

where  1\   Is  an  Eulerlan  variable  and  A   Is  a  Lagranglan 
variable;  both  e  and  f  are  fixed  functions  of  the 
Indicated  variables  and  dadp  =  da  dp  . 

One  Important  role  played  by  Interchange  stability  Is 
In  low  pressure  plasmas.   It  Is  frequently  stated,  even 
"proved",  that  at  low  pressure  Interchange  stability  Implies 
absolute  stability.   The  argument  In  support  of  the  Inter- 
change hypothesis  Is  as  follows.   The  second  variation  Is 
the  sum  of  a  magnetic  contribution  (non-negative  since  the 
vacuum  field  is  an  absolute  minimum)  and  a  small  plasma 
contribution.   For  any  given  variation  U  which  is  not  an 
Interchange,  the  plasma  contribution  will  be  dominated  by 
the  positive  magnetic  contribution  below  some  (3   (we  can 
consider  p  to  be  a  parameter  multiplying  f).   But  an  inter- 
change gives  a  positive  second  variation  by  hypothesis. 
Therefore  the  second  variation  is  positive  for  sufficiently 
small  p.   The  fallacy,  of  course,  is  that  the  critical  3 
depends  on  the  variation  U,  and  p   may  have  the  lower  bound 
zero  as  U  is  varied  over  the  admissible  class. 

The  important  question  is,  of  course,  not  whether  the 
argument  is  false  but  whether  the  conclusion  is  false. 
But  it  is  so  demonstrably  false,  with  so  many  very  familiar 
theoretical  and  experimental  counterexamples,  that  it  is 
an  interesting  psychological  question  how  this  belief  has 
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come  to  gain  such  wide  acceptance.   Experimentally,  any 
mlcro-instablllty  (well  known  to  exist  at  low  p)  is  an 
Immediate  counterexample.   Perhaps  the  reason  these 
experimental  "counterexamples"  have  not  been  faced  In 
this  context  Is  the  lack  of  appreciation  that  the  GCP  Is 
a  microscopic  theory.   But  the  belief  is  strong  enough  so 
that  when  one  author  is  faced  with  an  immediate  theoretical 
counterexample  to  his  "proof"  of  sufficiency,  he  resolves 
the  conflict  by  pragmatically  strengthening  the  "sufficient" 
condition.' 

In  the  GGF  and  the  GOP  it  is  very  easy  to  present 
counterexamples,  and  in  the  MH  theory,  more  subtle  ones. 
For  example,  one  must  always  require  local  stability,  since 
this  is  a  universal  necessary  condition.   Furthermore, 
local  stability  is  independent  of  interchange  stability  and 
is  not  consequent  upon  low  pressure.   In  the  GCF  we  can 
easily  exhibit  equilibria  in  which  Pp/p-^  is  unbounded  near 
the  mirror.   No  matter  how  small  is  the  value  of  p ,  local 
stability  will  be  violated  near  the  mirror.   In  the  GCP 

version  the  same  mirror  instability  is  excited  if  f   goes 

1/2 
to  zero  slower  than  (e.,  -  e)  '^   at  the  edge  of  the  plasma 

(e-j  =  I-lB^  ) .   But  we  have  much  more  flexibility  than  this  in 

the  GCP  theory;  the  velocity  gradient  instability  can  be 

excited  in  the  plasma  Interior  whenever  f  (e)  has  a  large 

gradient,  and  this  can  produce  instability  at  arbitrarily 
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low  p  and  with  bounded  p  /p, •   It  Is  even  possible  to 
give  unstable  counterexamples  where  VPo  ^-^^  curl  B 
go  to  zero  with  p.   And  there  still  remains  the  local 
instability  associated  with  f  not  monotone  in  e  which 
(especially  in  the  form  of  electrostatic  oscillations  as 
suggested  in  Sec.  10)  remains  at  arbitrarily  low  p.   These 
gaps  are  not  only  mathematical;  there  are  very  many  experi- 
mental instabilities  observed  at  very  low  p  and  under 
circumstances  which  are  probably  compatible  with  inter- 
change stability  (or  the  equivalent,  tied  ends).   Although 
the  experimental  identification  of  instabilities  is  not 
yet  a  fully  developed  discipline,  it  is  very  likely  that 
the  GOP  local  instabilities  lie  among  those  observed. 

A  further  gap  in  the  interchange  hypothesis  lies  in 
the  possibility  of  non-smooth  equilibria.   Even  in  the  MH 
theory,  a  shear  surface  across  which  the  tangential  com- 
ponent of  B  is  discontinuous  can  be  unstable  at  arbitrarily 
low  p.   We  recall  that  sufficiency  for  local  stability  in 
the  GGP  is  more  demanding  than  necessity  with  regard  to 
smoothness . 

Although  the  simple  form  of  the  Interchange  hypothesis 
is  false,  it  can  be  replaced  by  a  slightly  more  sophisticated 
version  which  is  much  more  plausible.   Local  stability 
involves  small  scale  disturbances.   Interchange  stability 
involves,  in  a  sertain  sense,  the  largest  scale  disturbances; 
the  scale  length  is  not  even  restricted  by  the  apparatus 
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dimension,  since  the  ends  are  free  to  move.   Tying  down 
stability  at  both  extremes  In  wavelength  would  seem  to 
give  a  good  rough  estimate  with  regard  to  absolute  stability. 
We  conjecture  that  the  combination  of  local  and  Interchange 
stability  Implies  absolute  stability  under  fairly  wide 
circumstances.   For  example,  there  is  plausible  evidence 
that  this  is  so  under  the  further  assumption  of  low  p,  or 
a  long  thin  device,  or  a  device  with  limited  curvature  to 
the  magnetic  lines,  etc.   The  hypothesis  of  interchange 
stability  can  be  dropped  in  a  system  with  tied  ends.   We 
have  already  given  several  examples  of  absolute  stability 
contingent  on  local  stability  and  tied  ends  (Sec.  13)-   By 
far  the  largest  and  most  diversified  part  of  the  current 
knowledge  of  GCP  stability  lies  in  this  combination  of 
interchange  and  local  stability  [30]. 

We  interpolate  a  brief  Introduction  describing  inter- 
change stability  for  the  scalar  pressure  MH  plasma.   We  have 
already  mentioned  the  restriction  that  p  be  a  function  of 
q  =  /  ds/B  in  order  for  the  equilibrium  to  be  compatible 
with  free  ends.   These  are  the  only  equilibria  for  which 
we  can  study  interchanges.    A  simple  analysis  gives  as  a 

necessary  and  sufficient  condition  for  stability  with 

Y 
respect  to  finite  or  infinitesimal  interchanges  that  pq 


This  restriction  on  interchange  equilibria  is  given 
in  [^2].   The  same  formula  for  a  torus  (at  low  p)  was  found 
earlier  [43].   In  a  torus  this  is  no  restriction;  it  is  a 
property  of  all  equilibria. 


130  - 


DO  a  monotone  function  of  q.    Writing  p  ==  p(q),  the 
criterion  is 

p'  +  yp/q  >  0  .  (1^-5) 

If  p  =  V^if)    and  q  =  Q('^)  are  expressed  as  functions  of  a 
flux  parameter,  this  takes  the  form 

(^+  T^)Q'  >  0  (U.6) 

(there  is  no  restriction  to  axial  symmetry  implicit  in 
this  description  in  terms  of  i) .      An  obvious  sufficient 
condition  is 

p'  >  0  ,   P'Q'  >  0  (1^.7) 

or  that  p  be  a  monotone  function  of  q.   For  a  contained 
plasma  In  which  p  is  monotone  in  physical  space,  P'  <  0, 
this  requirement  is  that  q  be  a  maximum  at  the  center  of 
the  plasma.   We  repeat  that  this  universally  cherished 
condition,  that  /  ds/B  be  a  maximum  within  the  plasma.  Is 
mathematically  neither  necessary  nor  sufficient  for  absolute 
stability.   It  Is  sufficient  for  interchange  stability  and 
therefore  may  be  sufficient  for  stability  at  low  p  or  in 
other  suitably  restricted  circumstances.   At  the  edge  of 


The  stability  criterion  (14.5)  is  also  given  in  [^5] 
for  low  p,  toroidal  equilibria  and  in  [^2]  for  general 
equilibria  and  finite  perturbations. 
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the  plasma,  if  both  p  and  p'  approach  zero,  the  more 
stringent  sufficient  condition  p'(q)  >  0  approaches  the 
necessary  and  sufficient  condition  (l4.5).   But  it  is  an 
unnecessarily  strong  restriction  within  the  plasma,  and 
it  is  not  required  for  a  plasma  which  terminates  at  a 
finite  gradient  (we  return  to  this  point  in  the  GCP  analysis) 

The  scalar  pressure  theory  is  applicable  only  to  a 
toroidal  system  or  one  where  the  plasma  ends  at  material 
walls.   In  the  latter  case,  the  shape  of  the  walls  exerts  a 
strong  influence  on  the  variation  of  q  from  line  to  line, 
and  therefore  on  stability.   We  note  that  GCP  stability  for 
an  isotropic  equilibrium  is  not  the  same  as  MH  stability 
even  though  the  equilibrium  theory  is  identical.   And  we 
reiterate  that  even  in  the  MH  case  where  it  is  quite 
plausible,  it  has  not  been  proved  that  interchange  stability 
implies  stability  at  low  p. 

We  now  turn  to  GCP  interchange  stability.   This  theory 
is  more  complex  than  the  MH  theory  in  that  an  entire  dis- 
tribution function  rather  than  a  number,  p,  is  assigned  to 
each  line.   No  simple  universal  necessary  and  sufficient 
condition  for  interchange  stability  such  as  (1^.5)  is  avail- 
able.  But  we  can  completely  analyze  several  special  classes 
of  distribution  functions  which  taken  together  seem  to  be 
wide  enough  to  give  a  reliable  evaluation  of  the  general 
situation,  and  we  shall  be  able  to  supplement  this  with  some 
sufficient  conditions,  not  restrictive  of  distribution 
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function,  which  support  the  same  general  conclusions.  A 
crucial  point  will  be  the  assessment  in  each  case  of  how 
representative  the  results  may  be. 

First  we  remark  that  the  special  case 

f  =  f*(J,M-)  ,  (1^-8) 

i.e.,  the  same  function  of •  J  and  I-l  on  each  magnetic  line, 
is  completely  neutral  to  all  interchanges.   This  is  evident, 
since  it  is  the  function  f  (J,i-i.)  which  is  carried  unchanged 
from  line  to  line  by  the  pessimistic  variation.   This  is  a 
mathematical  rather  than  a  practical  example  since  f  (J,|-l) 
cannot  be  contained  in  any  magnetic  field.   But  since  it  is 
neutral,  it  is  a  useful  basis  for  perturbations  and  for 
the  determination  of  a  transition  between  stability  and 
instability. 

The  next  special  class  of  distributions  is 

f  =  a(A)f(j,n)  ;  (1^.9) 

the  distribution  is  fixed  in  J  and  M-  except  for  the  amount 
placed  on  each  line.   This  case  can  be  examined  in  complete 
detail,  for  finite  as  well  as  infinitesimal  stability,  just 
as  in  MH.   The  condition  that  the  equilibrium  be  compatible 
with  free  ends  (j  =  O)  is  that  o(A)  be  a  function  of 

e(^)  =  /  e(J,M-,A)f (j,^L)dJdii 

(14.10) 

r   /I   ^     ,   '^   N   dS 

=  /  (2  Pi  +  P2)  -B"  ' 
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I.e.,  a  =  const,  must  be  placed  on  e(A)  =   const. 

A  necessary  and  sufficient  condition  for  finite  inter- 
change stability  is  that  a  be  a  monotone  decreasing  function 
of  £.   This  result  is  very  similar  to  the  scalar  pressure 
sufficient  condition,  p  '^  q,  but  the  analogue  is  only 
mathematical  since  the  parameters  a   and  e  are  quite  unrelated 
to  p  and  q. 

For  this  special  class  of  equilibria  there  always  exist 
surfaces,  e  =  const.,  which  are  neutral  with  respect  to 
interchanges.   Taking  ii   to  be  the  flux  parameter  on  the 
neutral  surfaces,  we  find 

■P2  SB  .  ,B^  .  __     „  x-1  _.   ^P* 

L 
O 


/,,,N    {.^2    Sb  ,  /B^  ,         x-1      P*  n  ds  ,  „    ,,,,  ,,x 
(^)  =  h-r  ^+  (—  +  P2  -  Pi)    Pl3^^^>°    (1^-11) 

<J  o 


as  the  criterion  for  stability  [based  on  a  monotone  profile 
a'('i/')  <  0].   The  derivatives  h/hii   are  to  be  taken  in  the 
direction  normal  to  the  ii   surface.   For  a  low  p  plasma 
(vacuum  or  force-free  field)  (l4.1l)  simplifies  to 


[ 


(Pl  +  P2)  If  %>  0  (1^-12) 

B 


which  is  a  simple  anisotropic  modification  of  the  maximum 
/  ds/B  criterion  since,  for  low  p,  one  can  verify  that 


Sq  _  _p  f  Sb  ds  n4  1^1 


We  shall  see  in  a  moment  that  the  validity  of  (l4.12)  depends 
on  the  precise  definition  of  low  p. 
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The  stability  criterion  (l4.12)  Is  very  frequently 
quoted,  and  was  originally  found  by  Rosenbluth  and 
Longmlre  [  9  ]  using  heuristic  arguments  which  hid  the  fact 
that  the  result  depends  on  a  very  special  assumption  with 
regard  to  the  distribution  function.   We  shall  also  find 
that  this  criterion  is  unnecessarily  pessimistic.   For 
example,  it  is  easily  verified  that  the  criterion  (1^.12) 
can  never  be  satisfied  near  the  axis  in  an  axially  symmetric 
mirror  field.   But  the  finite  pressure  formula  (l4.1l)  is 
much  more  easily  satisfied.   For  example,  in  any  plasma 
with  fixed  p,  increasing  the  length  to  diameter  ratio  will 
eventually  make  it  interchange  stable.   The  criterion  is 
roughly  that  stability  is  achieved  if  p  exceeds  some 
measure  of  the  line  curvature  compared  to  plasma  radius. 

If  both  a  and  a'(i/')  (i.e.,  p^  and  V  Po )  approach 
zero  at  the  edge  of  the  plasma,  then  the  low  p  formula 
(1^.12)  governs  interchange  stability  at  the  edge  of  the 
plasma.   This  condition  is  very  difficult  to  satisfy  without 
special  stabilizing  windings.   But  the  situation  is  quite 
different  if  a  approaches  zero  linearly  In  i/ .      In  this 
case  the  second  variation  must  be  carefully  evaluated 
together  with  boundary  terms  at  the  interface  between  plasma 
and  vacuum  which  represent  interchanges  between  vacuum 
lines  and  lines  bearing  plasma.   The  same  formula  (l4.1l) 
is  obtained  but  with  the  interpretation 


-  155  - 


,P2  hB    ,    Pi   ^P*^  ds 

,2 


o 
We  have  dropped  p^  and  Pp  compared  to  B^/p-  where  they  are 
not  differentiated.   The  derivatives  hB/'di'   and  Sp^Si/'  are 
evaluated  for  the  actual  equilibrium,  but  p.,  and  Pp  are 
evaluated  for  the  representative  function  f;  in  particular, 
P-,  and  Pp  are  not  zero  even  in  the  vacuum.   The  function 
e{i')    has  a  discontinuous  slope  Just  as  does  o{i').      The 
conclusion  is  that  the  plasma  is  stable  if  e'(i^)  >  0  on 
the  plasma  side;  a  change  in  sign  of  e^  (i/)    across  the  inter- 
face is  irrelevant  to  the  sign  of  the  second  variation. 
This  means  that  with  finite  f3  (i.e.,  finite  gradients  of 
pressure)  we  can  easily  stabilize  the  entire  plasma  out  to 
the  edge . 

But  the  problem  of  interchange  stability  with  the 
special  distribution  (l4.9)  is  amenable  to  nonlinear 
analysis.   If  there  is  a  change  in  sign  of  e'('iA)  across 
the  interface,  the  plasma  is  not  stable  to  finite  inter- 
changes.  Of  course,  no  plasma  is  stable  with  respect  to 
arbitrary  finite  perturbations;  the  most  stable  location 
for  the  plasma  is  at  infinity.   But  in  the  present  problem, 
the  amplitude  of  the  displacement  which  is  required  for  a 
release  of  potential  energy  is  proportional  to  the  distance 
of  the  magnetic  line  from  the  plasma  edge.   It  requires  an 
arbitrarily  small  interchange  close  to  the  edge.   With  a 
finite  amount  of  noise,  a  certain  thin  layer  will  be  unstable, 
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But  the  energy  released  is  also  arbitrarily  small.   It  is, 
for  example,  of  too  high  an  order  to  be  visible  in  the 
stable  second  variation.   A  simple  yes  or  no  answer  to  the 
question  of  stability  may  be  misleading.   Whether  the  linear 
or  nonlinear  conclusion  is  more  reliable  could  depend  on  the 
noise  level . 

One  can  also  point  out  that  any  smoothing  of  the 
profile  oiif)    to  make  a'(i^)  continuous  will  create  a  tech- 
nically unstable  situation  even  in  the  linear  analysis. 
But  every  stable  plasma  is  arbitrarily  close  to  another 
unstable  one  [e.g.,  by  introducing  an  arbitrarily  small 
change  in  f  (e)].   Also  the  instability  created  by  smoothing 
a(i^)  will  be  very  weak  as  well  as  localized. 

To  summarize  the  situation,  using  the  same  (linear) 
theory  as  in  almost  all  other  stability  analyses,  this 
configuration  is  stable.   In  a  practical  sense,  one  has 
reason  to  expect  a  weakly  excited  disturbance  in  a  thin 
layer  which  depends  on  the  noise  level  within  the  plasma. 
As  in  any  analysis  which  is  sensitive  to  a  small  region, 
the  conclusion  (and  the  instability)  may  be  wiped  out  by 
finite  gyro  radius  effects.   It  is  certainly  clear  that 
interchange  stability  is  enhanced  by  increasing  p  and  by 
a  sharp  rather  than  a  gradual  transition  at  the  edge. 

Next  we  turn  to  a  somewhat  more  flexible  distribution 
function 

f  =  af(j/T,  n)  =  af(j,^L)  (1^.15) 
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with  two  scaling  parameters  a(A)  and  t(?\).   This  time  we 
cannot  do  the  analysis  In  the  large  and  we  examine  the 
first  and  second  variations.   The  pessimistic  Interchange 
variations  are  given  by 

6|  =  -  [eC^-lx^  dJdM-dadp  (l4.l6) 

6^^  =.  -  [  (u..  ||-)(cc.|i)  dJdiidadp         (1^.17) 


where  S/S>v  Is  understood  to  imply  that  J  and  M-  are  fixed. 
For  simplicity  we  assume  that  the  contours  a(7\)  =  const, 
agree  with  t(7\)  =  const.   There  is  a  common  flux  parameter 
if    (unspecified)  such  that  o   =  o{i/),    t  =  t(-^).   This 
specialization  implies  the  existence  of  a  family  of  sur- 
faces {if   =  const.)  which  are  neutral  to  Interchanges   but 
is  not  restricted  restricted  to  axial  symmetry.   The  condi- 
tion that  the  first  variation  vanish  (equivalent  to  BCIl) 
is  that 

be  a  function  of  i/ ,   where 

£^(/\,t)  -  /  e(Tj,ii,A)f (j,^L)djd^i 
e2(^^'^)  =  -/  e(Tj,H,A)[j  Sf/^j]djd^i  . 


(1^.19) 


We  may,  for  example,  suppose  that  o{i')    and  t(i^)  are  given 
functions.   For  a  given  value  if  =  i/     vie   examine  the  contours 
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P(A,t/^  )  =  constant  in  the  A-plane  and,  assuming  that  they 

are  concentric,  choose  the  one  with  the  flux  value  i/     and 

o 

assign  to  It  the  value  if    .     Varying  i^  In  p{'\,i'    )   gives  a 
family  of  p -contours  each  with  a  value  i/    ,    thereby  specifying 
a  function  'i/'(a,p).  This  construction  can  be  carried  out 
provided  that  the  Jacoblan 


a(p,^)  _  a-  ^(^1-^)  ,  T-  ^(f21^ 


a(a,p)  -  ~  ^(a,p)   +  ~  a(a,f3)  (l4.20) 

does  not  vanish  (note  that  the  dependence  on  t  of  e.,  and  e^ 
can  be  Ignored  In  these  Jacoblans). 

The  second  variation  stability  condition  can  be  shown 
to  reduce  to 

—  e^  +  —  e,,  <  0  (14.21) 


where 


e^(VO  -  j  |§  f  cijdn 


(14.22) 

^4(*)  =  -jsfJsf -^J-i^^ 

by  using  the  property  that  p  In  (l4.l8)  Is  a  function  of 
i/   alone.   (Note  that  e^  and  e^   are  not  i/^ -derivatives  of  e-, 
and  Ep . ) 

We  are  primarily  Interested  In  equilibria  In  which 
a'(V')  <  0;  t(i/')  Is  at  our  disposal.   The  stability  condition 
e-,  <  0  which  results  for  t  :=  1  Is  exactly  the  same  as  the 
previous  result  (l4.1l).   In  an  axlally  symmetric  low  p 
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field,  ^e/hif   will  be  negative  at  the  center  of  a  magnetic 
line  and  positive  towards  the  mirrors.   Since  f  concen- 
trates Its  largest  values  at  the  center,  the  central 
negative  values  of  he/hif   tend  to  dominate,  and  (at  low  p) 
the  e^  term  will  frequently  be  destabilizing.   But  In  Su, 
f  Is  replaced  by  -j  Sf/Sj  =  -J  5f/Sj  which  concentrates 
more  at  the  large  values  of  J;  f  must  be  monotone  In  J 
for  local  stability,  but  this  Is  not  required  of  J  hf/bj. 
Thus  It  turns  out  to  be  possible  to  stabilize  even  an 
axlally  symmetric  vacuum  field  by  making  the  distribution 
f (J/t)  narrower  In  J  as  one  moves  out  radially.   On  the 
other  hand,  this  process  cannot  easily  be  carried  to  the 
edge  of  the  plasma  where  a^/o    Is  unbounded,  since  t  cannot 
be  decreased  below  a  certain  value  without  changing  the 
sign  of  £k  (one  can  arrange  to  shape  the  field  so  that  t 
can  approach  zero  with  a,  but  this  requires  quite  special 
fields) . 

In  summary,  we  have  described  two  stabilization 
mechanisms,  finite  (3  (which  can  even  be  carried  to  the 
edge  of  the  plasma)  and  scaling  f  to  be  narrower  outside 
(which  works  even  at  zero  p,  but  not  close  to  the  edge). 

We  now  turn  to  a  different  series  of  special  distribu- 
tions not  based  on  the  Interchange  neutrality  of  the  function 
f  =  f*(j,ii.)  .   Consider 

f(j,^i,A)  =  f°(e(j,H,A),n)  .  (14.23) 
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In  other  words,  place  the  same  distribution  of  e  and  I-l 
on  each  magnetic  line.   The  distribution  f  (J,iJ.)  cannot 
be  contained;  f  (e,[-L)  is  also  severely  limited,  but  there 
do  exist  special  containing  fields. 

The  local  stability  criterion,  Sf  /he   <  0,  is  seen 
by  inspection  to  yield  a  positive  second  variation  in 
(1^.17),  Sf/ax  =  {hf°/he)ihe/h'K)  .      All  special  distribu- 
tions f  (e,M-)  which  are  locally  stable  are  automatically 
interchange  stable.   The  condition  for  interchange  stability 
alone  is  weaker,  viz.,  that  the  2X2  matrix 

r  Se   Se   Sf° 


A,, (A)  - 


3X7  3X7  3F 

1   J 


dJdix  (14.24) 


be  negative  definite  on  each  line,  or  the  equivalent,  that 


1 


(|f)  i^dJd^  (1^.25) 


be  negative  for  every  directional  derivative  b/'bi/   in  the 
A-plane . 

Directly  from  their  definitions,  p. (b)  and  Pp(B)  are 
seen  to  be  independent  of  A.   For  confinement,  with  concen- 
tric pressure  contours  reaching  zero  at  the  outer  confines 
of  the  plasma,  we  must  have  concentric  contours  of  constant 
B  within  a  minimum-B  or  magnetic  well  field  configuration. 
Note  that  this  field  configuration  is  required  for  contain- 
ment, not  for  stability  which  is  automatic. 

The  next  example  in  this  sequence  is  the  more  general 


-  141  - 


function 

f(j,^i,A)  =3  f^(£(J,^i,?\),J,^L)  (14.26) 

This  class  of  functions  Is  too  general  to  yield  a  simple 

theory  of  Interchange  stability.   But,  as  In  the  previous 

example,  an  evident  sufficient  condition  for  interchange 
stability  Is 

T^  <  0  .  (14.27) 


And,  as  before,  the  weaker  (but  complicated)  condition 
(l4.24)  or  (14.25)  Is  the  more  precise  condition  for  Inter- 
change stability.   But  this  time  the  local  stability 
condition  Is  entirely  Independent  of  the  interchange 
condition.   Setting 

f°(e,[i,A)  =  f^(E,j(£,lx,A),^L)  (14.28) 

we  have 

^„o        hf-.         Sf-i  A-. 

for  local  stability.   The  factor  Sj/5e   (the  orbital  period) 
is  positive,  but  Sf-,/Sj  is  Indeterminate.   A  sufficient 
condition  for  local  and  interchange  stability  would  be 

f-<  0   ,  ^<  0  (14.30) 

together  with  the  inequalities  (13-1) • 
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The  class  of  distributions  (l4.26)  subject  to  (1^.30) 
Is  subject  to  the  same  restriction  on  containment  as 
f  (e,M-)j  viz.,  to  minimum  B.   If  we  denote  by  B-,  (J,ii-,A)  -  z/\i 
the  turning  point  specified  by  the  given  parameters,  we  see 
that  f-,  (i-lB-,  ,J,i-l)  is  constant  on  magnetic  contours  B-,  =  const, 
for  fixed  J  and  I-l.   If  f-,  is  monotone  in  e  =  M-B-,  ,  it  is 
monotone  in  B-,  for  fixed  J  and  \i.      For  containment  the  contours 
of  B-,  must  be  closed  and  monotone. 

The  significance  of  the  Interchange  sufficient  condition 
Sf-,/Se  <  0  can  be  seen  by  Introducing  a  particle  interchange 
(up  to  now  we  have  considered  line  interchanges).   Suppose 
we  consider  fixed  (J,M-)  in  f(j,|J.,A)  and  interchange  the  values 
of  f(A)  Incompressibly  in  the  A-plane  to  minimize  /  e  f dadp . 
The  minimum  is  obtained  by  transferring  the  largest  f-values 
to  smallest  e;  in  other  words,  f  is  a  monotone  decreasing 
function  of  £(A)  for  the  given  (j,[i),  Sf  /Se  <  0.   This  particle 
Interchange  violates  the  constraint  that  all  particles  on  a 
line  interchange  together.   The  minimum  is  too  low  and 
Sf-,/Se  <  0  is  sufficient  but  not  necessary.   The  stringent 
condition  Sf-,/Se  <  0  is  violated  by  most  of  our  previous 
examples  of  Interchange  stable  systems.   It  appears  to  be 
unnecessarily  strong  for  stability  and  unnecessarily  restric- 
tive of  containment.   But  this  criterion  would  seem  to  be  the 
one  that  should  arise  in  a  consistent  higher  order  "drift" 
GC  theory  in  which  particles  are  allowed  to  move  from  line  to 
line  [23],  [21]. 
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The  main  reason  for  introducing  the  special  class 
f-,(e,J,M-)  Is  not  its  simplicity  but  its  significance  for 
long  time  containment.   Some  elements  of  a  theory  of 
equilibrium  can  be  carried  out  for  a  long  time  scale  0(e). 
In  a  static  field  the  guiding  centers  will  drift  away  from 
magnetic  lines  on  this  time  scale,  but  will  remain  on  drift 
surfaces  J  =  const.   Thus  if  there  exists  an  unchanging 
equilibrium  configuration  over  this  time  scale,  the  distribu- 
tion function  must  take  the  form  f-,  (e,ii., j)  .   This  automatically 
satisfies  the  interchange  boundary  condition,  J  =0,  but  not 
conversely. 

There  are  several  reasons  for  carrying  out  the  stability 
theory  for  the  more  general  fast  equilibria  f  {e,\x,')\).      First 
there  is  no  apparent  simplification  in  taking  f-,,  and  all 
general  results  proved  for  f  apply  also  to  f.,  .   Secondly, 
no  set  of  equations  of  motion  has  yet  been  formulated  which 
is  accurate  on  both  the  short  and  long  time  scale.   There 
are  higher  order  drift  equations  which  uncover  certain  slow 
micro-instabilities  not  adequately  described  by  the  lowest 
order  GC  theory.   But  all  long  time  formulations  assume 
that  J  is  constant  and  are  thus  Incapable  of  describing  fast 
motions.   In  particular,  it  is  pointless  to  use  slow  "drift" 
equations  to  investigate  interchanges  which  are  strictly 
fast  instabilities.   As  we  have  already  mentioned  in  Sec.  2, 
the  treatment  of  interchanges  using  higher  order  drifts 
uses  an  energy  minimization  for  interchanges  which  properly 
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belongs  to  the  zero  order  GC  theory.   Finally  there  is 
the  important  practical  question  of  the  behavior  of  an 
actual  system  which  has  been  created  close  to  but  not 
exactly  in  equilibrium.   Assuming  that  there  is  an  approach 
to  a  long  term  equilibrium  f.  (e,J,ii.)  of  an  incorrectly 
created  plasma,  the  approach  will  be  on  the  long  time  scale. 
The  stability  of  the  fast  equilibria  met  along  the  way  is 
just  as  important  as  the  stability  of  the  ultimate 
equilibrium  state. 

The  relatively  complete  (necessary  and  sufficient) 
Information  which  we  have  described  for  interchange  stability 
of  the  special  classes  of  functions  af{j/x,[i)    and  f  (e,|-L), 
can  be  supplemented  by  a  sufficient  condition  which  is  valid 
quite  generally  [30], 


r  r  2/  ^      .-1  ^P*  ^Pi  ^  1  SB  ^P2  : 
J  I  (B  /^.^  +  P2  -  p^j     ^r~  w~  ^  bm  d^  ' 


[^  <  0  .   (14.31) 


Here  ^/hi/   is  a  derivative  normal  to  the  line  [as  in  (l4.ll)] 
and  d/di^  is  a  derivative  holding  B  fixed  [as  when  Pp(B,''^) 
is  given  in  equilibrium] 

If  there  are  neutral  surfaces,  i/   =   const.,  the  formula  is 
unambiguous.   If  not,  b/hif   is  a  directional  derivative  and 
(l4.3l)  must  be  satisfied  for  all  directions  (or  we  can 
write  the  condition  as  a  2  X  2  matrix  involving  derivatives 
h/^   and  h/b^)  . 
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We  confirm  qualitatively  from  this  sufficient  condition 
the  stabilizing  effect  of  finite  p,  of  slight  field  curva- 
ture, and  of  scaling  the  distribution  to  be  narrower  on  the 
outside.   For  example,  if  we  let  the  geometry  become  long 
and  straight  at  any  finite  p  (the  mirror  ratio  can  be  finite, 
so  long  as  the  curvature  is  small),  then  'dp^hif   becomes 
small  and  Sb/SV'  becomes  positive;  the  sufficiency  condition 
will  always  be  satisfied. 

There  are  a  few  results  which  have  been  proved  using 
the  full  second  variation  as  distinguished  from  interchange 
stability.   First  we  recall  the  stability  results  with  tied 
ends  that  follow  from  local  stability  (Sec.  13)-   With  ends 
free,  Taylor  and  Hastie  [22]  have  shown  that  for  the  special 
distribution  f°(e,|-L)  local  stability  is  not  only  necessary 
but  is  also  sufficient  for  stability  in  the  one-fluid  case. 
We  recall  that  in  this  special  case  local  stability 
trivially  implies  interchange  stability  (the  partial 
derivatives  Sf°/Se  and  Sf-,/B£  coalesce).   In  another  direc- 
tion, Taylor  and  Hastie  carry  out  an  ingenious  singular 
expansion  about  a  uniform  magnetic  field  and  low  p  in  which 
p  /p   and  Sf  /he   remain  finite  as  p^ ,  p^ j  and  f  become 
small.    Limiting  forms  of  the  local  necessary  conditions 
(for  low  p  and  one  fluid)  are  obtained,  also  a  sufficient 


The  result  is  claimed  to  hold  for  finite  P;  this  is 
somewhat  exaggerated  since  it  is  the  leading  term  in  an 
expansion  in  p . 
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condition  which  can  be  written 


1 


or 


/ 


(b'o/^  -  V^o)"'  ^  •  tS^-  (V^o- VBo)  ^J^-   <  0 


(14.34) 
Implicit  in  this  form  is  the  existence  of  neutral  surfaces; 
the  normal  derivatives  are  normal  to  the  unperturbed  uni- 
form field.   This  formula  is  quite  similar  to  (l4.3l)  but 

is  disturbingly  different.   To  compare  the  two  we  insert 

2 
p  «  p  «  B  and  the  perturbation  about  uniform  B   into 

(l4.3l)  [Sp^Sxj_  ~  Bb/Sx_|_  '-^  5pp/Sxj_  are  equivalent  in 

order],  to  obtain 

Examination  shows  that  (14.34)  is  a  much  more  stringent 
condition.   It  does  not  at  all  predict  the  finite  p 
improvement  in  stability.   It  is  very  difficult  to  satisfy 
(14.34)  without  a  magnetic  well;  (14.35)  can  be  satisfied 
in  every  long  thin  geometry,  almost  independent  of  the 
distribution  function  and  field  configuration.   Strictly 
speaking  the  two  criteria  are  mathematically  independent; 
(14.34)  is  a  sufficient  condition  for  absolute  stability, 
and  (14.35)  guarantees  interchange  and  local  stability. 
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which  are  only  necessary  for  absolute  stability.   Never- 
theless, barring  errors  In  computation,  we  must  conclude 
that  either  local  plus  Interchange  stability  at  low  p 
is  not  sufficient  for  stability,  or  the  sufficient 
estimate  given  by  Taylor  and  Hastie  Is  very  crude. 

We  conclude  with  a  qualitative  discussion  of  the 
practical  implications  with  regard  to  stability  of  the 
presently  available  GCP  theory.   The  possibilities  are  much 
more  elaborate  than  in  MH.   In  the  latter,  only  the  magnetic 
field  configuration  is  subject  to  significant  variation. 
In  the  GCP  there  Is  also  the  distribution  function. 
Stability  requires  the  proper  correlation  of  a  distribution 
function  with  a  magre  tic  field,  and  there  appear  to  be  many 
such  stable  combinations. 

Any  theoretical  analysis  is  severely  constrained  by 
the  criterion  of  mathematical  simplicity.   A  primary  prob- 
lem is  to  assess  whether  models  which  are  chosen  for  their 
accessibility  to  analysis  are  representative  of  the  general 
case.   The  simplest  choice  of  distribution  is  the  absolute 
Maxwelllan.   This  is  completely  stable  under  all  circum- 
stances, but  it  cannot  be  contained  by  any  magnetic  field; 
pressure  is  constant  everywhere.   Next  in  simplicity  is 
probably  the  isotropic  distribution  f  (e)  or  f  (e,A).   The 
former  Is  relatively  easily  stabilized,  but  like  the 
Maxwelllan  it  cannot  be  contained.   The  latter  can  be  con- 
tained only  in  a  toroidal  geometry,  and  we  shall  therefore 
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disregard  it  In  the  present  context.   The  simplest  dis- 
tribution that  can  be  contained  In  any  open-ended  system 
Is  f  =  f  (e,iJ.).   This  turns  out  to  be  relatively  easily 
stabilized  but  again  Is  quite  restricted  In  the  containing 
fields  which  cannot  be  other  than  magnetic  wells. 

But  the  special  assumption  f  {e,\i)    [and  the  more 
general  f-,  (e,i^,j)  with  5f-,/Se  <  0]  offers  no  clue  In 
distinguishing  a  well  as  being  more  stable  than  any  other 
field;   all  others  are  eliminated  by  lack  of  containment 
because  of  the  specialization  In  the  distribution  function 
To  Invert  the  problem,  establishing  containment  first  and 
then  examining  stability  Is  more  to  the  point  (and  more 
difficult).   The  Taylor  and  Hastle  sufficient  condition 
(1^.3^)  for  the  low  p  shallow  well  expansion  Is  much  more 
relevant  with  regard  to  the  value  of  magnetic  wells,  pro- 
vided that  Its  undue  pessimism  Is  not  selective.   But  the 
analysis  of  scaled  distributions  and  the  general  Inter- 
change sufficient  condition  (14.31)  seem  to  offer  a  much 
more  balanced  view,  especially  since  the  scaled  distribu- 
tions Investigate  the  neighborhood  of  the  borderline 
distribution  f  [J ,\i)  .      We  do  find  that  there  are,  roughly 
speaking,  fewer  restrictions  on  the  distribution  function 
In  a  magnetic  well,  but  It  Is  neither  true  that  a  magnetic 
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well  is  always  stable  nor  that  an  axially  symmetric 
geometry  cannot  be  stabilized  without  cusped  fields 
(both  are  opinions  which  seem  to  be  very  commonly 
accepted).   In  other  words,  modest  advantages  related 
to  cost,  ease  of  construction,  diagnostic  benefits,  etc. 
may  well  dominate. 

The  belief  in  magnetic  wells  as  a  panacea  is  based 
partly  on  the  available  evidence  and  partly  on  a  deep 
psychological  need  for  simplicity.   Certainly  a  result 
which  can  be  simply  stated  and  easily  understood  is 
greatly  to  be  preferred  to  one  whish  is  more  complicated, 
provided  only  that  it  is  not  dangerously  oversimplified. 
One  strong  psychological  support  comes  from  the  compari- 
son of  a  plasma  with  a  diamagnetic  body.    There  is  one 
exact  piece  of  mathematical  evidence  to  bolster  this 
analogue.   If  the  magnetic  field  is  entirely  excluded  from 
the  plasma  and  there  is  a  sharp  separation  surface  between 
the  plasma  and  the  containing  vacuum  magnetic  field,  then 
a  necessary  and  sufficient  condition  for  absolute  stability 
(even  in  the  large)  is  that  the  plasma  be  placed  at  an 
absolute  minimum  of  the  field  [^5]-   The  significance 
of  a  relative  minimum  at  a  finite  non-zero  field  value 
to  preserve  adiabatic  orbits  was  early  recognized  as  a 
qualitative  desideratum  [^6].   The  important  fact  that 
stability  can  be  achieved  by  the  GC  theory  as  well  as  by 
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the  MH  theory  Is  a  crucial  recent  discovery  by  Taylor  [20]. 
But  the  GCP  Is  much  too  complex  a  creature,  with  a  distri- 
bution as  well  as  a  field,  to  be  Interpreted  as  a  classical 
dlamagnetlc  solid.   Although  the  final  answer  Is  far  from 
known,  we  can  already  see  that  the  dlamagnetlc  analogue  Is 
a  qualitative  tendency,  not  a  good  answer. 

The  entire  question  of  Interchange  stability  evapor- 
ates with  tied  ends.   Even  the  question  of  injecting  a  proper 
equilibrium  is  enormously  simpler.   Any  Injection  which  is 
in  local  equilibrium  on  a  line  [i.e.,  f  (£,|-l,A)]  is  (subject 
to  existence  proofs)  an  acceptable  equilibrium.   With  ends 
free,  a  general  f  {e,\i,l\)   will  set  up  large  scale  motions 
(in  particular  Alfven  waves)  to  try  to  achieve  the  natural 
boundary  conditions.   With  ends  tied,  local  stability, 
together  with  mild  limitations  on  the  distribution  and  dis- 
tortion of  field  lines,  should  yield  stability.   For  example, 
the  expansion  introduced  by  Taylor  and  Hastle  yields,  for 
tied  ends,  the  conclusion  that  every  plasma  which  is  locally 
stable  is  absolutely  stable!    From  a  theoretical  viewpoint 
tied  ends  are  so  much  superior  that,  even  with  admitted 


■X- 

We  do  not  claim  this  as  a  proved  theorem;  nor  is 
Taylor  and  Hastle 's  result  more  than  formal.   One  explicit 
omission  in  the  simple  low  p  expansion  (which  leads  to  the 
"proof"  that  interchange  stability  suffices  at  low  p)  is 
remedied  in  this  singular  expansion.   But  it  is  not  known 
what  other  singular  behavior  may  appear  at  low  p. 
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technical  difficulties,  this  would  appear  to  be  as  promising 
a  pursuit  as  complex  field  configurations. 


Interchange  stability  has  been  studied  intensively: 
by  Taylor  in  [20]  for  the  special  distribution  f  (e^'^l)  and 
in  [21]  with  the  more  general  f-.{e,\i,j),    and  in  a  very 
large  number  of  Interesting  forms  and  approximations  by 
Andreolettl  [23] -[28].   In  all  of  these  Interchange  stability 
is  interpreted  as  absolute  stability  at  low  3,  Ignoring  local 
stability.   The  local  criterion  Bf  /he   <  0  is  implicit  when 
the  Kruskal  and  Oberman  5W  is  used;  but  it  Is  omitted  when 
the  Interchange  analysis  is  done  directly.   The  additional 
single-fluid  local  criteria  (7-1)  and  (7-2)  are  found  by 
Taylor  and  Hastie  [22]  in  the  special  forms  appropriate  to 
the  assumption  f  (e,M-)  on  the  one  hand  and  to  the  low  p 
almost  uniform  field  expansion  on  the  other.   The  general 
necessary  condition  in  two  fluids  and  for  an  arbitrary  con- 
tained equilibrium  as  well  as  its  recognition  as  a  micro- 
instability  condition  with  strong  dependence  on  the  shape 
of  f°  are  given  in  [29].   The  combination  of  local  and 
Interchange  stability,  in  particular  for  scaled  distributions, 
is  analysed  in  [30],  where  one  finds  somewhat  less  advantage 
to  magnetic  wells  and  many  possibilities  of  stability  in 
other  geometries,  based  on  what  seems  to  be  approximately 
the  same  available  evidence  as  in  the  other  literature. 
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